Superfluidity and phase transitions in a resonant Bose gas 
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The atomic Bose gas is studied across a Feshbach resonance, mapping out its phase diagram, 
and computing its thermodynamics and excitation spectra. It is shown that such a degenerate gas 
admits two distinct atomic and molecular superfluid phases, with the latter distinguished by the 
absence of atomic off-diagonal long-range order, gapped atomic excitations, and deconfined atomic 
TT-vortices. The properties of the molecular superfluid are explored, and it is shown that across a 
Feshbach resonance it undergoes a quantum Ising transition to the atomic superfluid, where both 
atoms and molecules are condensed. In addition to its distinct thermodynamic signatures and 
deconflned half- vortices, in a trap a molecular superfluid should be identiflable by the absence of an 
atomic condensate peak and the presence of a molecular one. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 
A. Background 

Remarkable experimental advances in manipulating 
degenerate atomic gases have opened a new era in stud- 
ies of highly coherent, interacting quantum many-body 
systems. One of the most striking advances is the 
ability to finely control atomic two-body interactions 
by tuning with a magnetic field the energy (detuning) 
of the molecular Feshbach resonance (FR) through the 
atomic continuumjiiS, This technique has led to a realiza- 
tion of a long-sought-after s-wave paired superfluidity in 
bosonio^i^ and fermionic atomic gases i^ii For fermionic 
atoms, it also allowed the system to be tuned between 
the BCS^ regime of weakly-paired, strongly overlapping 
Cooper pairs (famihar from sohd-state superconductors), 
and the EEC regime of tightly bound, weakly-interacting 
Bose-condensed diatomic molecules. 

Although this crossover has received considerable 
attention)^ i^°d^d^d^d^d^ because of the absence of qual- 
itative differences between the ECS and EEC s-wave 
paired fermionic superfluids, their equilibrium proper- 
ties are already qualitatively well described by early sem- 
inal worksi ^^d'^'d^ In fact for a narrow FR (unfortunately 
not realized by most current experimental systems), the 
crossover can even be computed quantitatively, as a per- 
turbation series in the ratio of the FR width to the Fermi 
energyji^ii^ In such narrow FR systems the crossover 
to EEC takes place when the FR detuning v (quasi- 
molecule's rest energy) ranges from twice the Fermi en- 
ergy 2ep (when it first becomes favorable to convert 
a finite fraction of the Fermi-sea into molecules stabi- 
lized by Pauli-blocking) down to zero energy, where all 
the fermions have become bound into Eose-condensed 
diatomic molecules. The complementary broad reso- 
nance regime of most experiments)^ particularly near 
a universal unitary point^ has been successfully stud- 



ied using quantum Monte Carlo ^^'^^i^^ and field theo- 
retic e-expansion2ii^ and l/iV-expansio n^^i^^ methods 
borrowed from critical phenomena. 

As was recently pointed out— 12^ and is the subject 
of this paper, the phenomenology of resonantly inter- 
acting degenerate hosonic atoms contrasts strongly and 
qualitatively with this picture J2 For a large positive de- 
tuning, molecules are strongly energetically suppressed 
and unpaired atoms (as in any bosonic system at zero 
temperature) form an atomic superfluid (ASF), exhibit- 
ing atomic off-diagonal long-range order (ODLRO)i^ In 
the opposite extreme of a large negative detuning, free 
atoms are strongly disfavored (gapped), pairing up into 
stable bosonic molecules, that then, at T = 0, form a 
diatomic molecular superfluid characterized by a molec- 
ular ODLRO. The MSF does not exhibit atomic ODLRO, 
nor the associated atomic superfluidity. Together with a 
gapped atomic excitation spectrum and correlation func- 
tions (characteristics that extend to finite temperature), 
these features qualitatively distinguish it from the ASF. 

In a trapped, dilute atomic gas the existence of these 
two qualitatively distinct superfluid phases should be 
most directly detectable through independent images of 
atomic and molecular density profiles. As illustrated in 
Fig. [Ija) , the atomic component should exhibit a EEC 
peak in the ASF phase, that is absent in the MSF phase, 
shown in Fig. [IJb). Both superfluid phases are distin- 
guished from the normal state by the EEC peak in the 
molecular density profile, as illustrated in the insets to 
these figures. 

Because of its paired nature, a complementary dis- 
tinguishing characteristic of a MSF are deconfined tt— 
(half-) vortices, topological defects that, in contrast, are 
linearly confined in the ASF state. Consequently, as 
illustrated in Fig. ^ a thermodynamically sharp quan- 
tum phase transition, at an intermediate critical Fesh- 
bach resonance detuning Vc, must separate the MSF and 
ASF phases. Each in turn is also separated by a finite- 
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FIG. 1: Atomic density profiles, ni(r) in (a) the ASF and 
(b) tlie MSF phases. These are distinguished by the presence 
and absence of atomic BEG peak, respectively. Each of these 
superfluid phases is distinguished from the "normal" (ther- 
mal) state by the BEG peak in the molecular density profile, 
n2(r), illustrated in insets. In the dilute limit, the width roo- 
((T = 1,2) of the BEG peak (set by the single-particle Gaus- 
sian ground state wavefunction) , and the extent rxa of the 
thermal part of the atomic cloud, are given by Eqs. (|5.33|) 
and (|5.34p . respectively. 



temperature transition from the "normal" (N) state lack- 
ing any order (i.e., breaking no symmetries). 

Experimental observations of these and associated pre- 
dictions have so far been precluded by a short lifetime 
of the vibrationally hot molecular statei^ The latter 
is believed to be limited by 3-body recombination and 
strongly enhanced atom-molecule scattering near the res- 
onance. In contrast to Fermi systems^''^'^ where Pauli 
exclusion greatly extends the molecular lifetime for a 
positive scattering length and stabilizes the Fermi-sea 
for negative scattering lengths by suppressing multi-body 
collisionS)^^ the resonantly interacting bosonic atomic gas 
is observed to be highly unstable in the negative two- 
body scattering length regimei^^i^i Viable proposals for 
surmounting these problems are currently being investi- 
gated. These include use of an adiabatic ramp of the 
detuning through resonance^^^ or a two-photon Raman 
transition to transfer the Feshbach molecular states to a 
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FIG. 2: The phase diagram (at fixed total density, n) for a uni- 
form condensate as a function of Feshbach resonance detuning 
u and temperature T. A curve of critical detuning Vc{n,T) 
separates the atomic (ASF) and molecular (MSF) superfluid 
phases by a phase transition, which is contirmous between the 
(T = 0) quantum critical point fdn, 0) and a tricritical point 
TCi. The section of the critical curve (gray) between the two 
tricritical points TCi and TC2 denotes a first order transi- 
tion boundary, that terminates the continuous MSF-N phase 
boundary at a critical end point Tco, where three phases meet. 
The critical temperatures T^, a — 1,2, correspond to the far 
detuned limits, v/kBT — > ±cx3. The dashed curve inside the 
ASF phase corresponds to a crossover line, Px{T), at which 
the molecules would condense on their own if there were no 
Feshbach resonance coupling them to the atoms. 



lower lying vibrational statci^ 

Although no direct evidence for an equilibrium Bose 
molecular condensate exists, observed resonant atomic 
loss in a stimulated Raman transition in *^Rb (Ref. [371 ) 
and time domain density oscillations in *^Rb (Ref. ^ are 
consistent with a coherent transfer of population from 
free bosonic atoms to diatomic molecules <^i22. It is not 
at the moment clear (at least to the present authors) 
whether current experimental difhculties of stabilizing 
bosonic atom-molecule mixtures near a FR are funda- 
mental or technical and system specific. One possible 
fundamental source of instability in bosonic atom sys- 
tems is the existence of Efimov bound states of bosonic 
atom tripletsi^i^ At least at a theoretical model level 
these can be suppressed by a sufficiently strong three- 
body repulsion. Even in the unfavorable scenario, where 
such phases of bosons near a FR are indeed metastable 
(as most states of degenerate atoms ultimately are) one 
expects that ideas discussed here should be important on 
sufficiently short time scales and for understanding of the 
associated nonequilibrium dynamics. 

The rest of the paper is organized as follows. The 
Introduction is concluded with a summary of the main 
results and their experimental implications. In Section 
im a microscopic two-channel model, that is believe to 



accurately describe resonantly-interacting atomic bose 
gas is introduced. The model is first used to compute 
the two-body s-wave scattering, showing that it cor- 
rectly captures the Feshbach resonance phenomenology. 
Matching the computed scattering amplitude to its mea- 
sured counterpart allows one to relate parameters ap- 
pearing in the Hamiltonian to experimental observables. 
In Section IIIII a general symmetry-based discussion of 
the expected phases and associated phase transitions in 
this system is presented. In Section IIVI by minimizing 
the corresponding imaginary-time coherent state action, 
the generic mean field phase diagram for the system is 
mapped out. In Sections FVl and FvH this Landau analy- 
sis is supplemented by detailed microscopic calculations 
of phase boundaries, spectra, condensate depletion and 
superfluid density for a dilute, weakly-interacting gas. 
The asymptotic nature of the ASF-MSF phase transi- 
tion is discussed in Section IVIII In Section IVIIII the 
mean field and perturbative analyses, performed within a 
two-channel model, are supplemented with a variational 
theory of a one-channel model. The latter is a better 
description of bosons in which the paired state is absent 
(i.e., there is no long lived metastable paired state with 
distinct internal quantum numbers) once the two-body 
attraction becomes too weak to bind atom pairs (which 
includes, of course, the more familiar regime of two-body 
repulsion). In Section llXl topological defects, vortices and 
domain walls, in the ASF are studied, and the ASF-MSF 
and SF-to-normal fluid transitions are characterized in 
terms of a proliferation of these topological defects. The 
paper is concluded in Section 1x1 



B. Summary of results 

In this paper a considerable elaboration and extension 
of predictions reported in a recent Letter— are presented. 
The primary results are summarized by the density pro- 
files in Fig. [1] and the phase diagrams in Figs. [5] and 
[TOl characterizing the phases and phase transitions of 
a resonant Bose gas. As illustrated there, it is found 
that Feshbach-resonantly interacting atomic Bose gas, in 
addition to the normal state exhibits two distinct low- 
temperature superfluid states. The first, appearing at 
positive detuning, is the more conventional atomic su- 
perfiuid, characterized by coexisting atomic and molec- 
ular BEC and their associated ODLROs, with finite or- 
der parameters \E'io and ^^20, respectively. The other, 
more exotic, MSF state, appearing at low temperature 
and negative detuning, is characterized by superfluidity 
of diatomic molecules, with a finite molecular condensate 
order parameter ^20- It is distinguished from the ASF 
by the absence of atomic ODLRO, i.e., inside the MSF 
phase v&io = 0. 

As illustrated in detail in Sec. lIVBl a finite ^^lo always 
implies a finite Vl'2o- In the presence of an atomic con- 
densate, V&io, the Feshbach resonance coupling allows a 
scattering of two Bose-condensed atoms out of the atomic 



BEC into the molecular BEC (i.e., ASF is really a su- 
perposition of Bose-condensed open-channel atoms and 
Bose-condensed closed-channel molecules) and therefore 
acts like an ordering "field" on the molecular order pa- 
rameter. This implies that a state in which atoms are 
condensed but molecules are not is forbidden by general 
symmetry principles!^ 

As noted above, a vivid signature of two distinct su- 
perfiuid orders should be detectable via time-of-fiight 
shadow images. In the dilute regime (described by a 
BEC approximation), the resulting images are schemat- 
ically illustrated in Fig. [T] At higher densities, where a 
local density approximation is more appropriate, it is ex- 
pected that for a range of atom number and detuning, 
phase boundaries as a function of chemical potential in 
the bulk system (see e.g., Fig.[TU]) will translate into shell 
structure^^iiiiiS which should also be observable experi- 
mentally in time-of-fiight shadow images. 

As for any neutral superfiuid, ASF and MSF are each 
characterized by an acoustic (Bogoliubov) "sound" mode, 
illustrated in Fig. [31 corresponding to long wavelength 
condensate phase fiuctuations, with long wavelength dis- 
persions 
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(1.1) 



where c^ (with a — ASF or MSF, or equivalently 1 or 
2) are the associated sound speeds with cmsf given by 
(|6.29p and casf given by (|6.46p in terms of the interaction 
parameters of the model (see Sec lVIBJ) . 

In the ASF state the gapless mode corresponds to in- 
phase fluctuations of the atomic and molecular conden- 
sates. However, as illustrated in Fig. [31 in contrast to 
ordinary supcrfluids, ASF and MSF also each exhibit a 
gapped branch of excitations. 
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with gaps Sfap gjygn explicitly by ([Oil) and ([QUI) , 
while the quadratic corrections b^ may be inferred from 
the general forms (|6.8p and (|6.2ip of the spectra. In the 
ASF the gap is controlled by out-of-phase fluctuations of 
the atomic and molecular condensates, and is set by the 
Feshbach coupling a. 

In the MSF, gapped excitations are single atom- 
like quasiparticles akin to Bogoliubov excitations in the 
paired BCS state, that however do not carry a defi- 
nite atom number. These single-particle excitations are 
"squeezed" by the presence of the molecular condensate, 
offering a mechanism to realize atomic squeezed states.^ 
We expect that these gapped, atomic quantum fluctua- 
tions associated with the presence of the molecular con- 
densate can be measured by interference experiments, 
similar to those reported in Ref. i47:. As detailed below, 
the low-energy nature of these excitations is guaranteed 
by the vanishing of the gap at the MSF-ASF transition, 
Vc, with £^m|f(j^c) = 0. 

In the dilute weakly interacting limit appropriate to 
atomic gases, the ASF-N and MSF-N transition tem- 




FIG. 3: Schematic low-energy excitation spectra character- 
izing ASF and MSF phases. In the ASF the acoustic and 
gapped branches of excitations correspond to in-phase and 
out-of-phase fluctuations of atomic and molecular conden- 
sates, respectively. In the MSF state the acoustic branch is 
the standard gapless Bogoliubov mode. The gapped branch 
corresponds to atom-like quasiparticle excitations that are 
squeezed by the presence of the molecular condensate. At 
the critical detuning Uc, the gap closes, signaling a quantum 
MSF-ASF phase transition. 



peratures, Tci{v) and Tc2{v), respectively, for a three- 
dimensional (3d) bulk uniform system are well approxi- 
mated by 
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2-5/3c2/3 andc 



(1.3) 
29/2^i/2/3cC(3/2), a2 = ai/8, 6i = c^/^, 62 = 
l-)-25/2 One sees that T^i > Tc2, with 
the asymptotic ratio 61/62 = 2^^^ set by the mass and 
boson number that both differ by a factor of two between 
the two phases. When interactions are included, for ly ^ 
the asymptotic nature of these thermal transitions is in 
the well-studied classical 3d XY universality class. 

As illustrated in Fig. [5] (see also Figs. [TUl and [TT|) . in the 
vicinity of the critical endpoint i> ~ 0, T ^ TcOi where 



J 



the three phases, N, MSF, ASF meet, a coupling of the 
molecular and atomic superfluid order parameters con- 
verts a section (between the tricritical points TCi and 
TC2) of the (otherwise) continuous N-ASF and MSF- 
ASF transitions to first order. ■^'' The resulting crossing 
point at Tco, that terminates a continuous N-MSF tran- 
sition is a critical endpoint (CEP). In the dilute limit, 
the CEP temperature is given by 
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(1.4) 



As illustrated in Sees. HVl and IVlIl the appearance of a 
first order transition (even in mean field theory) near N- 
MSF continuous boundary, is a generic feature resulting 
from the coupling of the ASF order parameter to the 
critical MSF order parameteri^l 

The corresponding transition temperatures in a trap 
(here distinguished from bulk quantities by a tilde) are 
also easily computed and in 3d are given by 
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(1.5) 

with a^ = 2C(2)/9cr\(3), and b„ = 2/ia^C{i)- The 
transition temperatures in the limit of asymptotically 
large positive (cr = 1) and negative {a = 2) detuning 
{\v\/kBTca 3> 1), and at the tricritical point [v = 0), are 
given, respectively, by 
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(1.6) 
(1.7) 



where ojq is the trap frequency. Comparing the first lines 
of (|1.3p and (|1.5p , note that the latter is now approached 
linearly with, rather than as the square-root of the re- 
duced detuning from either side. 

In the dilute limit the thermodynamics is also easily 
worked out. In the 3d bulk system, the condensate den- 
sities for the atomic and molecular BEC are given, re- 
spectively, by 



nw{T,v) = n 
n2o{T,v) = - 
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(1.8) 
(1.9) 



where C{3/2) ~ 2.612 and ga{x) — X^^i^"/"-" i^ the extended zeta function4^ 



As illustrated in Fig. [2] the MSF-ASF transition takes 
place at a critical value of detuning VciT, n) determined 
by the strength of atomic and molecular interactions, 
shifting it away from its noninteracting value of 0. At 
zero temperature this is a continuous quantum phase 
transition that for a d-dimensional system is in the (d+l)- 
dimensional classical Ising universality class^i^i^ with 



i/c(0, n) w -(.92/2 - gi2)n - 2a\/2n, 



(1.10) 



where gi, gi2, and 32 are, respectively, the atom- 
atom, atom-molecule and molecule-molecule interaction 
strengths, related in the standard way to the correspond- 
ing scattering lengthsfi^ and a is the Feshbach resonance 
coupling. The transition at Vc is characterized, upon ap- 
proach from the MSF side, by the vanishing of the single- 
atom excitation gap £'^gp(z/), and, upon approach from 
the ASF side, by the disappearance of the atomic conden- 
sate nio{i'). At zero temperature, in the critical region 
these are predicted to vanish according to 



nio{0,iy) r^ \v- Vc\ 
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<sf(o,^)- k-^cr^ (1.11) 



where /3/ and vi are, respectively, the order parame- 
ter and correlation length exponents for the (d + 1)- 
dimensional Ising model. One may hope that when long- 
lived molecular condensates are produced, nontrivial be- 
havior of E^'g^iv) and the full excitation spectra, E^{k) 
may be observed in Ramsey fringes^ and in Bragg and 
RF spectroscopy experiments^i^i^i^. 

At finite temperature, away from the critical endpoint 
Tco the transition is in the classical d-dimensional Ising 
universality class. Scaling, together with the relevance 
(in the renormalization group sense) of T at the quantum 
critical point, also implies a universal shape of the low T 
part of the MSF-ASF phase boundary 



I/c(n,T)-l^e(n,0) + a^l/'^^ 



(1.12) 



illustrated in Fig. [2l 

The dashed curve, Vx (T) inside the ASF phase of the 
phase diagram (Figs. [Hand [TO]) denotes a crossover (that 
becomes sharp with a vanishing Feshbach resonance cou- 
pling a) between ASF regimes with low and high values 
of the molecular condensate 7120. In the absence of the 
coupling, the molecules would condense on their own for 
v < Vc{T). For small a the weak symmetry breaking field 
generated by the atomic condensate smears this transi- 
tion into a ASF-AMSF crossover, and leads to small, but 
finite, 77,20 even for v > Vc{T). 

As for any superfluid, the ASF and MSF phases also 
exhibit interaction-driven condensate depletion Stiq = 
n — nio — 2n2o, quantifying the fact that, even at T = 0, 
not all atoms are in the condensate. At T = these 
are computed in Sec. IVIDI An interesting feature, illus- 
trated in Fig. [31 is that dnQ{i>) exhibit a cusp maximum 
at Vc, 
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FIG. 4: Schematic of the zero-temperature depletion 5no(0, v) 
as a function of detuning. The far detuned limits are given by 
Suq'^ = o-(8/3\/3)(ffl<T7i(To)'^ 5 where >, < correspond to ct = 
1, 2 respectively. The cusp peak is given by (|1.14|) . displaying 
the power- law form \i^ — Vc]'', with p = 1 in d = 3. 



associated with enhanced role of quantum fluctuations 
at the MSF-ASF transition. The maximum depletion is 
given by 
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(1.14) 
where 02 is the molecule-molecule s-wave scattering 
length. Outside the critical region one expects p — 1, 
crossing over to p = 1 — ai inside it, where ai is the 
{d + l)-dimensional Ising specific heat exponent. 

Another important qualitative distinction between the 
ASF and MSF phases is the nature of their topologi- 
cal excitations, namely vortices. The paired nature of 
the MSF allows for tt— (half-) vortices (as in a BCS su- 
perconductor), while the ASF, being (from a symmetry 
point of view) a standard "charge-one" superfluid, admits 
only standard 27r-vortices. However, pairing correlations 
present in the ASF lead to an interesting 7r-vortex exper- 
imental signature even in the atomic superfluid. Thus, 
in the ASF phase a seemingly standard 27r-vortex in an 
atomic condensate, ^10 > 0, will generically split into 
two TT- vortices (see Fig. [5]) confined by a domain wall of 
length 
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(5^0 (^) = Sno{vc) - c„\v - Vc\^- 



(1.13) 



that diverges as the ASF-MSF phase boundary is ap- 
proached from the ASF side. Sufficiently close to the 
transition, it is expected to track the associated correla- 
tion length. 

This confinement arises because in the large Feshbach 
coupling limit a 27r-vortex in the atomic condensate in- 
duces a 47r-vortex in the molecular condensate. Such a 
double molecular vortex is unstable to two fundamental 
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FIG. 5: 27r atomic condensate vortex in the ASF splits into a 
TT+TT vortex pair connected by a "normal" domain wall, whose 
length _Ro increases as the FR coupling a becomes weaker. 



27r molecular vortices that, in 2d, repel logarithmically, 
but are confined linearly inside the ASF phase. This 
provides a complementary formulation of the ASF-MSF 
transition as a confinement-deconfinement transition of 
TT (half) atomic vortices. 



experiments^. A resonant interaction is the key fea- 
ture special to a select class of atomic (fermionio^i and 
bosonio^) systems. A fully microscopic description of 
such resonant interactions is quite complex, involving a 
full set of internal nuclear and electronic spin degrees of 
freedom characterized by hyperfine states, mixed upon 
scattering by the interatomic exchange interaction. How- 
ever, in the vicinity of a resonance, two-atom scatter- 
ing in the "open" channel is dominated by hybridization 
with a two-atom molecular bound state in the "closed" 
channel, thereby allowing one to neglect all other off- 
resonant channels. As illustrated in Fig. [SI the two chan- 
nels are distinguished by the two-atom electron spins, 
with the open-channel an approximate spin-triplet and 
closed-channel an approximate singlet.— Consequently 
they have different Zeeman energies, allowing the cen- 
ter of mass rest energy v of the closed-channel molecule 
(bound state) to be tuned, relative to the open-channel 
two-atom continuum, via an external magnetic field. 
This yields an unprecedented tunability of the effective 
atomic interaction strength by varying a magnetic field. 



II. TWO-CHANNEL FESHBACH RESONANT 
MODEL 

The goal here is to model a resonantly interact- 
ing atomic Bose gas of the type realized in recent 



The two-channel model describing the resonant atom- 
molecule system is characterized by the following grand- 
canonical Hamiltonian: 
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where V'J. (r) , tAo- (r) are bosonic creation and annihilation 
field operators for atoms (tr = 1) and molecules (cr = 2). 
They are described by respective single-particle Hamil- 
tonians 



h„ 



Ira, 



-V2+Ai^ + T4(r), 



(2.2) 



with atomic and molecular masses m\ — m and 777,2 — 2777 
and effective chemical potentials fii = fi and /i2 — 2/x— t^q- 
The (bare) detuning parameter vq is related to the energy 
of a (closed-channel) molecule at rest, that can be exper- 
imentally controlled with an external magnetic field. In 
the ensemble of fixed total number of atoms N (free and 
bound into molecules), relevant to trapped atomic gas 
experiments, the chemical potential /i is determined by 
the total atom number equation 



N ^ dv 



(V31(r)^i(r))+2(^^(r)V;2(r)) 



(2.3) 



The positive local pseudo-potential parameters gi, 52: 5i2 
measure background (nonresonant) repulsive atom-atom. 



atom-molecule and molecule-molecule interactions, re- 
spectively, and in the dilute limit are proportional to cor- 
responding background 2-body s-wave scattering lengths. 
The Feshbach resonance coupling a characterizes the co- 
herent atom-molecule interconversion rate (hyperfine in- 
teraction driven hybridization between open and closed 
channels), encoding the fact that a molecule can de- 
cay into two open-channel atoms^S*^. The external po- 
tentials Vcr(r) describe the atomic and molecular traps, 
which for most of the paper will be taken to be a "box" , 
modeled (for convenience) using periodic boundary con- 
ditions. 

In principle it is possible to obtain the above Hamilto- 
nian H from a microscopic analysis of atoms interacting 
via a Feshbach resonance. ^-^^ However, its validity is ul- 
timately justified by the fact that for two atoms in a 
vacuum (for which one takes /i = 0) it reproduces the ex- 
perimentally observed Feshbach-resonance phenomenol- 
ogy. Namely, it predicts an atomic scattering resonance 
for positive detuning, a true molecular bound state for 
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FIG. 6; Above: Schematic illustration of a Feshbach reso- 
nance, modeled by two coupled channel interaction potentials 
(distinguished by two-atom electronic spin states) as a func- 
tion of inter-particle separation. The (so-called) "open" chan- 
nel is too shallow to support a bound state, while the other, 
"closed" channel supports a bound state or a resonance (in- 
dicated by a dashed line) that is tuned with a magnetic field 
via the Zeeman splitting between the two channels. Below: 
At resonance z/ = 0, when a bound state first appears, the 
s-wave atomic scattering length diverges according to (|2.4p or 



negative detuning (illustrated in Fig. [7]), and an s-wave 
scattering length of the experimentally observed form 



flbg 1 



B,, 



(2.4) 



B-Bo, 
Here, Obg is the background (nonresonant) scattering 



length, i?u, is the experimental width (not to be confused 
with the width of the Feshbach resonancei^) , and Bq is 
the value of the magnetic field at which the Feshbach 
resonance is tuned to zero energy. 

These properties follow directly from the s-wave atomic 
scattering amplitude fo{E) that for two atoms in a vac- 
uum can be computed exactlyJ^ii^i^ Focusing for sim- 
plicity on the resonant part of the interaction, (i.e., tak- 
ing go- = gi2 = 0) the scattering amplitude is given by 



ME) 



/n 



m E 



ivnE 



(2.5) 



with V the renormalized (physical) detuning and Fq a 
parameter measuring the width of the resonance. These 
are given by 



t'o 



h^ 



(Pp 1 



a^m^ 



" ~ 167r2?i6- 
The latter is related to an effective range parameter 



(2.6) 

(2.7) 



ro 



= -2h/^mro. 



(2.8) 



The integral in (|2.6p is implicitly cut off by the ultravi- 
olet (uv) scale A « 2TT/do, set by the inverse of the size 
do of the closed-channel (molecular) bound state [below 
which the point interaction approximation inherent in 
the Hamiltonian (|2.ip breaks down], so that 



1^0 






(2.9) 



relates the bare and physical detuning. 

The s-wave scattering length a^ = — /o(0) is then given 

by 



Fo^i 
m V 



(2.10) 



Thus, to reproduce the experimentally observed scatter- 
ing length variation with magnetic field (|2.4p . one fixes 
the detuning to be i/ « 2^j,b{B — Bq), with the approx- 
imate Bohr magneton proportionality constant set by 
the Zeeman energy difference between approximate elec- 
tronic spin-triplet (open) and spin-singlet (closed) chan- 
nels. Matching (|2.10p to (|2.4p also allows one to relate the 
Feshbach resonance coupling to the "width" B^, giving 



Fo 



Amnlal^Bl/h^. 



(2.11) 



Interpretation of the scattering physics in terms of an 
intermediate molecular bound or quasi-bound state fol- 
lows from the poles of fo{E), together with appropri- 
ate constraints arising from boundary conditions on the 
molecular wavefunction. From (j2.5[) the physical pole is 
given by 



E„ 



Er - iV/2, 



(2.12) 
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FIG. 7: The real part of the pole of the scattering ampli- 
tude fo{E), (I2.12|l . as a function of detuning, parameterized 
here by —l/us = {v/h)^/rn/ro, with — 1/ro = \/mro/(2ft). 
As discussed in the text, bound states and resonances must 
correspond to physical solutions of the Schrodinger equation 
with proper boundary conditions. The thin dotted line indi- 
cates asymptotic linear behavior of the bound state for small 
positive as- 



where E,. ^ v - Tq/2 and T = Va{Av/TQ ~ Xfl"^. For 
negative detuning, z^ < 0, the pole is purely real and 
negative, corresponding to a bound state with energy 



E- 



-z/VFo, |i^|«ro, 

V, \v\ > To. 



(2.13) 



and infinite lifetime. The corresponding wavevector fc~ — 
iy^2m\Ep\ lies on the positive imaginary axis, and cor- 
responds, as required, to a wavefunction ~ e~' p'^ that 
decays exponentially at infinity. As z^ — > 0+ one has 
Ep — > 0~, and the bound state coincides with the bot- 
tom of the continuum. 

For v > 0, one would like to be able to interpret the 
state created by V'2 as a metastable molecule with a finite 
decay time into two atoms. Such an interpretation makes 
sense only if the real and imaginary parts of Ep are posi- 
tive and negative, respectively, and Reii'p > |Im£^p| — Ep 
is called a resonance in this case, with the inequality be- 
ing the condition that a well defined resonance have a 
width that is narrower than its energy. However, as seen 
in Fig. [71 for a range of small positive v, rather than 
moving to positive values, Ep remains real and moves 
back to negative values. But this does not indicate a re- 
stored bound state because fc+ = —iy^2in\Ep\ now lies on 
the negative imaginary axis and corresponds to a wave- 
function ~ e"'"'''? '*" that grows at infinity. Thus, in this 
range of detuning the pole no longer corresponds to a true 
bound state and is often referred to as a virtual bound 
state<S 



Although for ro/4 < u < To/2, one has ImEp < 0, 
indicating a finite decay time, the real part of Ep remains 
negative. Only for i^ > To/2 is KeEp > 0, and only for 
v significantly larger than ro/2 does one obtain a true 
resonance, with a restored molecular interpretation for 
the ■02 field. 

This behavior is summarized in Fig. [7l where HeEp 
is plotted as a function of detuning. It should be em- 
phasized that the issue here is strictly one of physical 
interpretation of the microscopic scattering states. The 
model remains well defined, and is a valid description 
of experiments in the FR regime, over the full parameter 
rangci^. The thermodynamic phases that will be derived 
in later sections are also well defined for all parameters. 

In addition to a gas parameters n^'^a]^', correspond- 
ing to background scattering lengths associated with cou- 
plings gi, g2, and gi2 (that are constant in the neighbor- 
hood of the Feshbach resonance) , the two-channel model 
(|2.ip is characterized by a dimensionless parameter 



To 



ksTBEC 



TT 72 



1/3 



'^ol 



(2.14) 



that measures the effective strength of the Feshbach- 
resonant interaction relative to the kinetic energy 
ksTBEC — {2TTh^ /ra)'n?/'^ set by the total atomic density 
n. As long as these are all small, i.e., the gas is dilute 
with respect to background scattering lengths and the 
Feshbach resonance is narrow, the description of phases 
(i.e., properties away from any phase transitions) can be 
accurately given by a perturbative expansion in these di- 
mensionless interaction parameters . ^^'^^ i^ This will be 
verified through explicit calculations in Sec. IVIDI Phys- 
ically, this narrow resonance limit, 7 <SC 1 corresponds to 
molecules that are predominantly in the closed-channel, 
i.e., have a long lifetime before decaying into two free, 
open-channel atoms. 

In the opposite, broad resonance (7 S> 1) limit the 
molecular wave function is strongly hybridized with the 
open-channel, and is characterized by a high density 
of continuum states above threshold. For such a sys- 
tem, although for negative detuning a bound molecu- 
lar state still exists, no resonance remains for positive 
detuning — the physics of this regime can no longer be 
interpreted in terms of populations of coexisting atoms 
and (metastable) molecules. In this limit the disper- 
sion of the bare molecular field can be neglected and 
"02 can be adiabatically elinrinated (integrated out, ig- 
noring the subdominant atom-molecule and molecule- 
molecule density interactions) in favor of two open- 
channel bosons. ^-'^^ The resulting single-channel broad 
resonance model Hamiltonian is then given by 



Ih 



1 



Ih 



Hi_,h ^ I dY\ ^l{r)hi,jji{r) + ^.g0{(r)20i(r)2 
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(2.15) 



where g is the effective atom-atom interaction couphng, 
approximately given hy g = gi — o? jv, and a stabiliz- 
ing (against collapse) three-body interaction has been 
added, with a coupling w > 0. In contrast to the 7^1 
two-channel model, due to the divergence of 5, the one- 
channel model is strongly interacting when the resonance 
is tuned to low energy, 1/ « 0, and the s-wave scattering 
length exceeds the inter-particle spacing, i.e., nc?^ ^ 1. 
Consequently, in this regime predictions derived from a 
perturbative analysis of the one-channel model can only 
be qualitatively trustworthy. Moreover, as \g\ and w in- 
crease, quantitative predictive power would require one 
to include higher than three-body interactions, so (|2.15[) 
really only makes sense when \g\\s not too large. 

Considered as a function of [i\ and g. the single channel 
model (|2.15p also exhibits the same three N, MSF and 
ASF phases (with MSF requiring g < 0). So long as |g| 
is not too large (i.e., so long as the associated scattering 
length obeys njagp <C 1), a variational BCS approach 
can be used to accurately compute the thermodynamics. 
This approach is discussed in Sec. IVIIII 

With the model defined by i^,^° Eq. dm, the ther- 
modynamics as a function of a chemical potential [i (or 
equivalently atom density, n), detuning v and tempera- 
ture T can be worked out in a standard way by computing 
the partition function Z = Tr[e"'^^] (/3 = l/ksT) and 
the corresponding free energy F = —kBTlnZ. The trace 
over quantum mechanical states can be conveniently re- 
formulated in terms of an imaginary-time (t) functional 
integral over coherent-state atomic (cr = 1) and molecu- 
lar {a — 2) fields ipa{i^, t), 



Z 



D^aDfpa e 



-s/h 



where the imaginary-time action is given by2i 



S^ 



I3h 



dr / dr 



Y^rMr^. + H{r,,rl:,) 



(2.16) 



(2.17) 



The total atom number constraint (|2.3p that allows one 
to eliminate fi in favor of N is then simply given by 






(2.18) 



III. 



SYMMETRIES, PHASES AND PHASE 
TRANSITIONS 



Before delving into detailed calculations it is instruc- 
tive to first consider the general symmetry-based charac- 
terization of phases and transitions between them. Since 
both atoms (V'l) and molecules {1P2) are bosonic and 
can therefore Bose condense, the system's thermodynam- 
ics is determined by two Bose-condensate order param- 
eters, ^10 and ^20, respectively. As usual, microscopi- 
cally, these label thermodynamic averages of the corre- 
sponding field operators, or, for weakly interacting sys- 
tem, equivalently, are single-particle wavefunctions into 



which all bosons (atoms and molecules, respectively) 
Bose-condense. 

The condensate fields ^10, ^20 are legitimate order 
parameters that uniquely characterize the nature of the 
possible phases. Naively one would expect four phases: 
(i) Normal (N) (^10 = ^20 = 0), (ii) (*io ^ 0, ^20 = 0), 
(iii) MSF (^10 = 0, ^20 7^ 0), and (iv) AMSF (^10 ^ 0, 
^20 7^ 0), corresponding to four different combinations of 
vanishing and finite order parameters. However, a finite 
Feshbach resonance interaction explicitly breaks C/(l) x 
C/(l) symmetry of the a = Hamiltonian down to C/(l) x 
Z2 . Physically, this corresponds to the fact that only the 
total number of atoms 



N ^ Ni+ 2N2, 



(3.1) 



((^I^i) + 2(V'^V2))V^ 



(V is the system volume) is conserved in the presence 
of Feshbach resonant scattering (break up) of a molecule 
into two atoms, rather than a separately conserved num- 
ber of atoms, Ni and a number of molecules, N2. This 
is why only a single chemical potential /i is introduced in 
H, Eq. (|2.ip . Consequently, a Feshbach resonant inter- 
action requires a condensation of molecules (ordering of 
-02) whenever atoms are Bose-condensed ("01 is ordered) 
and thereby forbids the existence of the state (^10 7^ 0, 
^20 = 0). As a result, the system of resonantly inter- 
acting bosonic atoms exhibits only three distinct phases: 
N, MSF, and AMSF-^ since the atom-only condensate is 
impossible, for brevity of notation the AMSF state will 
often be referred to as simply ASF, using the two names 
interchangeably 

In addition to the fully "disordered" normal state 
that does not break any symmetries, the above dis- 
tinct thermodynamic phases are associated with differ- 
ent ways that the f/(l) x Z2 symmetry is broken. ^^ Bose- 
condensation of molecules breaks the U{1) subgroup and 
corresponds to a N-MSF transition to the MSF phase. 
Since it is characterized by ordering of a complex scalar 
field, "02, this transition is in an extensively-explored 
and well-understood XY-model universality class^. The 
low-energy excitations in the MSF phase are gapless 
Goldstone-mode phase fluctuations of the condensate 
^20j associated with the broken [/(I) symmetry. In ad- 
dition there are gapped excitations associated with the 
magnitude fluctuations of ^20- 

The breaking of the remaining Z2 symmetry is associ- 
ated with Bose-condensation of atoms, tpi, in the pres- 
ence of a molecular condensate. As will be shown explic- 
itly in Sec. IVIIi the corresponding MSF-AMSF transi- 
tion is associated with the ordering of a real scalar field, 
and one would therefore expect the MSF-AMSF transi- 
tion to be in the well-explored Ising universality class. 
However, as will be seen, a coupling of the scalar order 
parameter to the strongly-fluctuating Goldstone mode of 
the MSF phase has a nontrivial effect on the Ising tran- 
sition, quite likely driving it first order sufficiently close 
to the transition i^°i^^ 



10 



Bose-condensation of atoms (ordering of ipi) directly 
from the normal state breaks the full U{1) x Z2 symm- 
metry and corresponds to a direct, continuous N-AMSF 
phase transition. Since it is associated with the ordering 
of a complex scalar field, one expects (and finds) it also 
to be in the well-studied XY- model universality class, 
and to exhibit a single Goldstone mode corresponding 
to common (locked) phase fluctuations of the condensate 
fields ^01 and ^02- 

There is an instructive isomorphism of this descrip- 
tion, in terms of two complex scalar order parameters 
$10 and Vl/20, to that in terms of a two-dimensional 
rank-1 vector (Mq), together with a rank-2, traceless, 
symmetric tensor (Qo) order parameter. The latter de- 
scription is well known in the studies of ferroelectric ne- 
matic liquid crystals. ^^ There, ordering of Qo describes 
the isotropic-nematic transition, where the principal axes 
of mesogenic molecules align macroscopically, breaking 
the 2d rotational symmetry modulo tt rotation. The 
latter remains unbroken in the nematic phase. This is 
isomorphic to the N-MSF transition discussed above. 
For polar molecules, at lower temperature this isotropic- 
nematic transition can be followed by vector ordering (of, 
e.g., the molecular electric dipole moments) of M. In 
the non-rotationally invariant, quadrupolar environment 
of the nematic phase, this corresponds to spontaneous 
breaking of the remaining Z2 symmetry. Therefore, the 
subsequent nematic-polar (ferroelectric) transition corre- 
sponds to a spontaneous selection between two equivalent 
and TT orientations of M relative to the molecular princi- 
pal axes characterized by the nematic Qq order. Clearly, 
this latter transition can be identified with the MSF- 
AMSF transition in the atomic system. The mathemat- 
ical mapping between the two descriptions is elaborated 
on in more detail in Appendix \K\ 



IV. MEAN FIELD THEORY 

The first goal is to determine the nature of the phases 
and corresponding phase transitions exhibited by the res- 
onant bosonic atom-molecule model introduced above. 
To this end one must evaluate the free energy, which in 
the presence of interactions and fluctuations can only be 
carried out perturbatively. However, away from contin- 
uous phase transition boundaries, i.e., well within the 
ordered phases, fluctuations are small. The functional 
integral in Z is then dominated by field configurations 
'0(j(r, t) ~ ^0-0 (r) that minimize the action S*, and there- 
fore can be evaluated via a saddle-point approximation. 
For time-independent solutions that characterize ther- 
modynamic phases, the order parameters ^o-o(r) equiva- 
lently minimize the variational energy functional iJ[^cro], 
in which one substitutes the classical order parameters 
for the field operators in (|2.1[) . For the case of a uniform 
bulk system with Va (r) — (that is the focus of this sec- 
tion) one expects that -ff[^CTo] is minimized by a spatially 
uniform solution (see however Ref. 1651) ^o-o = |^cro|e'^"- 



A mean field analysis then reduces to a minimization of 
the energy density 



Tir, 



i/[*10,*20]/T^ 



(4.1) 



= -ml*ior + y l^ior - M2|'J'2or + y l*2o| 

+5121*10^*201^ - aRe[*;o*?o]. 

Total atom number conservation ensures a global t/(l) 
symmetry with respect to uniform, cr-independent phase 
rotation. The Feshbach resonance interaction 



aRe[*;o*2^] = a|*2o||*ioPcos(20i-e2) 



(4.2) 



clearly locks atomic and molecular phases together, anal- 
ogously to two Josephson-coupled superconductors, with 
the energy minimized by 



62 = 26li(mod27r). 



(4.3) 



where without loss of generality a is taken to be positive; 
for a < 0, the molecular phase 62 is simply shifted by tt. 
The corresponding saddle point equations are given by 



= 



dH 



mf 



5*10 



(4.4a) 



a*io*20 + *io (-Ml + 5i|*io|^ + 5i2|*2 



5H„ 



5*20 



(4.4b) 



-a*?o + *20 (-Ai2 +.92|*20p +312!* 
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For a trapped system with a fixed total number of atoms 
appropriate to atomic gas experiments, these equations 
must be supplemented by the total atom number (|2.3p — 
given by N/V = |*ioP + 2|*2oP within the mean field 
approximation — so as to map out the phase diagram as 
a function of atom number N and detuning v. However, 
it is simpler to instead treat the atomic (^i) and molecu- 
lar (/i2) chemical potentials as independent variables and 
first map out the phase behavior as a function of /ii and 

Ai2- 



A. Vanishing Feshbach resonance coupling, a = Q 

To complete a minimization of Timf it is instructive to 
first consider a special case of vanishing Feshbach reso- 
nance coupling, a = 0, for which saddle-point equations 
reduce to 



= *io(-A*i+.9i|*ioP+5i2|*20P) (4.5a) 
= *2o(-A*2+52|*20P+ffi2|*ioP)- (4.5b) 

In this special case the model corresponds to an easy- 
plane (XY-model) limit of two energetically-coupled fer- 
romagnets, a model that has appeared in a broad variety 
of physical contextsi^i^ 
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In the a = limit, the model exhibits four different 
phases, corresponding to the four different combinations 
of zero or nonzero order parameters ^lo, ^20- 

For /ii < and /i2 < 0, Timf is convex with a unique 
minimum at 



1*1 



1*2 



0, for ^i,/i2 < 0, 



(4.6) 



corresponding to the normal state. 

For /ii > and /i2 < 0, the minimum continuously 
shifts to 

912 

(4.7) 



l*iol = V/^1/51, |*2o| = 0, for fii > 0, ^2 < Ail 



corresponding to a continuous transition at 

Ml = 0, Ai2 < (Normal-ASF transition hne) (4.8) 

from the normal state to the atomic superfluid (ASF), 
where atoms are Bose condensed but molecules are not. 
Substituting this solution into the second saddle-point 
equation, (j4.5bll . one confirms that *20 = is indeed a 
minimum so long as fi2 < ^1312/51- 

In a complementary fashion, for ^2 > 0, *20 becomes 
nonzero, while ^ig continues to vanish so long as fix < 
^2512/32- Hence the normal state undergoes a transition 
to the molecular superfluid (MSF) along the line 

Ml < 0, M2 = (N-MSF transition hne). (4.9) 

The MSF phase is characterized by order parameters 

|*2o| = \/m/92, l^iol = 0, for H2 > 0, Ml < M2 — , 

52 

(4.10) 
i.e., Bose condensed molecules but uncondensed atoms. 
Along the transition line 

^2 = Mi^— , Ml > (ASF-AMSF transition line) 
51 

(4.11) 

the ASF phase becomes unstable to development of a 

nonzero molecular Bose condensate. Conversely along 

the line 

^1 = /i2 — , M2 > (MSF-AMSF transition line) 

52 

(4.12) 
the MSF phase becomes unstable to a finite atomic Bose 
condensate. It is clear that as long as (7^2 < 5i52 the re- 
gion defined by above two boundaries is finite and corre- 
sponds to atomic-molecular superfluid (AMSF) in which 
both atoms and molecules are condensed, with 



*io| = 



/ 52M1 -512M2 
5152 - gh 



l*20l = 



r ^ n 512 ^ , ^ 51 

for Hi, M2 > 0, < M1/M2 < • 

52 512 



/ 51M2 -512M1 

5152 - gh 

(4.13) 



For (7^2 > 5152 this fourth AMSF phase is absent 
because energies —(x^^/2ga of the ASF and MSF min- 
ima cross before either becomes locally unstable. Conse- 
quently, instead of continuous transition to AMSF, the 
system undergoes a first order transition between the 
atomic and molecular superfiuids at 



for gl^ > 5152, (first order ASF-MSF line). 



Within this mean field analysis, for this range of param- 
eters all transitions above are second order. 



Ml _ 

M2 

(4.14) 
where the ASF and MSF minima become degenerate. In 
this case, the lines mi = M2512/52 and M2 = M1512/51 are 
spinodals, beyond which the phase on the opposite side 
of the first order line becomes locally, not just globally, 
unstable. As usual with a first order transition, for a 
fixed total atom number (relevant to trapped atomic gas 
experiments), along this first order line corresponds to a 
coexistence region where the system phase separates into 
coexisting atomic and molecular superfiuids. The two 
possible phase diagrams for a = are illustrated in Fig. 

E 



B. Finite Feshbach resonance coupling, a 7^ 

As is clear from the Hamiltonian (|2.ip . and its mean 
field form (|4.ip . the full system is characterized by a finite 
Feshbach resonance coupling a > 0. If MI1M2 < it is 
easy to see from (|4.4p that the normal phase *io = *20 = 
remains a local minimum of Timf , and a nonzero a does 
not affect the boundaries of the normal phase, with (14. 8p 
and (14. 9p remaining valid. 

A key physical consequence of a finite a is that, in 
a phase where atoms are condensed, a finite Feshbach 
coupling [that mathematically acts like an ordering field 
on the molecular condensates; see (|4.4bp ] scatters pairs 
of condensed atoms into a molecular BEC. Equivalcntly, 
it hybridizes states of a pair of atoms and a molecule, 
and as a result a finite molecular condensate is always 
induced in a state where atoms are condensed. Conse- 
quently (much like an external magnetic field eliminates 
the distinction between the paramagnetic and ferromag- 
netic states), as anticipated in Sec. IIIH a finite a elim- 
inates the ASF phase, replacing it by the AMSF. A fi- 
nite Feshbach coupling thereby converts the ASF-AMSF 
transition into a crossover, 

^2 = Ml — , Ml > (ASF-AMSF crossover), (4.15) 
51 

between two regimes of AMSF with low and high density 
of a molecular condensate, with the quantitative distinc- 
tion and crossover between these becoming sharp in the 
small a limit. This is illustrated in Fig. [51 For simplicity 
of notation, from here on, both of these regimes will be 
referred to as simply ASF. 

For small a the value of the atomic and molecu- 
lar condensates throughout the ASF phase can be esti- 
mated from the saddle point equations (|4.4p . As for the 
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FIG. 8: The two possible i-ii-fj,2 mean field phase diagrams 
for a special case of a vanishing Feshbach resonance coupling, 
a = 0. (a) For gig2 > 312 there are four distinct phases 
separated by four second order transition lines, meeting at a 
tetracritical point, (b) For gig2 < 512 there are only three 
distinct phases, meeting at the bicritical point. The uniform 
AMSF phase, exhibiting both atomic and molecular super- 
fluidity, is unstable and is replaced by a direct first order 
ASF-MSF transition denoted by the hatched double line. On 
this line the dominant atom-molecule repulsion drives the sys- 
tem to phase separate into ASF and MSF regions. The other 
two transitions remain continuous. In (b) the dashed lines 
fJ-i = /i25i2/(?2 and fi2 = A'1512/yi denote spinodals. 



case of a = 0, for fn > one has ^10 ~ \/\f^i\/9i, 
which (through the atom-molecule repulsion 312) acts 
to shift the effective molecular chemical potential to 
fj,2^ = A^2 — ^J^l9l2/9l■ Vanishing of jjl^ defines the 
ASF crossover line (|4.15p . To the left of and above this 
crossover boundary, the effective molecular chemical po- 
tential is negative and the molecular condensate is small. 




AMSF crossover ASF 



FIG. 9: Schematic illustration of the crossover from small 
molecular condensate fraction (in a region in which n20 = 
|*I'20p would vanish exactly for q = 0) to a large molecular 
BEC. The horizontal axis represents a path, parameterized 
by /xi and /i2, which intersects the crossover line in Fig. 1101 
The estimates for the exhibited sizes of n20 deep in the AMSF 
phase, and in the vicinity of the crossover line, follow from the 
discussion surrounding equation (|4.20p . 



It is induced to be finite, via (|4.4b[) . only by virtue of a 
finite Feshbach resonance coupling to the atomic conden- 
sate. This gives 



*20 = 



a*?o 



2(-M2 + gi2|*io| 



+ 0(a2) 



(4.16) 



2(512/^1 -51M2)' 



tor /ii > 0, /i2 < ^1 . 

eff 



On the other hand below the crossover line, ^2 > 0, a 
would-be spontaneous (for a = 0) molecular condensate 
is only weakly modified from its a = value 



I* 



20 1 



52 

51M2 



512M1 



5152 



for /ii > 0, /i2 > Ml 



(4.17) 

512 
51 ' 



In the intermediate regime, in the vicinity of the crossover 
line itself, one expects ^20 still to vanish with a, but 
more slowly than linearly. Thus, if |vl/2o| ^ Q^/512, then 
(|4.4ap becomes /ii — 5i|*ioP ~ 5i2|*20p- Substituting 
this into the second term on the right hand side of (|4.4b[) 
one obtains to leading order the relation 



A = t|^ 



20 1 



I* 



in which 



A = 



a/ii 



20 1 



512M1 -51M2 



2(5152 -5^2)' 



5152 



5?2 



(4.18) 



(4.19) 



are the scaled Feshbach coupling and deviation from the 
crossover line, respectively. The solution may be ob- 
tained in the scaling form 



|^2o| = A'/'^(r/A2/3) 



(4.20) 
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in which the scahng function X{x) is the solution to the 
cubic equation 



l = xX{x)+X{x)^ 



(4.21) 



with X{Q) — 1. Close to the crossover line, where 
|r|/A^/^ — 0(1), one sees that \l/20 is of order 
(al^PioP)^/"^. For large positive r/A^/^ one enters the 
linear scaling regime where |vl/2o| ~ A/r, which is con- 
sistent with (|4.16p so long as the constraint |^2o| ^ 
a/512 is obeyed. This leads to the condition ((712/^1 — 
911^2)/ 912 fJ'i <C 1 on the deviation from the crossover 
line. The crossover behavior is sketched in Fig. [S] 

It is clear from (|4.4p that ^10 = is still a solution, so 
that the N-MSF transition line (14. 9p is not modified by 
the Feshbach resonance coupling. Then equation (j4.4b|) 

still gives 1^20 1 = \/ 1^2! 92 in the MSF phase. However, 
the subsequent MSF-ASF transition boundary is modi- 
fied, as a finite a shifts the effective chemical potential 
of the atomic condensate (in addition to the shift due to 
the atom- molecule repulsion, gyi) to 



PI 



/^i -gi2|*2 



aW-}. 



(4.22) 



with the MSF-ASF transition located by [I'f = 0. Using 
1*20 1 = vaWss of the MSF one obtains an estimate of 
the MSF-ASF transition boundary 



Ml 



512 

52 



-M2 



(MSF-ASF transition line). 



(4.23) 
For small [ii (specifically, < /i2 <C c^g^lg^-i), the sec- 
ond term on the right hand side dominates, and the 
boundary displays a sharp square-root singularity into 
negative values of [i\ (near the origin preempted by a 
first order transition: see below) illustrated in the phase 
diagrams for 9\gi > 5^2 and 5152 < 512 in Figs. [10] and 
nn respectively. In the opposite limit (/i2 ^ 0^^52/512)1 
the boundary asymptotes to the a = phase boundary, 

Another important consequence of a finite Feshbach 
resonance coupling is that for small chemical potentials it 
drives the N-ASF and MSF-ASF transitions first order. 
The somewhat technical calculation of the corresponding 
first order phase boundaries, illustrated in Figs. [TO] and 
im are relegated to Appendix [BJ Here a more approx- 
imate, but more transparent, analysis is presented. To 
this end, one can use an approximation to (j4.4bp 



* 



20 



2|M2 
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IOj 



(4.24) 



valid for sufficiently negative ^2 and fii, to eliminate ^20 
from Timf ['I'lOi *2o] in favor of ^iq. The resulting energy 
density in the normal state is well approximated by 



n 



mf ' 



-Mil* 



101 



-^1.91 



2|M2| 



I* 



101 



"'^"1*10^. 



4^1 



(4.25) 
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FIG. 10: Mean field phase diagram in the /ii-/i2 plane for 
a finite Feshbach resonance coupling a 7^ and gig2 > gi2- 
In contrast with the a = limit shown in Fig. HJa), a fi- 
nite Q eliminates the distinction between the ASF and AMSF 
phases, converting the ASF-AMSF transition, indicated by 
the dashed line, into a crossover. Feshbach resonance scat- 
tering also strongly modifies the MSF-ASF phase bound- 
ary, and for small chemical potentials drives the N-ASF 
and MSF-ASF transitions first order (indicated by hatched 
curves) , with the first order section terminated by two tricrit- 
ical points. The point where the three phases meet, and the 
continuous N-MSF phase boundary terminates at the first 
order boundary, is a critical endpoint. 



Clearly, for /ii < and sufhciently large negative effec- 
tive quartic coupling, u^ = gi — a^/2|/i2| a secondary 
minimum at ^^lo ^ develops, that can compete with 
the normal state 'I'lo — 'I'20 = minimum. It is easy 
to show that the corresponding ASF state minimum be- 
comes degenerate with the normal state at a critical 
value fj,1 = — 3m|/16w6 {uq = a^5i2/4/^| is the effective 
l^iol^ coupling), that translates into a first order N-ASF 
boundary 



Mf 
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IM2I51 



(4.26) 



as illustrated in Figs. [TO] and [TT] On the other hand, 
for sufficiently small /i2 > 0, the relation between ^^20 
and '5 10 following from (j4.4bp becomes (keeping only the 
52|*20p term inside the parentheses on the right hand 
side) , 
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20 
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(4.27) 
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FIG. 11: Mean field phase diagram in the /ii-/i2 plane for a 
finite Feshbach resonance coupling q 7^ and gig2 < 312- In 
contrast with the a = limit in Fig. [8jb), finite a eliminates 
the distinction between the ASF and AMSF phases. Fes- 
hbach resonance scattering also strongly modifies the MSF- 
ASF phase boundary, and for small chemical potentials drives 
a segment of the N-ASF and MSF-ASF transitions first order 
(indicated by hatched curves), with the first order section ter- 
minated by two tricritical points. The point where the three 
phases meet, and the continuous N-MSF phase boundary ter- 
minates at the first order boundary, is a critical endpoint. At 
large positive p's the MSF-ASF transition retains its a = 
first order character, separated from the continuous section of 
this transition by another tricritical point. 



which, when inserted into Timf, Eq. (|4.1[) . determines the 
first order MSF-ASF phase boundary, in a way detailed 
in Appendix [Bl 



V. DILUTE BEC LIMIT 

Thus far, the phase diagram in the /ii-/i2 plane, has 
been studied by treating the atomic and molecular chem- 
ical potentials as independent tuning parameters. As 
seen, within a mean field approximation, the tempera- 
ture then plays no apparent role. 

However, to make a direct contact with trapped de- 
generate atomic gas experiments, where it is the total 
number of atoms N and the detuning v that are varied, 
one needs to eliminate the chemical potentials fi^ in favor 
of the atom density n, detuning i> = 2/11 — fi2, and tem- 
perature T. For an interacting system, this is a nontrivial 
change of variables that can usually only be carried out 
perturbatively However, in the dilute limit, appropriate 
to atomic gas systems, the transition out of the normal 



state can be treated by ignoring weak atomic interactions 
(with corrections in powers of na^ and 7) , thereby reduc- 
ing the problem to an easily calculable BEC limit i^ The 
system then reduces to two independent ideal Bose gases, 
coupled only through the overall constraint of fixed den- 
sity n^ni + 2n2, Eq. (PTTj) . 



A. Bulk N-ASF and N-MSF BEC transitions 

In the noninteracting limit, for a bulk (uniform) system 
the free energy and atom density in d spatial dimensions 
are easily calculated and are given by^ 



f'^w-^H 
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-P{eka—fJ.cr) 



PA. 



k,cr=l,2 
T cr=l,2 



(5.1) 



o-=l,2 



A^ E ^ 






k,cr=l,2 



(5.2) 



where e^ = ft fc /277ia- (with mi — m and 7712 = 2m) are 



the single particle energies, Zc 



— pPt^a 



are the fugacities. 



At = h/\/2TTmkBT is the (atomic) thermal de Broglie 
wavelength, and 



9a{z) 



1 



r(a) Jo 



^Q-l 



dx- 



1 



y- 



(5.3) 



is the extended zeta function. 

For positive detuning, v > Q, atoms (being less energet- 
ically costly than molecules) condense first at a critical 
line /ii = 0, where ^2 — 2//i ^ v = —v < Q. The corre- 
sponding critical temperature Tci{i') for N-ASF transi- 
tion is easily determined by the fixed density condition 



A, 



Tel 



C(d/2)-t-2('^+2)/2^^^^(^-/5.i.^ 



(5.4) 



with C(q^) = 5q(1)- As usual, the transition exists only 
for d > 2. In the far detuned limit, v/ksTd ^ 1 
(z2 ^1), the molecular population is exponentially sup- 
pressed, the second term in (|5.4p may be neglected and 
one obtains the standard BEC result 



rpac 



2'KmkB 



C{d/2) 



2/d 



(5.5) 



that is approached exponentially as e '^'=1". In the oppo- 
site limit, < v/kBT <^ 1, the expansio n^^i^° 



g„(e--)=r(l-a)x"-i+£ 



C(a - "•) 



{-xT (5.6) 
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may be used to obtain 



T, 



cl 



Tco 



1 



2(d+6)/2p (^-d-\ 



did - 2) [1 + 2(''+2)/2] C(d/2) yksT, 



(5.7) 



valid for the range 2 < d < 4 of interest to us, where 
. n r ^ 2/rf 



TcO — 



27rTO^ 



1 + 2('^+2)/2 ({d/2) 



(5.8) 



is the transition temperature at zero detuning, v = 
corresponding to the critical endpoint in Fig. [21^^ 

Similarly, for negative detuning, v < 0, molecules are 
energetically less costly and therefore condense first. The 
corresponding N-MSF critical line, given by /i2 = 0, with 
fii — V jl < 0, using the fixed density condition translates 
into a Tc2{v), determined implicitly by 



n = L 



-d 



2('^+2)/2c(d/2)+g,/2(e''^^''/')l. (5.9) 



are noteworthy. The normalized transition temperatures, 
Tccr{v)/Tco, give the corresponding phase boundaries (a 
function of h'/ksTco) displayed in the phase diagram in 

Fig.m 

The thermodynamics of this dilute Bose gas mixture 
above the transition temperature (i.e., inside the normal 
state) can be obtained by using the atom number con- 
straints (|5.2p to express the chemical potentials fii = fi, 
/i2 = 2/i — 1/ as functions of temperature and detuning. In 
the neighborhood (above) the N-ASF and N-MSF tran- 
sitions this can be done analytically using (|5.2p and (|5.6p . 
To leading order, for 2 < d < 4 and |/i| <C i' one obtains 
near the N-ASF line: 



T 



d/2 



r(2^)(/3|M|)('^-2)/2 



'~C(d/2) + 2('^+2)/25,/2(e-/3'')' ^^- -"> 
while in the neighborhood of the N-MSF line, defining 
5iJL = jl — v/2 ^ \v\, one obtains 



In the far detuned limit the N-MSF transition tempera- 
ture approaches 



T'O 



/l2 



AirmkB 



2C(d/2) 



2/d 



(5.10) 
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T 



d/2 



2'^r(^)(/3|<5/i|)( 



2('^+2)/2C(d/2)+5,/2(e'3-/2)- 



Finally, for ly — one finds 



(5.14) 



exponentially as e ""'^''^'^ while in the small detuning 
limit li'l/kBT ^ 1 it approaches Tco according tc^i 



lc2 



2(6-d)/2p (^4^^ 



d{d - 2) [1 + 2('i+2)/2] ({d/2) \kBT, 



The resulting ratios 

rpOO 
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Tco 

T'CXD 



-^ j^ 2(d+2)/2 



2/d 



1 + 2-(''+2)/2 



Tco 

Tc2{,-v) - TcQ 



2/d 



2'^, 0< 



IcO 



(5.11) 



ksTco 



< 1. (5.12) 



Tco 

T 



d/2 



(i + 2'^)r(^)(/3H)('^-2)/2 

C(rf/2) (1 + 2('^+2)/2) ' ^^^^' 



Therefore, for a = 0, 1, 2, corresponding to i^ = 0, z^ > 0, 
and v < 0, respectively, one obtains 



ksTa 



An 



T 

J- err 



-|2/(d-2) 



- 1 



(5.16) 



where T^a are transition temperatures evaluated at the 
given values of i/, n, the chemical potential deviations are 
given by 5^i — J/io — /i, 5^2 — Sfi, and the amplitudes 
are 



A2{l^) 



j d{d-2) 

WW) 

{ did -2) 
\ 2d+^T .' ' 



2/{d-2) 



C(d/2) + 2('^+2)/2.g,/2(e-/'=«^-) 



2/(d-2) 



4-d\ 
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4r(^) lT2d~ 
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(5.17) 
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For d > 4 the exponent 2/{d— 2) sticks at unity, so that Sfi varies hnearly with the temperature deviation/? 

For d — 3, using r(l/2) = y^ and C(3/2) — 2.612, all of the above results reduce to those quoted in the Introduction 
and summarized by the phase diagram. Fig. [H 

The rest of the thermodynamics in this noninteracting limit now follows in a standard fashion, leading to the familiar 
Gaussian model critical behavior. For example, the atomic and molecular condensate densities below their respective 
normal-superfluid transition temperatures Tca{v) are easily computed. For v > Q the atomic chemical potential 
/ii = /^ vanishes before the molecular one ^2 = 2/i — ;/ and for T < Tci{i') an atomic condensate nio(T, v) = |^ioP 
develops with 



nio 



tV" C(rf/2) + 2^g,/2(e-A^^) 
Tci) C(d/2) + 2^g,/2(e-'Vfc-^-) 



, for v>Q, T <Tci{y), 



(5.18) 



as the gas transitions to the ASF in the BEG limit. Glose to Tci{v)^ the atomic condensate (|5.13p grows linearly with 
reduced temperature 



nio(T, :/) ^ nil 



T 



Tci{v) 



for i/>0, T^T-^{v), 



consistent with the expected order parameter exponent (3 — 1/2. 

Similarly, for v < Q the molecular chemical potential fi2 — ^fi ~ v vanishes before the atomic one fii 
T < Tc2{v) a molecular condensate n2o{T,i/) — |^20p develops with 



n2o 



—n 
2 



1- 



T 



d/2 



2^Cid/2)+g,/2ie-/^''-^) 
I 



for u <0, T <T^2{v), 



(5.19) 
/i, and for 

(5.20) 



as the gas undergoes a transition into a molecular BEG 
(MSF). Again, close to Tc2{v), the molecular condensate 
(|5.14p grows linearly with reduced temperature 



n2o{T,v) ^n\l 



T 



Tc2{y) 



for V <Q, T^T-^iv) 



(5.21) 
consistent with the same order parameter exponent f3 = 

1/2. 

Glearly, the above expressions for the weakly interact- 
ing limit are quite close to standard ones for a single- 
component Bose gas, reducing to a purely atomic and 
molecular BEG for the far detuned cases, Ivl/ksT ^ 1. 
There are, however, nonstandard contributions to Tca{v) 
and n(jQ{v,T) arising from the contribution of the sec- 
ondary, off-resonance, bosonic component that is gapped 
out for V ^ Q. For example, for 1/ < 0, upon warming to- 
ward Tc2{i'), the molecular condensate is reduced due to 
both the conventional mechanism of thermal excitations 
of molecules out of the molecular condensate, as well as 
the depairing of molecules into thermally excited bosonic 
atoms, with the latter special to a Feshbach-resonant sys- 
tem. 

Because of the suppression of Tccr{v) near v — Q ior 
TcQ < T < T^{i/), the gas is expected to undergo a 
sequence of ASF ^ N ^ MSF transitions upon lowering 
of v (see Fig. [5]). For T < Tco the transition is a direct 
ASF — > MSF one, that for this noninteracting limit is 
first order. The condensate densities are undefined right 
on the critical line i^ = 0, T < Tco, and the noninteracting 
approximation becomes particularly questionable there. 



B. N-ASF and N-MSF BEC transitions in a trap 

The above results are straightforwardly extended to 
the experimentally more relevant case of a harmonic trap. 
The modifications due to the trapping potential can all be 
incorporated through the change in the density of states. 
For an isotropic harmonic trap (easily extendable to an 
anisotropic trap) the single particle energy spectrum 



fihJoini 



n2 



nd) 



(5.22) 



is linear in n = ^^ rii and exhibits a well-known degen- 
eracy that, for large quantum numbers of interest to us, 
in the macroscopic limit ksTca- 3> hcoo, is given by 



D{s) 



1 



{d-iy.{hLJoy 



(5.23) 



Note that £„ is actually cr independent, i.e., the trap 
frequency uq is the same for atoms {a — I) and molecules 
(a — 2). This is a good approximation in the physically 
relevant limit of the size do of the closed-channel molecule 
being much smaller than the trapped cloud size. 

In the thermodynamic limit the sum over single- 
particle states appearing in (|5.ip and (|5.2p can be re- 
placed by integration over energies e weighted by above 
density of states, giving 
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(5.24) 
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Paralleling the above calculations for the uniform system, 
from this one obtains all the relevant quantities for the 
trapped system. Specifically, the transition temperatures 
Tcaii^) are implicitly given by 



N = 
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(id) + 2gd 



2C{d)+gAe 
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(5.25) 



.(5.26) 



These can be solved in the asymptotic regimes of small 
and large detuning, giving 
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Tr, 
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\v\ < fcsTco 
\v\ > fcsTco, 



(5.27) 

with a„ = 2C{d-l)/'Sa^dCid) and b^ = 2/a^dC{d). The 
transition temperatures T^ in the limit of asymptotically 
large positive (cr = 1) and negative {a = 2) detuning 
{\v\/kBTca 3> 1), and at the tricritical point Tcq (v = 0), 
are given by 



Tco = hujQ 



N 



add) 

N 



Kid) 



l/d 



1/d 



(5.28) 
(5.29) 



The latter is approached linearly with reduced detuning 
from either side, in any dimension d > 2. 

As in the bulk case above, and in the well studied single 
component trapped Bose gas, here too one can easily 
compute the number of condensed atoms and molecules 
below the transition into the ASF and MSF states. This 
is determined by the extension of the total atom number 
constraint (|5.24p to include the condensates No-q: 
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ksT 
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(5.30) 



The analysis of these equations closely follows that 
for the bulk BEC of the previous subsection. Below 
the transition into the ASF and MSF phase one can 
straightforwardly compute the number of atoms N^a — 
J (ir|^o-o(r)P in the corresponding condensate. For i^ > 
the atomic chemical potential fii = n vanishes before the 
molecular one /12 = 2/i — i^, and for T < Tci(^) the molec- 
ular condensate iV2o = and a finite atomic condensate 
Nio{T,h') develops, given by 



N, 
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C{d) 


+ 2gd{e 


-v/kBT\ ' 


fc 


(id) 


+ 2gd{e- 


t^/ksTai) 



for v>Q, T <fci{v). (5.31) 

For V < the molecular chemical potential /Lt2 = 2/i — j/ 
vanishes before the atomic one /zi = /i, and for T < 



Tc2{v) the atomic condensate iVio = and a finite molec- 
ular condensate N2o{T,v) develops, given by 
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2C(d)+5d(e'^/2'=BTe2) 
for t/ < 0, T < fc2(j^), (5.32) 



Just below the transition temperatures Tea the conden- 
sate growth is of the expected linear in T form, charac- 
teristic of the order parameter exponent (3 ~ 1/2. Also, 
for a far detuned gas, \v\/kBT ^ 1, the above results 
reduce to the standard single component BEC behavior. 
The advantage of a trapped system is that, as in the 
case for an ordinary single-component trapped conden- 
sate that exhibits a striking narrow BEC peakf22ii2i here 
too we expect ASF and MSF condensates in a trap to 
display clearly identifiable BEC peaks. As discussed in 
the Introduction and illustrated in Fig. [H provided that 
atoms and molecules can be imaged separately, the ASF 
should be easily identified by atomic and molecular BEC 
peaks,— while the MSF is identified by the presence only 
of a molecular BEC peak. In a harmonic trap at low 
temperature, T ^ T^a-, the density profile of the cloud 
is dominated by a narrow Gaussian tr-condensate peak, 
with the width given by the quantum oscillator length 
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(5.33) 



This should be easily distinguishable from the high- 
temperature, classical Gaussian density profile (coming 
from the Boltzmann distribution) with the much wider 
width set by the thermal oscillator length 



rcrT 



\m„u^l 



ksT 2kBT 

= raO\ —i > 7V0- 



hw(\ 



(5.34) 



The full atomic (whether free or bound into molecules) 
density profile nir) at arbitrary temperature is easily cal- 
culated for a nonintcracting gas. It is given by 



*W ^E^^-^W' 



(5.35) 



consisting of atomic and molecular contributions, in the 
BEC limit tied only by a common chemical potential 
/i determined by the overall particle number constraint 
(I5.30p . As derived and analyzed for a single Bose com- 
ponent in App. [ni these in turn are given by 
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where (p-aai'^) are harmonic oscillator eigenstates and 
p°^'^(r, r;/3?iwo) is the diagonal element of the single- 
particle density matrix for a harmonic oscillator with 
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mass m^. In 3d it is given by 



pr-(r,r;/3?^o) = 
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hujo/ksT > 1, 
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(5.37) 

(5.38) 
(5.39) 



is the finite-temperatur e "oscillato r length" that reduces 
to the quantum one ^ h/ {nia-ujo) , Eq. (|5.33p . at low T, 
and the classical (thermal) one, Eq. (|5.34p . at high T. 

The spatial profile of the cr-density, n^{r), is deter- 
mined by the ratio of the chemical potential /Xo- to the 
trap level spacing Hluq, with former in turn determined 



by the temperature through the total atom number con- 
straint. At high T ;:g> Tea (where the gas is nondegen- 

erate), such that < -^^ w -kBT\Ti[(j^^ N] w 

3kBTln{T/Tccr) ^ ksT, the result is a purely classical 
thermal (Boltzmann) distribution, 
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(5.40) 

with only of order unity occupation of the lowest oscil- 
lator n — state and a vanishing "condensate" density 
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As T is lowered further, approaching Tea from above, 
the magnitude of the chemical potential drops below T 
(remaining negative) and the boson density profile devel- 
ops a small r non-Boltzmann peak structure even above 
T ■ 
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The linear in r cusp is rounded on the length scale below ro-oVTv/^^- 

Finally at an even lower T < Tea, \^-a\ drops below the level spacing, |/icr| ^ fv-oo, and the density profile changes 
dramatically, developing a bimodal distribution na{r) = naT{f) + '^cro(^) (see Fig. [1}, that consists of a broad (width 
To-t) thermal part 
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with a small r cusp (rounded by rao) and large r Gaussian 
tails, together with a narrow (width To-o) condensate part 
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In dESD, 
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has been defined, while in (|5.44p 
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is the number of condensed bosons, given by ()5.3ip and 
(I5.32p when the total atom number constraint is taken 
into account. 



VI. ELEMENTARY EXCITATIONS 



Having established the approximate nature of atomic 
and molecular supcrfluids, consider next the study of 
their excitations. On general grounds, as required by the 
Goldstone's theorem, one expects one collective gapless 
(sound) mode in each of the ASF and MSF phases, asso- 
ciated with spontaneous breaking of global U{1) charge 
(phase- "rotation" ) symmetry. In the MSF it is associ- 
ated with the phase 02 of the molecular (two-atom) con- 
densate, ^20 J while in ASF it corresponds to in-phase 
fluctuations of the phases 6i and 02 of the atomic and 
molecular condensates. 

In addition, there are three2^ gapped excitations in 
each of the superfluids. In the MSF these are asso- 
ciated with atom-like (squeezed by Feshbach resonance 
coupling to the molecular condensate) quasiparticle ex- 
citations (accounting for two modes, ipi,tpi) ^nd molec- 
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ular density fluctuations (fluctuations in the order pa- 
rameter magnitude |^2ol)- In the ASF one gapped mode 
corresponds to out-of-phase fluctuations 29i — 62 of the 
atomic and molecular condensates (gapped by the Fesh- 
bach resonance coupling a), and two others are atomic 
and molecular condensate densities (fluctuations in the 
order parameter magnitudes |^io| and |^2o|)- 

As will be seen below, the MSF-to-ASF transition is ac- 
companied by closing of the gap for atom-like quasiparti- 
cle excitations. However, this mode remains gapless only 
at the MSF-ASF critical point, and is replaced by an- 
other gapped mode (associated with out-of-phase phase 
fluctuations of the two order parameters) that emerges 
inside the ASF. As discussed in Sec. IIIIl this is consis- 
tent with the Goldstone theorem as (due to the Feshbach 
resonance coupling) it is only a discrete {Z2) symmetry 
that is being broken at the MSF-to-ASF transition and 
as such leads to no new gapless modes. 



One obtains 
H2 



E 

k.cr 



eko-aj^o-'^ko- + -z (AcrflkCTa-ko- + h.^ 



^ (tialj^ak2 + ^2akia^k2 + ^■ 



where the coefficients are given by 

Eki = £ki -Mi + 25i|*iop+5i2|*2or 
ew2 = £k2 -Ai2 + 252|*20P+5i2|^ioP 



Ai = 51*1 



a^ 
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A2 = 52*2 



h = 512*10*^0 -"^10' *2= 512*10*20, 

with single particle energies Eka — h'^k^/2ma- 
1. MSF phase 



(6.4) 



(6.5) 



A. Bogoliubov diagonalization 

Bogoliubov theory provides an asymptotically exact 
description of the low energy excitations in a dilute Bose 
fluid, not too close to the transition lines. Focusing on 
quadratic fluctuations, it ignores interactions between 
quasiparticles and, among other things, misses the possi- 
bility for their decay. The method proceeds by expanding 
the field operators about the mean field solution (equiv- 
alently, a coherent state of k = fields labeled by *(jo): 



Consider first the excitations in the molecular super- 
fluid, characterized by a finite *20 and vanishing ^iq. 
As a result, the cross terms ti and t2 vanish, and the 
atomic and molecular terms can be diagonalized inde- 
pendently. From (|4.10p . the mean field order parameter 
is given by ^20 = ^^2/52 = \/(2^ - v)/g2 (chosen real 
and positive for simplicity — more generally any phase 62 
can be absorbed into the operators via the redefinition 
flkcr -^ e~*'^^^/^akcr). It is straightforward to verify that 
the Bogoliubov canonical transformation 



i/'<x(r) = *„o + ?!'a(r). 



(6.1) 



Ok^ = Wk^7k^ - Wko-7lk^ 

7ko- = Wkcrflkcr -f Wkcra„k<T 



(6.6) 



and keeping terms in the Hamiltonian only to quadratic 
order in the small deviations (pa. In the molecular super- 
fluid state *io = 0, so ^1 = -01. Substituting (|6.ip into 
(|2.ip one obtains 



H = H,^i + H2 + 0{4>l,c^'i) 



(6.2) 



in which H^i = _ff ['J'o-o] is the mean field approximation 
(|4.ip to the ground state energy. The absence of terms 
linear in excitations <}>„ is guaranteed by the condition 
that ^ctO is an extremum of the mean field free energy 



aO 



0. To quadratic order, this is equiva- 
lent to the requirement ((/)cr(r)) = Oj ^^i^^ For a homo- 
geneous system [generalization to the trapped case may 
then be accomplished through a local density approxima- 
tion (LDA)] the quadratic Hamiltonian, H2 governing the 
dynamics of fluctuations, can be represented in terms of 
momentum space operators 



0a (r) 



E 

k 



ak.e'''•^ cj^liv) 



E 



"kcr^^ 



(6.3) 



to new bosonic creation and annihilation operators 7^^^, 
7k(T, with real, positive coefficients given by 

1 / £k(T 



'■kff 



pMSF 
^kcr 



1 + V- 



kcr 



2\E, 



kcr 



kCT 



IA.P, 



leads to the diagonal form 
SHmsf = E ^i 



MSF 
kcr 



k,o- 



[ll 



7k<T - t'ka 



(6.7) 
(6.8) 

(6.9) 



The diagonalized Hamiltonian, Ji^MSFj governing exci- 
tations in the MSF naturally separates into "atom-like" 
((T = 1) and "molecule- like" {a — 2) contributions, 
with corresponding (explicitly positive) excitation ener- 
gies E^^^ and condensation energy 

cpMSF _ V^ pMSF,,2 _ \^^(B pMSFN 

O^cond = - Z^^-ka ^kcr,- "2^ 2*^ '"^ ~ kcr )■ 

k,(7 k,(T 

(6.10) 
The latter lowers the energy of the MSF below that given 
by the mean field condensation energy value, -ffmf • 
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In the normal phase, ^20 = 0, one obtains \a — 0, 
Eko- — ^ko-— MtT) yielding Wko- — and Uko- — 1- One there- 
fore recovers the original atomic (aki) and the molecular 
(ak2) operators as true (to quadratic order) excitations in 
the normal state, with corresponding free single-particle 
spectra eku — /iu- 



2. ASF phase 

It is clear from the structure of the Hamiltonian SH asf 
in the ASF phase (most notably the finite values of the 
ti and i2 couplings), that in addition to the usual Bo- 
goliubov mixing between particles and holes, a true ex- 
citation is also a mixture of an atom and a molecule. 
Physically, this is a reflection of a coherent scattering (by 
the Feshbach and atom- molecule density interactions) of 
atoms and molecules mediated by their respective con- 
densates. The Bogoliubov theory for the ASF phase is 
handled most simply by first converting from creation 
and annihilation operators to corresponding "position" 
and "momentum" operators (canonically conjugate "co- 
ordinates" , that are Fourier transforms of Hermitian field 
operators))^ 

1 



flko- 


- ^(gk^+«Pka) 


«L 


- /2^^-'^'^ *P-k<T 


g-ko- 


= (zL, P_k. =pL 



(6.11) 
with the only nonvanishing commutation relations being 

[qka,P-k'a']=-iSkk'Saa'- (6.12) 

By substituting (|6.12p into (|6.4p one obtains 



6H 



ASF 



E 



^^k'^-^(^~ki+£k2) 



in which the 2x2 matrix structure is defined by 



Qk = 

Pk = 

Qk = 



(6.13) 



9kl 

9k2 


j , Pk = 


( Pkl 


Ekl 


-Ai 


tl 


-t2 


h 


-t2 


£k2 


^\2 , 


Ekl 


+ Ai 


ii 


+ t2 


h 


+ t2 


£k2 


+ A2 



(6.14) 



In deriving (|6.14p the symmetry e_kcr = eko- has been 
used, and ti, t2, Ai, A2 have been taken to be all real (or, 
equivalently, their phases absorbed into redefinitions of 

Skcr). 

One seeks a (real) linear transformation 



Pk = -4kPk, Qk = SkQk 



(6.15) 



which diagonalizes SH^^^. The canonical requirement 
that the transformation preserve the commutation rela- 
tions (|6.12p . i.e., that 



[Qk,PL] 



[qk,Pk'] ^j^k'i, 



implies that 



Bi = A-, 



(6.16) 



(6.17) 



Thus, the transformation Bk should simultaneously di- 
agonalize V^^^ and Qk- Without loss of generality this is 
equivalent to demanding that 



^k^k'^k^a, SkSk^k^fk, 



(6.18) 



in which fk = diag[(i?^f^)^, (i?^^^)^] is diagonal, con- 
taining the squares of the Bogohubov energies (see be- 
low). It follows that 



^w 'PkQkSk = fk, 



(6.19) 



so that £k is obtained by diagonalizing PkQk- The 
squared energies are therefore solutions to the eigenvalue 
equation 

det[{E^f)H - T'kQk] = 0. (6.20) 

The solutions to the resulting quadratic equation in 
(i?ASF)2 are 



(£^i^r)^ = ek± 



dl 



J1)J2) 
^k "^k 



(6.21) 



in which the upper sign corresponds to ct = 1, the lower 
sign to a — 2, and the various parameters are defined by 



T'kQk = 



Gk + dk 



ek 



Ml 



(2) 
Ck 

ek - dk 

2 I ~2 



Al + ^k2 ^ A2) + tl — ^2 



(6.22) 



"^k = 2'^^ki - Al -ek2 + A2) 

4'^ = (tl-t2)(ekl+Ai) + (ti+t2)(£k2-A2) 
Ck^ = (tl+t2)(ekl-Ai) + (ii-t2)(£k2 + A2). 



It is easy to check that the MSF results (|6.7p and 
are recovered when <i = <2 = 0. 

The columns bka of Sk = (bkibk2) are the eigenvec- 
tors of T^kQk and take the form 



1 



JkCT — 



.(2) 



A^ka V ek + rfk - £^k<T 



(6.23) 



in which the normalization A/kcr is chosen so that 
bL^k 'bk. = 1. 

The quadratic Hamiltonian takes the form 

6HaSF = lY.^P-UaPka + {E^rfQ-).aQka 



E 



^ASF-t 



1 



^ka 7k^7k<T + 2(^k^ 



£ko 



- Ekcr) 

(6.24) 
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in which the new bosonic raising and lowering operators 
are given by 



7ka = ^= I E^l^Qka 



7L 



V2 ^ 



E, 



1/2 



ko- W-kcr 



E, 



1/2 

ko- 



k(T 



P-v 



(6.25) 



These may be reexpressed in terms of the original raising 
and lowering operators via 



1 



ak = 2V^kfk 
1 
2 



;(Skfk^''-Su"4/l7lk 



71/4R-1 



'1/4rT 



7k = 2(CX-^+^k"Xjak 

+ ^(^k''Sk^-^k^/X)a^-k (6.26) 

in which ak,a_j.,'yk,'7lk ^^^ ^-H column vectors defined 
in the natural way, consistent with (|6.14p . 



B. Acoustic and gapped modes 

Consider now the Bogoliubov excitation spectra, 
and (|6.2ip in more detail, ft will be shown that in both 
phases there is indeed one acoustic mode and one gapped 
mode (in addition to two other less interesting gapped 
modes^^), as required by general principles discussed in 
the beginning of this section and in Sec. IfffI As previ- 
ously indicated, the MSF phase the acoustic mode cor- 
responds to long wavelength fluctuations in the phase 
of ^20 J while the gapped mode is associated with pair- 
breaking fluctuations of molecules into two atom-like ex- 
citations, with spectral gap corresponding to a renormal- 
ized molecular binding energy. In the ASF phase the 
acoustic and gapped modes correspond to in-phase and 
out-of-phase fluctuations of \l/io and \l/2o, respectively, 
with the gap in the latter governed by the Feshbach res- 
onance coupling a. 



1. MSF phase 

The MSF quasiparticle spectrum (|6.8p appearing in 
(|6.9p (and summarized in Fig. [3]), may be written in the 
form 



E),a = V(ek<T+£^-)(ek<T+£^+) (6.27) 

in which the (positive) energies ea-± are given by 

ei± = -iy/2 + (.gi2-.g2/2)|*20P±a|*2o| 

£2+ = 252|*2o|', £2- =0. (6.28) 



The molecule-like branch {a = 2) is gapless (consistent 
with Goldstone's theorem), having an acoustic spectrum 
Hc^^^k at small k with sound speed 



ijMSF 



MSF 



= |*20|Vff2/m2, 



(6.29) 



corresponding to collective, long wavelength oscillations 
of the molecular condensate. The spectrum crosses over 
to a particle-hke EMp^ w ek2, for k^^l^ > 1, where 



fMSF 
>coh 



2m2C^ 



MSF 



(6.30) 



is a coherence length beyond which superfluid behavior 
sets in: the collective superfluid response-''^ dominates 
disturbances with wavelength longer than ^coiu while the 
microscopic single molecule response^^ dominates those 
with shorter wavelength. This length ^^^^ oc l/|'I'2o| ex; 
1/^7*2 diverges as the normal phase boundary, 112 — 0, 
is approached. 

In contrast, the atomic- like branch has a gap 



E. 



MSF 



E, 



MSF 



V^i+ei- 



(6.31) 



which closes with increasing ly precisely on the ASF-MSF 
transition line, the latter being equivalent to the condi- 
tion ei_ = 0. This leads to the critical detuning 



= -(52- 2512)1*20!' -2a|^ 



20 1 



2/i 



2^- 



92 

52 /i' 



a^ + 2/Z512 ± a\/aM^4/i5i2 



512 ^ 0, (6.32) 



where the second line requires /i > — a^/4(7i2, and fol- 
lows by substituting |*20p = (2^* — '^)/ff2 and solving for 
z/. The existence of two solutions reflects the reentrant 
behavior as a function of chemical potential seen in Fig. 

m 

At low temperature and for weak interactions, the con- 
densate depletion is minimal, ri20 = |*20p ~ n/2, and 
the critical detuning for the quantum MSF-ASF transi- 
tion is given by 



MT = 0) 



2(52-2312)71 



/2n 



(6.33) 



The behavior of i^ciT) for high temperature [as well as the 
corresponding temperature dependence of the condensate 
"-20 (T) at fixed density n] , illustrated in Fig. [21 will be 
discussed in Sec. IVIDI below. 



2. ASF phase 

In the ASF the extremum conditions (|4.4p allow Sko- 
to be reduced to the forms 

^ki = ^ki + Ai + ai 

ek2 = ek2 + A2 + a2, (6.34) 
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and from (16.51) one has 



t2-ti^ as, 



(6.35) 



with the definitions, 



ai 



2a|*2o| 

1 

— ( 
2 



^2 = ^a|*iolVl*2ol 



as = a|*io| = V"i"2- (6.36) 

Substituting ([O^ - fOS)) into (lO^ one obtains, 

ek + rfk = -a3(ii +^2) + (^ki +ai)(eki + ai +2Ai) 
ek-dk = -a3(ii +^2) + (£k2 + a2)(ek2 + ^2 + 2A2) 



.(1) 

-k 

,(2) 



(£k2 + a2){ti + 12) - a3(£ki +ai+ 2Ai) 



k = (£ki +ai)(ii +^2) -Ck3(£k2 +^2 + 2A2 



(6.37) 



At k = it is easy to verify that 
eo + do = 

and therefore that 




21* 



l^iol 
l^iol 

2|*20| 



201 (1) 



„(2) 



2 _ ,2 , (1) (2) 
So — "0 + Cq Cq . 



(6.38) 



(6.39) 



Substituting these results into (|6.2ip one obtains the ex- 
citation energies at zero momentum, i.e., the gaps: 



pASF 



/^^ 



E, 



ASF 
02 



= V"i(ai + 2Ai) + a2(a2 + 2A2) - 203(^1 + 12) 
= 0, (6.40) 



which confirms the existence of one gapped and one gap- 
less mode2^ in the ASF state. From (|6.5p and (|6.36p one 
sees that a2, a^, ai -I- 2Ai, and hence eo, vanish on the 
MSF-ASF phase boundary where ^I^io = 0. The gap 
therefore closes on the transition line, as expected. Note 
also that if a = one has eq = 0, and the atomic- like gap 
remains closed throughout the ASF phase, as expected 
from the additional spontaneously broken U{1) symme- 
try (separate atom and molecule number conservation) 
and associated Goldstone modes, as discussed at the be- 
ginning of this section and in Sec. IIIII 

The small k (low-energy) behavior of the excitation 
spectra are now examined in the ASF phase, |*io| > 
0, and in the neighborhood of the MSF-ASF transition, 
where |*io| -^ 0, but ^20 remains finite. To this end. 



the 1*10 1 and k dependencies are isolated by writing 
eo = 2a|*io|^ (ge + a- 
do = 2a|*ioP (gd-a- 



*20| 



'20 






.(2) _ 



*20 



4a|*io||*2o| (gf^ +a 
Sek = ek - eo 

= 1*20 kki ( 7. 

5dii = dk- do 

= |*2okkl ( 7d 



'lOl 



^10 I r , 5£kl 



T^d 



1*20 I 

3eki 

!|*20 



.(1) = Jl) 



(1) 



K' ^ c\:'-c)," = \^io\e^.iT 



.(1) 



Sc, 



2) 



.(2) 



.(2 



c\:'-c'o'-\^io\e^.ij^c 



(6.41) 



in which the coefficients 



7, 



5e 


= (51 -512+32/4)1*20! + a/2 


gd 


= (ffi -52/4)1*20! 


(1) 

c 


= (2.gi - 512 )!*2o!+ a/2 = 5e+.gd 


.(2) 
c 


= (52/2 -5i2)!*2o! + a/2 = 5e-5<i 


7e 


= 52!*2o!/2 + a 


Se 


= 5i!*2o! + a/4 


Id 


= -52!*2o!/2 + a 


Sd 


= 5i!*2o! -a/4 


(1) 

c 


= 5i2!^2o!-3a/2 


(2) 
c 


= 25i2!*2o!-3a/2 



(6.42) 



are all finite for !*io! = and fc = 0. 

a. At the ASF~MSF critical point, !*io! = 0; On 
the ASF~MSF transition line, !*io! = (and also for 
a vanishing Feshbach resonance coupling, a ~ 0, when 
the order parameter phases are decoupled), the zero mo- 
mentum coefficients — the first four lines of (|6.41[) — vanish 
identically and one obtains two gapless spectra 



"'*^2 _ s^^ ± ^/Sdl + S(^^\ 



(Kf) 



which lead to two acoustic critical modes, i^J"* 



(6.43) 



hcT^k 



at small k. For !*io! = 
sound speeds are given by 



0, (5c^ , fci. vanish and the 



,crit\2 



(cT) 



!*2 



2?Tii 



-(7e±7rf) 



-k '' 



;^l*2o!, 
^^!*2o!^ (T = 2. 



cr = 1 



2mi 



(6.44) 



b. In the ASF phase, !*io! > 0; As found above, 
Eq. (|6.40p , in the ASF phase the spectrum of out-of-phase 
excitations, labeled by cr = 1, is gapped, while that for 
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a = 2 excitations, corresponding to in-phase fluctuations 
of the two condensates is given by 



ASFn2 



/pASi-'N 
{^k2 ) 



J1)J2) 
''k ^k 



-^kl 



(6.45) 



1 
2^ 



(eo + do) (5ek - 5di^ 



I* 



101 



21*2 



■do 



(2) 



+ (eo - do) ( ^ek + Sd^ + -r-j- — -Scl. 



'101 



in which (|6.38p has been used. As expected from the 
general symmetry arguments discussed in Sec. IIIII and 
earher in this section, the in-phase excitations are acous- 
tic, E^2^ « hc2^^k at small k, with sound speed given 

by 



(c: 



ASFn2 
2 ) 



I* 



201 



Aeo+do)f^'^ + {eo-do)£^ 



ieomi 



(6.46) 



with the constants /q defined by 



i'y) 



/^'^ ^ 52 1*20 1 



+ (5i2|*2o|-a/4) 



'lOl 



'201 



/(') EE 2(7121*201 



a-(-2.gi|*2 



'lOl 



'201 



(6.47) 



Note in passing, that, as expected, for a vanishing Fes- 
hbach resonance coupling, a ~ 0, both in-phase and out- 
phase modes become acoustic, with sound speeds 



2micl^^ 7r52|*20p+5i|*ioP 



(6.48) 




7:52|*20p-.gi|*i 



2.g22|*io|2|*2o| 



that are real and positive for gig2 > <?i2- 

c. Scaling form for small k and |*io|.' It is easy to 
check that the |*io| -^ limit of ()6.46p is very different 
from (|6.44p . which therefore does not commute with the 
k —f limit. In order to show more carefully the distinc- 
tion between these two limits, a scaling form that is valid 
when k, |*io| are both small, but have arbitrary ratio, 
is derived. By keeping only leading terms in |*ioP and 
eici, one obtains 



EtF = 2a|* 



^kCT 



101 



ge + ley 



(6.49) 



± Vt+7^y?+^ffF^2^I^Mff^ 



in which the dimensionless scaling variable is 

|*2o|£ki 



y 



'lOl 



(6.50) 



It is easily checked that for large y (I6.44p is recovered, 
while for small y (I6.46P is recovered. 



C. MSF paired ground-state wave function 

The zero temperature molecular superfluid ground 
state is constructed by requiring that it be the quasi- 
particle vacuum: 

7ka|MSF) = 0, for aU k^O,a. (6.51) 

The additional constraint 

ao2|MSF) = v1^*2o|MSF), (6.52) 

where ao2 — V^^-'^^^ /(in/'2(r), ensures that the MSF is 
a coherent state for the lowest single particle trap state 
k = and thereby has the correct amplitude *20 corre- 
sponding to molecular superfluid order. 
Using the commutation relations 



[a,e^"^] = Ae^"\ \a,e 



^a+bti ;^^teAatSt 



(6.53) 



where a, b are any two independent harmonic oscillator 
operators, it follows that the state 



|MSF) = exp f *2oVv42 - J E 



Xka4^a^_k<T ) |0) 



k^O.CT 



indeed obeys (|6.51[) and (|6.52p with the choice 



Xka 



Vka _ £kcr — Ey.^ 



Uko 



A. 



(6.54) 



(6.55) 



The factor of 1/2 in front of the sum is required because 
each term actually appears twice, once for k and once for 
-k. 

The quantity Xko- rnay be identified as the Fourier 
transform of the atomic {a = 1) and molecular {a — 2) 
pair wavefunctions with zero center of mass momentum. 
The asymptotic long-distance behavior of its Fourier 
transform Xa(r), which is now computed, is governed by 
the singularity of Xka nearest fc = 0. Since Xko- depends 
only on the magnitude k, one may use the Bessel function 
identity 



dQ^e 



:k-r 



27r'*/2 



fcr 



J d-2 (kr) 



(6.56) 



to perform the d — 1-dimensional angular integration in 
its Fourier transform, yielding 



Xa{r) 



dk 



„iU.-i 



Xka, 



(27r)'' 

°° 2k'^-^dk / 2 , . ^ , 

r^ka( — ) Ji^ikr). (6.57) 



(477)^^/2 



kr 



Since the right hand side of (|6.56p is an even function of 
fc, the integration may be extended to the full real line, 
avoiding the branch cut along fc < by shifting the con- 
tour an infinitesimal distance into the upper half plane, 
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iir(d-l) 



and simultaneously dividing by the factor 1 + e 
Since fc''~^(2/A:r)('*~^)/^ J_((;_2)/2(fcr) is analytic through 
the origin, and an odd function of k, its integral vanishes, 
and one may write Xa (r) in the form 



Xa{r) 



oo+iri j,d-l 



dk 



rXko 



H^^iikr) 



(6.58) 
in which 77 is a positive infinitesimal and Hi, (x) is a 
Hankel function of the first kind. 

From (|6.27p , one observes that E^J'^ has finite branch 
cuts along the imaginary k = |k| axis over the inter- 
vals ±i(fccr-, fc(T+), where ka-± = yflme^/h. To evaluate 
Xa- (r) the integration contour is deformed into the upper 
half plane to run down, around the origin, and then back 
up the imaginary axis, avoiding the upper branch cut. 
Since the Hi (x) decays exponentially in the upper-half 
plane, one can close the contour and then shrink it around 
the upper branch cut of E^^^ . Because the integrand is 
finite near the branch points, the infinitesimal circular 
parts of the contour integral, and the complete integrals 



of the analytic parts of Xko-i both vanish. The remaining 
parts on the left and right sides running along the branch 
cut double up, giving 



X<T(r) 



7r(27r)'^/2A^TO, 



2-d 
-f 2 



dKK'^/'^K;l-2{t 






X ^{K^~kl_){kl^~K^) (6.59) 



with Ki,{z) = {■Ki/2)e''''''^H'i^\iz) the modified Bessel 
function. For large ka+r the integral is dominated by the 
region near k^-, and one may safely (with exponential 
accuracy) extend the upper limit to infinity and approx- 
imate the square root factor by the form 



(6.60) 
the lower relation being especially required for a = 2 
where fc2- = 0. One therefore obtains 



J 



Xa(r) 



2X„n 






2A„m„ 



d+l ^d+1 5 



d/2 



d-f2 J rZcr—T ^ i 



(6.61) 



r 



in which the asymptotic form Ki^{x) w ■\/7r/2xe ^, \x\ 3> 
1, has been used to obtain the first line, and the identity''* 



f 

Jo 



x-ir.(x)d. = 2--ir f i±4±^ V ^' + '^ " " 



(6.62) 
to obtain the second. 

It thus follows that in the MSF phase the relative 
atomic wavefunction decays exponentially according to 

Xi(r) ~ e^*"/^", with a decay length 



ft-^ 



/ci_ -y/2TOiei_ ■ 



(6.63) 



reflecting the confinement of (gapped) atomic excita- 
tions, and the corresponding absence of atomic long- 
range order inside the MSF. Since ei_ ^ h' — i^c, ^1 ^ 
(y — v^~^l'^ has a square root divergence as the ASF 
phase boundary is approached. On the other hand, since 
fc2- — 0, the molecular wavefunction X2(r) ^ Xjr'^^^ has 
a power law decay, reflecting the existence of molecular 
long-range order inside the MSF. 

Note that the ground state (|6.54p . in addition to be- 
ing a molecular coherent state, is also an (atomic and 
molecular) 'pair coherent state. It thus makes explicit 
that within the molecular superfluid state, a molecular 
k = condensation, ^20 7^ Oi is accompanied by a 



nonzero BCS-like atomic pairing at finite relative k, with 
an anomalous correlation function, 



avp 



(akifl-ki; = -WkiWki = 



20 



2-Eki 



(6.64) 



Exactly the same branch cut structure as described above 
applies to the right hand side of (|6.64p . and its Fourier 
transform, the BCS-type atom pair correlation function, 
falls off exponentially at the same rate e"''/'^^ . The cor- 
relation length ^1 (that is finite inside the MSF, but di- 
verges as the transition into ASF is approached) char- 
acterizes the size of the virtual cloud of atom pairs sur- 
rounding each closed-channel molecule (whose size, doj 
characterized by the microscopic range of the interatomic 
potential, remains finite throughout). 

On the other hand, the molecular anomalous pair cor- 
relation function 



(ak2a-k2) 



-u\av\a 



92^ln 
2Ei,2 ■ 



(6.65) 



exhibits a 1/fc divergence near the origin [on top of the 
^2o^(27r)''5(k) condensate contribution due to the long- 
range order] , so that its Fourier transform approaches the 

^ power law decay. This 



^20 asymptote via a slow l/r"^ 



is a signature of quantum fluctuations in the low energy 
molecular Goldstone mode. 
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D. Thermodynamics 



MSF phase 



As is clear from (|6.9p and (|6.24[) . within the Bogoh- 
ubov approximation a superfluid (be it MSF or ASF) 
is a coherent state with excitations described by a gas 
of noninteracting bosonic Bogohubov quasiparticles, 7kCT, 
respectively given by (j6.6p and (|6.26p . Thermodynamics 
is therefore easily computed in a standard way. 



The free energy density in the MSF consists of the 
ground state condensate energy Tirai = '^[^20]: plus a 
contribution from the noninteracting Bogoliubov quasi- 
particles, governed by dHyiSF: Eq. (|6.9p . A standard free 
boson computation gives 



/mSf[/J2] = -M2I* 



20 I 



7^321*20!'' -Rc[/7^;*2o] 



E 



{2tj:Y 



ilnd 



-f}El 






(6.66) 



r 



where a complex molecular "source field" hi (that van- 
ishes for a physical system) has been included. As usual, 
derivatives of /msf[^2] with respect to /12 generate cor- 
relation functions of the molecular field. In interpret- 
ing this quantity, it is important to emphasize that ^20 
here (in an unfortunate abuse of notation) is the mean 
field order parameter, an explicit function of the Hamil- 
tonian parameters fi, v, hi, etc., that docs not include any 
fluctuation correctionsi^l The leading Bogoliubov correc- 
tions are provided by the hi derivatives of fyiSFlhi]. For 
the molecular condensate order parameter, corrected by 
quantum and thermal fluctuations this gives: 



*9 = * 



20 



5'^') 



a/] 



MSF 



dh% 



(6.67) 



'12=0 



in which hi enters through its explicit appearance in 
the first line of (|6.66p as well as implicitly through ^20- 
The extremum property of ^20 with respect to i^mf [^20] 
therefore gives 



5*2 = 



<9/msf \ 



9* 



20 



dh% 



11-2=0 



*20 
2^11 



[Id,i{n,iy) + 2Id,i{f^,j^)], (6.68) 



where the mean field longitudinal susceptibility is 

idno/dh*i)h2=a = l/(-2M2 + 6g2|*20p) = 1/4^^2- Con- 
sistency requires that the original forms (|6.5[) be used for 
the '520 dependence, and 



id.i 



Id,2 — 



dk 

^2^ 
rfk 

{2ttY 



n\^i 



nk2 



1\ 2ffi2£ki ~ a^ _ 

2) i^MSF 

1^| 2g2£k2-g||*2 
^MSF 
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(6.69) 



where riko- = {e^-^'^" — 1)~^ are the standard Bose occupa- 
tion factors for Bogoliubov quasiparticles. The number 
density to this same order is 



5/1 



MSF 



„ , =2|*2p + <^riMSF(T,i/), (6.70) 

where the density of bosons not condensed into the lowest 
k = single particle state (i.e., the condensate depletion) 
is given by 



SnusFiT, v)=^(J j ——J [vl„ + (uL + f'D^kff] . 

(6.71) 
The depletion density (J^msf {T, v) comes from the ex- 
plicit /i-dependence in E^^^ and remains finite even at 
zero temperature due to the interaction-induced zero- 
point contribution v^^i^ The remaining implicit jjl- 
dependence entering through the condensate ^20 gives 
rise to the term |*2p in (|6.70p . in place of the mean field 
condensate density |*20p- 



Evaluating (|6.7ip at T = and v = Vc one obtains 
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fcMSF(0,.e) = E|/(|^ 



Eko- +£(T+/2 



- 1 



in which ei+{i'c) — 2a|'I'2o|, £2+{vc) = 2(72|'J'20pi and the coefficient is given by 



Bri f2^mY'' ^^a,, ^ 2(,+2)/2^^|2) (g 72^ 



2r(d/2) \ K^ J ^^'+ 



B, 



d«z;(^-2)/2 



i; + l/2 



' t('-1\v^'-' 



d^/n 



(6.73) 



In d = 3 one finds B^ = 1/3 and Eq. p.l4p quoted in the Introduction immediately follows. Since Sj^^ ex n2Q , this 
"correction" term becomes much larger than n2o close to the MSF-N transition line. This is a sign of the breakdown 
of the mean field description of criticality, and (|6.73p ceases to valid in this nontrivial critical regimei^ 
For 1/ < I'c (i-e., inside MSF phase) one obtains: 



in which 5 — £\^(y) j e\j^{v) , and 



Bd ( lixm \ 



d/2 



2r(d/2) V h'^ 



oo ^(d-2)/2 



j {[l + bM]e{': + 2^'^''"4':} 



v + {l + 5)/2 t; + l/2 

^{v + 5){v + l) ^v{v + 1) 

I 



(6.74) 



(6.75) 



Of interest is the behavior of this integral near v^.^ i.e., 
for small 5. The singular behavior can be obtained by 
first computing the derivative 



dS 



l-S 



vid~^)/^dv 



4Brf Jo {v + S)y^V^7TT 

For d < 3 the integral diverges as (5 ^ 0, and one obtains 
dbd l-S 



d5 



4.B, 



-[l3d,sS^''-^^/^ + (3d,i + 0{S)], (6.77) 



where the singular coefficient f3d,s is obtained from the 
small V part of the (infrared divergent) integral by scaling 
out 6 via the change of variable u = v/S:'^^ 



Pd,: 



°° u('^-2)/2^^ 



2 /d 



3-d 



/o (7/ + 1)3/2 ^- \2J' \ 2 

(6.78) 
The linear term is obtained by first subtracting this small 
V singular (in S) part of the integral, and then letting 



Pd,i = 



\Vv + l 



-1 



1 



=r 



3\^A4- 



(6.79) 



On the other hand, for 3 < d < 4, dbd{0)/dS is finite^^ 
and one finds the leading term (3d,i simply by setting 6 = 
0. Related to this, the singular term no longer diverges, 
and it is obtained by first subtracting the (3d^i (the (5 = 0) 
term, and then again simply scaling 6 out of the integral. 



One may verify that the final results for both coefficients 
are identical to (p?75)l and ^7^. Integrating (pTT)) with 
respect to 5, one finally obtains 



(6.76) bd{S) 



1 
45^ 



d-1 



Pd,sS^'-'^/' + Pd^S 



[1 + 0{S)]. 



(6.80) 

In d = 3 both (3d.s and Pd,i separately diverge. However 
the sum is finite, giving rise to a logarithmic dependence 
on S: 

baiS) = -^S{\ii{l/S) + 41n(2) - 1} [1 + 0{d)]. (6.81) 

This same result also follows from a direct asymptotic 
evaluation of the integral (|6.76p in d = 3. 

Defining a T = critical exponent a via SriMSpiT — 
0,(5) ^ S^^"' (i.e., the zero-temperature quantum transi- 
tion analog of a specific heat exponent), one finds 



^, d^3 
(log), d = 3 



(6.82) 



This result will be modified by critical fiuctuations 
sufficiently close to the MSF-ASF quantum phase 
transitionsi^ The resulting behavior of the condensate 
depletion SnyiSFiT,!^) is illustrated in Fig.[31 

Before ending this subsection, the order of magnitude 
of the MSF zero-temperature depletion (|6.74p is exam- 
ined in light of the identification in Sec. [ll] of the small 
parameter 7, Eq. (j2.14p . In order for the fluctuation 
correction (|6.74p to be accurate, it is necessary that it be 
much smaller than the mean field value n2o (and the same 



27 



should be true of the correction S'i'2 relative to ^20) 
ing forms (|6.28p for the energy gaps, one obtains 



Us- 



n2o 






ex 7' 



d(4-d)/4 



(6.83) 



{2me2+/n^y/^ 



"20 






(6.84) 



in which 32 oc {h'^/2m)a'^~^ and a ex (;iV2m)|ro|(''-'')/2 
relate the Hamiltonian parameters to the molecular scat- 
tering length and effective range [see (|2.7p and (|2.8p ] in 
rf dimensions. In (j6.83p it has been assumed that ei+ is 
of the same order of magnitude as it is on the MSF-ASF 
phase boundary, where ei+ — 2ay/n2o- 

It is seen that the two terms in (|6.74p are very different 
in character. The second term, estimated via (|6.84p . is 
the standard result for a monatomic Bose gas, and is (for 
d > 2) small in the dilute hmit, n^/'^a2 < 1. The first 
term, estimated via (|6.83p . is small (for d < 4) only if 
7^1. However, this requires n^''^rQ >> 1, which places 
a lower bound on the density. The expansion about mean 
field theory presented in this paper therefore requires a 
sufficiently narrow Feshbach resonance (small a) such 
that the separation of scales tq :^ 02 exists, and its va- 
lidity is limited to densities in the intermediate regime 



1 



1 






(6.85) 



This confirms, within an explicit perturbation calcula- 
tion, the claims made in Sec. Illl 



2. ASF phase 



The Bogoliubov corrections to the mean field order pa- 
rameter are then found in the form 



with 



S^a 



>r- / a/ASF \ /avI','o\ ,„„„. 



in which, due to the mean field conditions ()4.4p . only 
the non-mean field part of /asf actually contributes to 

Self-consistency, via the mean field equations, but with 
hi/2, /12/2 replacing the zeroes on the left hand sides 
of (|4.4ap . (|4.4bp . respectively, determines the h^r and fi- 
depcndcnce of ^o-q. The four (complex) equations for 
{d'^a'o/dha-)h^=o dccouplc into a pair of separate equa- 
tions for id'Sao/dhi)hi=h2=o and (9*cto/9/i2)/ii=/j2=o: 



{d^2o/dK)h„^o 



A B 
C D 



(6.89) 



where 
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D 
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2*10 (2ffi2*20 - a) 
2*io(2gi2*io - a) 
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*2n 



(6.90) 



and ID = 



Jj and|2)=(^5 

The free energy derivatives in (|6.88p , taken at constant 
values of the Hamiltonian parameters fi^, ga,ct, gi2, are 
given by 



Next consider the i^ > Vc case where both vJ/^q ^ 0. 
Computations in the ASF phase are most conveniently 
performed by taking hcrQ and ^o-o real and positive at the 
outset. Expressions for thermodynamic quantities are 
quite long and involved, and they will only be sketched 
here. 

The Bogoliubov free energy density in the ASF is given 
by 



/asf = 



^ f dk 



(27r)° 



■ln(l 



-PEi 
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( pASF - N 



(6.86) 



where 7imf[^iOi ^20] takes the form of (14. ip . but with, as 
in (|6.66p . additional ordering field terms — J^a R-^i^^o-^o-o] 
now included. 
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2 9*^0. 
(6.91) 

with the two energies, and the parameters entering them, 
given by ^^, ((OT|) . and (|0^ . These will not be eval- 
uated any further here, except to note that, as in the 
MSF phase, at zero temperature the leading behavior of 
the integrand at small k is proportional to 1/Eu.aj while 
at finite temperature it is proportional to ^/E^^. As ex- 
pected, the acoustic mode therefore generates divergent 
fluctuation corrections for d < 1 at T = 0, and for d < 2 
for T > 041^. 

The total density is given by 



dfASF 

-1-2*10(5*1 



T,v 



2m 



4*90(5 'J'? 
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SnASFiT,iy) 
(6.92) 
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where, in the second Hne, the two terms Unear in S'^a 
subsume the implicit dependence of ^o-o on /it. This re- 
sult follows from the fact that, via (|4.4p . d'^ao/dji obeys 
^^Wji . but with 2*io|l) +4*2o|2) replacing \a) on the 
right hand side. The depletion 5nxa^Y{T,v) may be de- 
rived either from the derivative of the non-mean field 
part of (j6.86p with respect to the explicit ^-dependence 
(which appears only additively in £ko-) — yielding a form 
identical to the right hand side of (|6.9ip . but with d/dfi 



(performed at constant ^cro) replacing d/d'^ao — or as 
the total number of uncondcnscd particles. 



SnASFiT,iy) = y^<^{ai 



WaO'ka) 



k,CT 



(6.93) 



in which (j6.26p connects the d and 7 operators. Using 
either approach, one obtains 
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di 
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(6.94) 



The terms involving ci 



are of order t{,t\,tit2. 



and 



therefore vanish along the MSF-ASF transition line. Be- 
cause these expressions involve a number of parameters, 
such as scattering lengths and the Feshbach resonance 
coupling, their final integrated expressions are not very 
enlightening without additional (e.g., experimental) in- 
put. Thus these predictions are not explicitly evaluated 
further. It is noted only that estimates similar to ()6.85|) 
(but now involving all three scattering lengths on the 
right hand side) may be derived for the range of validity 
of (I6J4I) . 



E. Superfluid density 

The superfluid (number) density Ug is a measure of 
the stiffness T^ of the order parameter against a long- 
wavelength spatial gradient in its phase 0(r), defined by 
the corresponding change in the free energy 



1 



AF, = -T,(r) / dr|V0| 



(6.95) 



Expressing the free energy in terms of the superfluid ve- 
locity Vg = {h/m)Wd, one obtains AF^ = ^Ug Jdrjvsp 
with the standard relatior 



,83 



ft2 



(6.96) 



For a two-component Bose (atomic and molecular) 
gas that is considered here, at long length scales the 
two phases are locked by the Feshbach coupling to be 
02 = 29i, i.e., 9a = a9. As discussed previously, this is 
obvious in the ASF state, where the gas is a superfluid 
with respect to both atoms and molecules and out-of- 
phase fluctuations 9i — 62/2 are gapped. It is also valid 
more generally, coming from the requirement that the im- 
posed velocities of atoms {a = 1) and molecules {a = 2), 



VsCT = {fi/ma)'\79a, relative to a stationary boundary are 
the same. 

As outlined above Us is calculated by computing the 
free energy change AFs{ko) in the presence of a uniform 
phase gradient ^o(r) = ko • r, corresponding to a super- 
fluid component with uniform velocity vq = hliQ/m (and 
stationary normal component). To this end one imposes 
phase twist boundary conditions on the field operators: 

Mr + Ln)^e'^''"Mr), (6.97) 

where L is the system length along a chosen direction n 
of the phase gradient and ^o = ^o(-^n) — ko • hL. From 
this definition, one obtains^^ 



2^2 



dkl 



(6.98) 



where fe^ is the free energy density in the presence of the 
twist ^0: a,nd L —f 00 includes the thermodynamic limit 
and is to be taken here at fixed ^q. 

Reexpressing the Hamiltonian in terms of periodic field 



operators ipa - 



■cr'k.o-r 



jpa, one obtains 
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2mi 



vo-P, 



(6.99) 



where 



Pa = ~ih dr^i(r)VV;<,(r) 
Na = /dr^i(r)^,(r) 



(6.100) 



are momentum and number operators for component a 
andP = E.P.- 
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From the above form for ff, and defining equations 
(|6.96p . (|6.98p for 12^ and T^, one observes that the super- 
fluid density is also given by 



djs 



dihko) 



ko=0 



or equivalently defined by the relation 



lim j^ = TigVo, 



(6.101) 



(6.102) 



where the supercurrent density js is the expectation value 
of the (number) current density operator 



1 dH 
V dhko ' 

V \ m 



(6.103) 



To compute Ug one expands Hlip^] to quadratic order 
in the fluctuations (j)^ — ips — ^ctO and diagonalizes it 
at a finite /cq. Because it is odd under k — > — k the 
new momentum term remains diagonal under the ko = 
Bogoliubov transformationjil 



Pa = ^ hkal^Uka = Y^ ftk7^^7k<T- 



k 



(6.104) 



Thus the BogoUubov transformation at finite kg is un- 
changed from that at fco = 0, except for a shift in the 
chemical potential fi -^ jl — fi — fi?k'^/2m. The spec- 
trum, however, does change, but in a simple way 



-Bko- = £^kcr + vo • ?ik, 



(6.105) 



that is in accord with a general requirement for a 
Galilean-invariant system. Computing the expectation 
value of js in (|6.103p and using (|6.10ip [or, equivalently 
computing the free energy and using (|6.98p ] one finds 



ns{T) = n- n.n{T) 

2 v-^ /■ dk 



dn 



ko- 



dE'ko 



(6.106) 



where clearly rig < n, i.e., the normal fluid density 
n„ > 0. At zero temperature all excitations are expo- 
nentially suppressed and the normal fluid density n„ van- 
ishes, giving ns{T — 0) = n independent of interactions, 
as required by Galilean invariance. In the normal phase, 
where -Eko- = Sko- — Mo-, an integration by parts yields 
n„ = n, and rig vanishes as expected for a normal fluid. 
In contrast, in the MSF phase, despite the absence of 
atomic long-range order, at finite temperature there is a 
nontrivial atomic contribution (a = 1) to the superfluid 
density (though not to the condensate) . The correspond- 
ing reduction in n^ is due to thermally excited, unpaired 
atoms with a gapped spectrum E^^^ (due to Feshbach 
coupling to condensed molecules), that is not simply the 



free spectrum Ski of the normal state. In the weakly 
interacting limit, away from both T — and T = Tc, 
ns{T)/n w 1 - {TlTcY'"^ is well approximated by the 
ideal gas form, which in turn coincides with the con- 
densate fraction n2o{T)/n. As usual, sufficiently close 
Tc this result must be modified by critical fluctuations 
which strongly modify (|6.106P ]^^i^^ On the other hand, 
deviations near T = 0, where the reduction in Us is dom- 
inated by gapless molecular excitations that are sound- 
like with E^^^ ~ \/£fc2£2+) are accurately described by 
(|6.106p . which implies that n„(T) ~ T'^+^. The low tem- 
perature crossover from T'^'^^ to T''/^ takes place when 
the temperature is high enough that excitation of the 
higher energy quasiparticles with quadratic dispersion, 
Ewa ~ £kcr, dominate the thermodynamics.— For spe- 
cific model parameters, the full detailed form of UgiT) 
can be straightforwardly evaluated numerically. 



VII. MSF-ASF PHASE TRANSITION 

As has already been seen in Sec. IIIIl many of the prop- 
erties of the phase transitions appearing in the phase di- 
agram. Fig. [5] can be deduced based on the nature of 
the underlying symmetry that is spontaneously broken 
in the MSF and ASF phases. In particular, there it was 
argued that because the MSF exhibits a discrete residual 
■01 — > —01 (global phase rotation by tt) symmetry associ- 
ated with the diatomic nature of the molecule, the MSF- 
ASF transition at the level of mean field theory is of the 
Ising type.-®^ However, it is important to stress that the 
Z2 symmetry that is broken at the MSF-ASF transition 
does not necessarily imply that the critical properties of 
the transition (beyond a mean field approximation) are 
of Ising type. One condition for this is that, in addition 
to the local interactions of Ising symmetry, the derivative 
terms associated with kinetic energy (spatial gradients) 
and (in quantum theory) the Berry phase ("|0p9t6'" time 
derivative) terms must reduce to a standard rotationally- 
invariant (in Euclidean space) d -\- 1 dimensional gradi- 
ent term. The other condition is that additional fields 
(if any) coupled to the Ising order parameter must not 
modify the Ising critical behavior, i.e., must be irrelevant 
in the renormalization group sense. 

Below these issues are explored in more detail. It will 
shown that although the first condition is indeed sat- 
isfied (i.e., for the quantum MSF-ASF transition the 
scalar Ising order parameter indeed has a standard gradi- 
ent d + 1-dimensional Lorentz-invariant "elasticity" ) , the 
existence of the Goldstone mode (the phase of the molec- 
ular order parameter) that couples to the Ising field can 
have nontrivial effects and (as was first pointed out by 
Lee and Lee^^ based on an earlier study by Frey and Ba- 
lents in a different context^) likely drives the MSF-ASF 
transition first order. For the extremely dilute gases of 
experimental interest, this first order behavior is weak 
and may only be visible very close to the transition. 

Focusing on a homogeneous trap (a box), the T — Q 
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and finite T MSF-ASF transitions will now be studied in 
more detail. This can be most easily done working with 
the coherent-state action, S, Eq. (|2.17p . corresponding 
to the two-channel Hamiltonian (PTTj) . The MSF-ASF 
transition will be studied from the MSF side, where the 
molecular field ^^2 = |V'2|e'^^ exhibits massless Goldstone 
mode phase fiuctuations in 02, and small, gapped fluc- 
tuations in the magnitude |V'2| about the molecular con- 
densate (■02) = ^20- Integrating out the latter leads to a 
standard superfluid hydrodynamic action^ 



<S'2[6'2 



] = \^sj'^ ^drjdr [c^lAdr92f + (^9,)'] 



(7.1) 



that controls the (acoustic) fluctuations of 02, with 
sound speed cmsf given by (|6.29p . and helicity modu- 
lus/superfiuid density Tg given by (|6.96p and (|6.106p . 



The atomic contribution to the action, together with 
the key Feshbach resonant atom-molecule coupling, is 
given by 



^1 



f3h 



dr / dr 



^{hdri^i - ^1 ( ^V^ +^1)^-1- aRe (|*2| 



'V'lV'i) 



s, 



nonlinear; 



(7.2) 



in which S'noniinear contains the conventional quartic scattering terms. As discussed earlier, the latter locks the molecule 
and atom phase fluctuations such that low energy excitations are governed by 9i = ^2/2. Thus it is convenient to 
define "dressed" atomic fields V'l according to 



^] 



= „ie2/2 



which leads to 



(3h 



1 ~ n^ 



Si= I dr I dr^ ri^dr^i + i^\i^ivme2 + —ri 



■01, 



_iV + -V6'2 



(7.3) 



Ml 



V'l - aRe(|*2|V'lV'l) \ + SnonM 



(7.4) 

Straightforward analysis^^ shows that near the MSF-ASF transition the minimal coupling to the induced gauge- 
like field V02 above is irrelevant near a Gaussian fixed point. Dropping this subdominant contribution and writing 
"01 = '/'-R + ^V'/ in terms of its real and imaginary parts one finds 



Si^ dr I dr 



^ii'% + V'l)^9,02 - 2i^iWr^R - 4>R (^V2 + flRJ tjjR - ^Pj i^V + m ] ^I 



S, 



nonlinear 



(7.5) 



where the shifted chemical potentials are given by 
//fl = At + 2a|*2|, /i/ = M-2a|*2|. 



(7.6) 



This form of Si makes it clear that in the presence of the 
molecular condensate, |^2| > 0, positive a reduces the 
0(2) = C/(l) symmetry down to Z2, and with im > /// 
results in ipR reaching criticality before 4>i- Because the 
canonically conjugate field tpj remains "massive" (non- 
critical) at the MSF-ASF critical point [defined by where 
the coefficient fiR oi^j^ vanishes, consistent with (jl.lOp ]. 
it can be safely integrated out and leads to a d -f- 1- 
dimensional (Lorentz-invariant) action which is even in 
the scalar order parameter (f> = ipR, and whose relevant 
part is given by 



Ses[02, 0] = Ssf[02] + SlsinM + ^int [^2 , 



(7.7) 



i./3h 



in which 

SsF[e2] - 2^^ J '^^ j '^'^ V'^sAdre2? + (V02)'] 
-Slsing [(!>] = 



ph 



^KAdr(t>f 



dr I dr 



fi2 



SnA02,<l>] = 



/3h 



dr / drcj) hdr02, 



(7.8) 



represent separate superfluid hydrodynamic and Ising ac- 
tions, together with a Berry phase-like term that cou- 
ples them. The coefficient K^ « l/(a|\E'2|) to lowest 
order in 1/a, and the leading (f>* nonlinearity comes from 

'-'nonlii 



^nonlinear- 
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Thus as advertised, if the coupHng of the Ising order 
parameter </> to the molecular Goldstone mode 62 is ne- 
glected, the T = MSF-ASF transition (near ^ir = 0) is 
indeed in the {d-\- l)-dimcnsional Ising universality class; 
at finite T it crosses over to the d-dimensional Ising tran- 
sition. The Ising transition is well studied, and leads to 
the following predictions;— At T = 0, for d = 3, up 
to logarithmic corrections, the mean field theory derived 
above remains an accurate description. On the other 
hand in d = 2, the MSF-ASF exponents are nontrivial 
but are well-known. For example, standard scaling argu- 
ments predict: 



nw 



\2Pi 



gap 



(7.9) 



where /?/ « 0.31, zj — 1, and I'j w 0.63 are 3d clas- 
sical Ising exponents. These, together with the rel- 
evance of T at this quantum critical point, also im- 
ply a universal shape of the MSF-ASF phase boundary 
Vc{n,T) - Vc{n,Q) + aT^I"' , as shown in Fig. H One 
may hope that when long-lived molecular condensates 
are produced, nontrivial behavior of EgJ^iv) and the full 
excitation spectra may be observed in Ramsey fringes^, 
and in Bragg and RF spectroscopy experiment o^^'^'^i^'^'^^ . 

However, as first emphasized and studied by Lee 
and Lee^i^ the existence of O2 fluctuations can mod- 
ify this conclusion sufficiently close to the MSF-ASF 
transition — intuitively this follows from the fact that if 
one integrates out the superfluid fluctuations, a long- 
range power law cjp'-cjp' interaction (highly anistropic in 
space-imaginary time) is generated. Indeed, at T = and 
d+ 1 < 4 the (p^drd2 coupling term is relevant around the 
Gaussian fixed point, scaling like 6'-^"'^^/^ with increasing 
renormalization length scale b, and thus competes with 
the Ising cj) nonlinearity. The resulting theory embod- 
ied in S'cff has, in fact, a form very similar to that of an 
Ising model on a compressible lattice^ with 62 playing 
the role of a phonon u and the 62 — <j) coupling scaling 
similarly to the magneto-elastic coupling (fP'V ■ u.^^ The 
latter model (as well as its Heisenberg generalizations) 
have been extensively studied. ^i The conclusion of that 
work is that for a/ > (with a/ the specific heat expo- 
nent of the uncoupled model) the magneto-elastic cou- 
pling leads to runaway flows that has traditionally been 
interpreted as a signature of a fluctuation-driven first or- 
der transition. ^^ For a/ < Goldstone mode fluctuations 
(lattice elasticity) is in fact irrelevant and the transition 
is in the universality class of the usual (elastically rigid) 
Ising model. 

Based on these results, since the d-dimensional Ising 
specific heat exponent is positive for d > 2^^^ one thus 
concludes that here too, sufficiently close to the MSF- 
ASF critical point, the transition is driven first order. It 
should be emphasized that many results in this paper, 
namely those that refer to thermodynamic and elemen- 
tary excitation properties of the different phases, not too 
close to the transition lines, remain valid and are unaf- 
fected in any way by this issue. 



VIII. BOSE-BCS MODEL 

As is shown in this section, the analysis of the Bose 
atom-molecule system via the two-channel model, pre- 
sented in the previous subsections can be complemented 
by Bose-BCS variational approach of the one-channel 
model (|2.15p . Similar analyses have been presented 
previously,^-'^^ but it is worth presenting, and general- 
izing them somewhat, here in a form that can be com- 
pared to the results of the two-species model (|2.ip that 
has so far been the focus of this paper. A description 
entirely in terms of Bose atomic constituents (single- 
channel model) lends further physical insight into the 
microscopic nature of the phases and phase transitions, 
and facilitates comparisons with BEC-BCS crossover in 
Fermi systems i^i^ii^ii^ii^iiiii^ Furthermore, the atom- 
only model should be more appropriate for describing 
the case of a broad resonance. 

The phase diagram of the system is explored as a func- 
tion of the gas parameter n^/^a, where n is the atomic 
density and a is atomic scattering length. As discussed in 
Sec. ini experimentally a is controlled by magnetic field- 
tuned proximity to the Feshbach resonance (diverging at 
the resonance) as well as by the atom specific background 
scattering length, a^g- However, within the one-channel 
model (|2.15p . this is encoded into the tunable pseudo- 
potential amplitude gi, characterizing atom- atom micro- 
scopic interaction and connected to experiments via the 
atomic scattering length, 0(51)1^ Focusing on a dilute 
gas, the gas parameter away from the Feshbach resonance 
(where a — > Obg ~ do, with do the microscopic range of 
the interatomic potential) is taken to be small. 

For detuning far below the resonance, the attractive 
interaction is strong enough to lead to a deep two-body 
molecular bound state, that corresponds to the appear- 
ance of a dilute gas of strongly bound compact molecules. 
In addition, for a large atomic scattering length the sys- 
tem is known to exhibit Efimov states of trimerS)^ leav- 
ing the questions of the stability and the nature of the 
condensed state open.^^ However, for a sufficiently large 
repulsive three-body interaction, and/or away from the 
Feshbach resonance (short scattering length), the system 
is expected to be stable. In this case, at low energies 
it will be governed by an effective low energy s-wave 
molecule- molecule scattering length am, corresponding 
to a repulsive interaction characterized by a molecu- 
lar pseudo-potential (72- Thus, the appropriate effective 
Hamiltonian is given by (j2.ip with all terms containing 
the atomic field ipi dropped. The theory is identical to 
an atomic theory of a dilute gas of composite bosons with 
mass ?7i2 = 2?7i, and chemical potential /i2 = 2/z — i?b, 
where i?b is the binding energy. Thus a Bogoliubov anal- 
ysis provides an essentially exact description. In par- 
ticular, at T = the system is a vacuum for ^2 < 
and Bose condensed for /i2 > with order parameter 
*20 = \J l^2lgi^ density ^2 = ^20 = |^20p and low 
energy acoustic excitation spectrum £'(k) = C2?ifc with 

C2 = 1/71252 /'712 = ?Ia//X2/to2- 
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In the opposite limit of repulsive interactions, no 
molecules are present and the system is, conversely, de- 
scribed by the Hamiltonian (|2.ip with all terms contain- 
ing the molecular field -02 set to zero. The phenomenol- 
ogy is again that of a dilute, single component (this time 
atomic) Bose gas as described above, with constituents 
of mass m, chemical potential /ii and interaction gi. 

The focus here is on the interesting intervening re- 
gion around the transition between these two atomic 
and molecular superfluid phases. Thus, the behavior 
of the single species fluid is explored from the conven- 
tional gi > (positive scattering length) atomic BEC 
limit, through the small gi < (such that the molecular 
bound state, or the resonance, is of spatial extent com- 
parable to intermolecular separation) Bose BCS limit, to 
the larger gi < (indicating a two-body bound state 
that is much smaller than the intermolecular separation) 
molecular BEC limit. As shown in Fig. [SI over this re- 
gion the atom-atom scattering length goes from positive, 
to negative, and then back to positive. 

The one-channel grand canonical Hamiltonian is given 

by 



H 



rfr|-^t(r)(^|!.V2+M)^(r) 



Ur.\3, 



+ -gijj\ry4,{ry + -wTp^ryi>{r) 



where for simplicity the subscript "1" on the atomic fields 
has been dropped. Though g can have either sign, and 
the new physics of primary interest here enters for g < 0, 
w remains positive to ensure thermodynamic stability. 
The relation of (18. 1|) to the more general two-channel 
model (that reduces to it in the wide resonance limit) 
was summarized in Sec. |lTl and discussed in detail in 
Ref. M. 



A. Variational mean field approximation 

In the dilute limit, and away from any phase tran- 
sitions, the variational approach to be presented, es- 
sentially equivalent to the Bogoliubov approximation, 
provides asymptotically exact results. In a standard 
treatments^ the approach relies on the inequality 



F < F„ = Fi 



MF 



(H-H, 



MF/MF, 



(8.2) 



between the true free energy F ~ — A;BTlnTr[e ^-^] 
(where H is the system's full interacting Hamiltonian) 
and the variational free energy Fy defined in terms of an 
arbitrary Hamiltonian Hmf and its corresponding free 
energy Fmf = ~kBT\nTr[e~P ""■"'']. Here (•)mf is the 
thermodynamic average with respect to Hmf- Since Fy 
is an upper-bound for F, Eq. (|8.2p guarantees that the 
better the choice of Hmf the closer one can approximate 
the true free energy with Fy. On the other hand to take 
advantage of the variational method one needs to pick 



a simple enough Hmf that the thermodynamic averages 
appearing in Fy may be calculated explicitly. Thus, Hmf 
is chosen here to be a quadratic Hamiltonian 



Hmf — 



E 



1 



A^kakak + -Fk(aka_k + H.c 



At At 



- ^QoWiao+al), 



(8.3) 



with variational parameters A'k, Pk, Qo (chosen real, by 
absorbing any extra phase factors into the Bose operators 
if necessary) to be selected to minimize Fy . 

The linear term (required to deal properly with the 
possibility of an atomic condensate) is removed via a zero 
momentum shift 



0(r) = ■0(r) - V'o, Ok = flk - V VV'o'^k.o 



with 



V'o 



Qo 



2{No + Po) 
Following this, the Bogoliubov transformation 



(8.4) 



(8.5) 



flk 

4 



Mk7k - I'kTlk 

Mk7k ~ «k7-k 



(8.6) 



with the choices 



Eu 



l+vt 



2 I Fk 



Nl - ^k^ 



■7) 



leads to the diagonal quadratic form 



Hmf — 



E 



k '- 



A,U 



1 



i?k7^7k+2(^k-i?k) 



+ [{No + Po)^Po ~ Qo^o]V. 



(8.8) 



The two-point averages are easily computed, and are 
given by 



iVk 



(4ak)MF = TT ( "k + ^ ) - ^ 



(8.9) 



/;;t;;t 



Pk 



k/MF — Vfl-kftk/MF — —-^T- I "k + X 



Eu 



1 



in which rik 



p/3iSi 



■^ — 1) ^ is again the Bose occupation 



factor. The atomic number density and molecular order 
parameter defined by Hmf are therefore given by riMF = 
-00 + "MF, $MF = V^o + "^MF, in which 



"-MF = (V' V')mF = 
^MF = (iplp) MF = - 



dk 

(27r)° 



iVk 
Eu 



dk Pk 

{2ttYK 



1 
nk+2 



nk+2 



(8.10) 
(8.11) 
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Using (|8.4p to shift the operators in _ff , together with 
Wick's theorem and Eqs. (|8.10p . (|8.1ip . the variational 
free energy density (|8.2|) takes the form 



/. 



//71 C 
— ^ < ksT [nk In(nk) - (1 + "k) In (1 + ?ik) 



1\ 1 



^k 



^o("MF,*MF,l/'o), 



.12) 



in which 



/o = -finuF + 2.9(2"-MF + ^mf) 



-wnMF(2nMF + 3$^ 



MF^ 



fl = fi- g{2nMF + $mf) 

3 , _ 

- -w(2nMF + 2nMF$MF + $mf) 

g = .g + u)(3nMF + 2I>mf)- 



(8.13) 



In minimizing (|8.12p one may treat the ratio Rk — 
Pk/Nk, nk, and -00 as independent variables. The deriva- 
tive with respect to i?k yields the simple result 



Rk — 



7mf 



Ck — Mmf 



.14) 



in which 



Mmf 



dfi 



MF 



7MF 



^ - 2gnMF - -w(2nMp + $Mp) 

- [2g + 3u;(2nMF + *mf)]V'o " ^wt^t , (8-15) 

a$MF 

= (g + 3ii;riMF)$MF 

+ [5 + 3w(nMF + $mf)]0o + ^t'V'o- (8-16) 

All k-dependence therefore resides in the energy denom- 
inator of (Umi). 

The derivative with respect to nk yields 



1 , / "k 
■ m 



(i Vl + "k 



Ex. 



(fk - MMF)A^k - TMF-Pk 
i?k 



(8.17) 

which, upon substitution of (18. 7|) and (|8.14|) . yields the 
single particle excitation spectrum 



^k == y (ck - Mmf)2 - 7mf' (8-18) 

with an energy gap 

-E'gap = £^k= 



One identifies in addition, 



-^k = Ck — MMF, ^k = 7MF- 



.20) 



By substituting these results into (|8.10p and (|8.1ip . one 
finally obtains the self-consistency conditions. 



«MF 



dk 

{2'kY 



$ 



MF 



-7mf 



£k — AtMF 
Sk 

dk 1 

J2^E^ 



nk+2 



nk 



(8.21) 
(8.22) 



In the normal phase, Vo — 0, $mf = 0, Eq. (|8.22p 
is automatically satisfied, and (|8.2ip determines nMF- 
In the MSF phase, Vo = but <i>MF t^ 0, giving 
7mf = (3(7 -I- wnMF)^MF and leading to a self-consistent 
("gap" -like) equation ([5?^ . that together with (jOTjl . 
determines 'I'mf and nMF- Since the integral in (|8.22p is 
positive, clearly, there is a nontrivial solution, $mf 7^ 
only if 7mf < 0, i.e., for sufficiently attractive atomic 
interactions {g sufficiently negative). 

Finally, the derivative with respect to -00 yields 



1 



= -fi + gi/jo + -jW%l)Q, 



.23) 



that, together with (I8.2ip and (I8.22p . self-consistently 
determines tpQ inside the ASF phase. The limit '00 ~* 
then yields the ASF-MSF phase boundary at the crit- 
ical value jlc — 0. At this point (and only at this 
point) /iMF = 7mf, and the energy gap (|8.19p there- 
fore vanishes at the ASF-MSF transition. This is consis- 
tent with the results of the two-channel (atom-molecule) 
model [see (|63T|) and ([QO]) ]. where in both the ASF 
and MSF phases the atomic branch _Eki (that corre- 
sponds to spectrum E\^, above) remains gapped except 
at the MSF-ASF transition point. Because, in the one- 
channel model, molecular excitations do not explicitly 
appear in the Hamiltonian [though molecular superfiuid 
order clearly does appear, via anomalous averages (|8.1ip ]. 
neither does the gapless spectrum of the corresponding 
(molecular) superfiuid mode. However, as will be seen in 
the discussion in Sec. IVIII Bl the presence of these gap- 
less molecular modes will appear in the calculation of the 
superfiuid density. 

With the above substitutions, the free energy density 
([5TT^ simplifies to 



/. = 



rfk 



{2ttY 

fik — MMF 



fcerinfl-e-^^'' 



7mf „ 
Ek 



■nk 



Ck — Mmf 
^^k 



•?^o(^MF, 'I'MF, V'o) + A^MF^MF- 



(8.24) 



The similarity of (|8.24p to the two-species form (|6.66p is 
evident. 



=0 - Y^MF -7mf- (8-19) [df^/duMF)^ 



The extremum conditions (9/t,/9$MF)nt 



0, 



0, in which the derivatives include 
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all dependence in -Ek, MmFiTmf, yield precisely the con- 
straints (|8.2ip . (|8.22p . As a consequence, one also has 

the relation niviF = —{dfv/dfiMF)nMF;'S'MFJ i^^ which the 
derivative includes only the dependence from the combi- 
nation Ek — Mmf (appearing especially in Ei^). 

Equations ([OT|) . ((O^ and ((05)) are the fundamen- 
tal results of this section, providing a set of closed re- 
lations to be solved for timf, *&mf and ipo as functions 
of /i, T. The stabilizing three-body repulsion w plays no 
essential role here: the equations remain perfectly well 
defined for w = 0. This is because the form (|8.3p for 
the variational Hamiltonian already precludes the type 
of real space system collapse against which w is intended 
to stabilize. It should be kept in mind, however, that it 
is only in the presence of w that the type of superfluid- 
order considered here would actually occur. Thus, w > 
motivates the form of the variational ground state, but 
once this form is adopted, w effectively disappears from 
the calculation. As observed in experiments, where w is 
generally quite small and system collapse does eventu- 
ally occur, such states are expected to be dynamically 
metastable in the dilute limit even when they do not de- 
scribe true equilibrium. 



the odd part is invariant under the transformation. 



Y,Mkdi~ak = Y,Mkil% 



k 



.27) 



the Bogoliubov transformation is independent of kg . Af- 
ter the diagonalization, Hmf takes on the following form: 



Hmf — 



k 



i^,t„-,. 



1 



iEk + Mk)%% + -iNi,-Ei,) 



+ [No - Po]iJo ^ QoA]V. (8.28) 

The atom number and molecular order parameter densi- 
ties follow in a form identical to (|8.10p and (|8.1ip . but 
with the Bose occupation factor given by 

1 



nw 



al3{E^+M^) _ I ' 



(8.29) 



that includes the odd function Afk- The variational free 
energy (|8.2p follows in the form (|8.12p . but (i) with a 
single additional term 



AF„ 



dk 



(2^) 



^vo • ?iknk 



(8.30) 



B. Superfluid density 

One can infer the existence of gapless molecular exci- 
tations in the single species model from the existence of 
a nonzero superfluid (number) density, ris = {m/h^)Ts 
defined (as before) in terms of the change in the free en- 
ergy, AF^, (|6.98p associated with imposition of twisted 
boundary conditions on ip{r). 

As in (|6.97p - (|6.99p . one expresses H in terms of the 
periodic field operator ip{r) = e~'^^°''^ip{r), with the result 



H[^P] = H[ij] +eoN + ^^o■P 



(8.25) 



where Eq = fi?kf^/2m, vq = ftko/ni, and iV, P are given 



by (|6.100p with ■0 replacing ^^ and dropping a summa- 
tions. In the presence of twisted boundary conditions one 
generalizes the variational Hamiltonian to the form 



H^ 



MF 



E 

k 



(A^k 



Mk)akak 



:Pk 



4«-k 



H.c. 



V[~ao 



4 



(8.26) 



in which Ok is the Fourier transform of tp. All variational 
coefficients depend on the twist wavevector ko, but A'k, 
Pk are even functions of k as before, while A/k is an odd 
function of k. In fact, as will seen shortly, Mk = fi.Vo • k, 
but this is not assumed at the outset. The shifted Bogoli- 
ubov transformation (|8.4p - (|8.7p is performed in an iden- 
tical fashion, with ak replacing Ok, and itk, I'k depending 
only on the even functions A'^k, Pk- As in (|6.104p . because 



arising from the vq • P term in (|8.25p . and (ii) with /j, 
replaced by fi — so everywhere, arising from the eqN term 
in (|05|) . 

One performs the variational minimization as before, 
treating nk, i?k — Pk/A^k, and -00 as independent vari- 
ational parameters. Minimization over Pk yields (I8.14p . 
but again with /x replaced by /i — Eq everywhere in (I8.16p . 
Minimization over rik yields 



Pk + Ai"k = Y (ek - A^mf)^ - 7mf + '^o ■ ^-k, 
and one immediately recovers (j8.18p for Pk, and 
Mk = ^vo • k, 



.31) 



(8.32) 



as promised. Minimization over ipo recovers (|8.23p , again 
with /i replaced hy fj, — Eq everywhere. 

The superfluid density, Ug, Eq. (|6.96p is proportional 
to the second derivative of /„ with respect to fco, Eq. 
(I6.98P . Since there is an implicit AiQ-dependence through 
all terms in the free energy, the computation of Tg ap- 
pears overwhelming at first sight. However, two observa- 
tions simplify it enormously: (1) All single /co"derivatives 
of even functions of fco are odd functions of fco , and hence 
vanish at kg — 0, and (2) the variational property implies 
that all single derivatives of the free energy with respect 
to Pk, n^ of "00 vanish identically. Thus, for example, 
(1) implies that ((9nMF/9fco)fea=o = {d^MF/dko)ko=o = 
{d4'o/dko)ko=o = 0, and hence that cross terms 
such as {d"^ fv / dhMFd^MF)idnMF / dko){d^MF / dko) van- 
ish in the limit fcg -^ 0. Similarly, (2) implies 
that {dfJdni,)idW/dk^) = (9/./9Pk)(9'Pk/aA:2) = 
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The result is that there are only two contributions to 
the superfluid density. The first comes from the fi ~ Eq 
combination, and yields a term 



dfv d^eg 



— nMF, 
m 



(8.33) 



where timf = —df^/dji = V'o + "-mf is the number den- 
sity [the variational conditions again imply that the only 
contributions to the //-derivative come fr'om the explicit 
dependence in Tq — see (|8.13|) ]. The other contribution, 
interpreted as the normal fluid density, comes from the 
term (|8.30p and together these give: 



K^ 



T, 



lim 

2 
d 



nMF - n, 
dk 






dk 



ek- 



dn-L. 



{2tj:Y dEi, 
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in which kg is the unit vector along kg, and kg has been 
set to zero inside rik in the last line. The resemblance to 
(|6.106p is clear. As promised, the superfluid density is fi- 
nite even though £'k is gapped in both the ASF and MSF 
phases, indirectly indicating the presence of the gapless 
molecular excitationsi^ 



C. Solutions to the variational equations 

In what follows /imf will be treated as the indepen- 
dent control parameter, and 7mf viewed as fixed. If one 
wishes, one may use solutions to (|8.16|) and (|8.17p ob- 
tained in this way, together with (|8.12p . to solve in the 
end for the behavior as a function of /i at fixed g, w. For 
simplicity only T = where rik = will be considered 
here. 



1. Onset of molecular superfluidity: vacuura-MSF 
transition 

In the MSF phase the order parameter constraint 
([OS)) reduces to 



.9m F 



dk 1 

(2^:^l'"^"^2 



.35) 



in which ^mf 
ond line of (H 
at T 



iwriMY, and /^mf is given by the sec- 

The onset of molecular superfluidity 

takes place at the value /x = /xq at which parti- 



= 9- 
13. 



cles first begin to enter the system, i.e., at zero density. 
Letting timf, 'I'mf —> 0, (|8.35p yields 



2 
9 



dk 



1 



where the molecular size is much larger than the diam- 
eter do of the attractive potential. Comparing (|8.36p to 
()C3P (in the case where {"(k) = 1 — note that an effective 
cutoff /cA ~ 7r/do is required in this case to regularize the 
integral for d > 2) one sees that /xq = —E/2. Therefore 
the onset of molecular superfluidity indeed occurs pre- 
cisely when the chemical potential rises above the molec- 
ular binding energy per particle. This confirms that TYmf 
correctly captures this limit. 



2. MSF-ASF phase boundary and the closing of the single 
particle gap 

More generally, the MSF phase single particle spec- 
trum (|8.18p has a gap 



E, 



gap 



^k=0 



A^mf 5mf^mf' 
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representing the minimum energy required to create a 
single atom excitation. ^^ 

As /x increases from /xg this gap shrinks and vanishes 
at a critical value /Xc such that /x^p = gyiF^M-p- This 
point identifies the MSF-to-ASF transition, and at T = 
(|8.35p reduces to 



d'^k 



1 



(8.38) 



5MF J i^TvY y^ek(ek - 2/i^p) 
Subtracting ((05)) (at T = 0) from (jOS)) one obtains 



= 



d'^k 
(2^ 



v/ek(ek - 2/x 



MF^ 



1 



(8.39) 



in which the integral is now fully convergent for 2 < d < 
4, and the short-scale uv-cutoff may be dropped. The lat- 
ter is now effectively subsumed into a nonuniversal value 
of the critical chemical potential, /i^F' while critical be- 
havior, which depends only on deviations of /xmf from 
this value, remains universal. 

The density integral (|8.2ip is already convergent for 
d < 4, and one obtains the critical value 
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-2m ^l 



MF 



h^ 



rf/2 



(8.40) 



which exhibits a simple power-law relation between n£jp 
and /x£jp. The substitution u^ = —h?k'^/2mA^J-Mp ^^^ 
been used to define a dimensionless constant 



/d = 
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{2TrY 
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VuHi 
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/^ N^ (8-36) 

{2ttY ek - Mo ^ ' 

In App. [C] an equation is obtained for the ground state 

energy E oi a single molecule in the weak-binding limit tion va 



The quantity ri\i — ^j2m[jiyYpjW may be interpreted 
as the background molecular diameter. Thus (|8.40p 
shows that the condition for closing of the atomic gap 
corresponds to a criterion of the mean atomic separa- 



l/d 

*MF 



reaching tm- The MSF-ASF transition 
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therefore takes place in a regime in which pairs begin 
to strongly overlap, which is the condition under which 
atoms can begin to hop from one molecule to another as 
they themselves become delocalized in a sea of extended 
pairs (see Fig. [T2)) . 



from criticality, all of which will turn out to be positive 
in the MSF phase. Thus, in fact, although fj, increases 
as the density increases, /Umf decreases due to the extra 
terms in (|8.15p . In terms of these one obtains from (|8.2ip 
and (|835)) : 



3. Critical behavior near the ASF phase boundary 

The approach to ASF-MSF critical point may be an- 
alyzed by considering small deviations 
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U2 + 1 



m2 + 1 _ r 



v/u2(m2+2) ^(u2 + 1 - t)2 - 1 + / 
I 



(8.43) 



The detailed analysis of (|8.43p is relegated to Appendix 
IeI The results may be summarized as follows. 

The critical behavior is found to change dramatically 
at dimension d = 3, which is to be expected because 
the upper-critical dimension for a quantum Ising model, 
controlling the zero-temperature MSF--ASF transition is 
df,,. = 3. For d > 3 the critical behavior is Gaussian, 



I{r) = r f/« + /(^VC^-^)/^ + 0(r) 



P{r) 



Pd 



(1). 



-H p|;^V('^-3)/2 + 0(r)l , (8.44) 



where the d-dependent coefficients are given in (jE5[) . 
(IE6[) . and one may expect TYmf to provide an asymp- 
totically exact description of the MSF-ASF transition. 

For d < 3 the critical behavior is nontrivial. The vari- 
ational theory predicts 

/(r) = r [/f + /(^M3-'^)/('^-i) + 0(r); 

Pir) = rWP+pfr^'-'^'''-''+0{T)\, (8.45) 



[see (|E12p . (jE13p ]. but this approximation must break 
down sufficiently close to the transition point, determined 
by a Ginzburg criterion that can be worked out in a stan- 
dard way. For a dilute gas, relevant to the atomic gas 
experiments considered here, the size of the Ginzburg re- 
gion should be be very small, and therefore it is unlikely 



that the resulting asymptotic critical behavior (which, 
as discussed in Sec. IVIIi is likely actually a fluctuation- 
driven first order transition) can be observed. The ap- 
parent exponent singularity at d = 1 is a signature of the 
lower-critical dimension below which the ordered phase 
(ASF), and therefore the transition to it, is destabilized 
by quantum fluctuations. 

In d = 3 there are logarithmic corrections: 
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lnln(To/T) ' 
ln'(ro/r) . 

lnln(ro/T) ' 
ln^(To/T) _ 



,(8.46) 



with amplitudes given in (jE17p and (jE18p . 

In all cases, there is a leading analytic dependence, 
linear in r, followed by a subleading singular contribution 



with "quantum speciflc heat" exponent 

-^, d>3 
(log), d = 3 

d-1' "■ < '^■ 



..47) 



Note that a here differs from that in (|6.82p which, al- 
though computed for a slightly different quantity, never- 
theless reflects the same universal energy singularity. The 
change is not due to a change in the universality class, but 
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FIG. 12: Schematic illustration of the single species MSF-ASF transition, accompanying the BEC-BCS crossover, (a) For 
sufficiently strongly attractive interactions, the extent of the molecular wavefunction (^o(x, x') is much smaller than the inter- 
molecular separation d, and is essentially determined only by two-body physics. Interactions between molecules are weak, and 
only rarely do they overlap sufficiently to exchange atoms. At low temperatures the molecules Bose condense, but coherent 
intermolecular hopping that would lead to atomic Bose condensation is suppressed, (b) The MSF-ASF transition takes place 
when the size of the molecular wavefunction (which is strongly renormalized by many body effects) becomes comparable to d 
and molecules begin to significantly overlap. Coherent atomic hopping occurs over ever larger distances as the overlap increases, 
and the divergence of this distance signifies the MSF-to-ASF transition, (c) Deep in the ASF phase, the size of the molecular 
wavefunction is much larger than d, and essentially loses its physical meaning: coherent hopping of atoms over the entire system 
leads to atomic Bose condensation, and one can no longer identify any given pair of atoms with a particular molecule. Instead, 
the physics is described by a BCS-type many body wavefunction, embodied in the Hamiltonian (|8.3|l . which maintains strong 
nonlocal, pair correlations in the absence of identifiable molecules. 



rather to the extra constraint embodied in the first line 
of (|8.43p . which essentially reflects the existence of only a 
single species. It is well known that such constraints lead 
to a "Fisher renormalization" a — > — a/(l — a) whenever 
a > 0, while leaving it unchanged if a < (the con- 
straint therefore always leads to a negative specific heat 
exponent) i^ If one were to enforce the total density con- 
straint n = 2n2 + fii, the condensate depletion (|6.72p 
would also display the Fisher-renormalized exponent. 

One may use a similar analysis to extend these results 
into the ASF phase, and to positive temperatures. How- 
ever, rather than exploring (difficult to probe) critical 
behavior, the aim here is mainly to demonstrate the ex- 
istence of the same three phases (N, MSF, ASF) illus- 
trated in Fig. [21 the role of the detuning v is played here 
by the s-wave interaction parameter g, or more properly 
the corresponding dimensionless measure of scattering 
length, i.e., the gas parameter n^/^a. Despite the com- 
plete loss of molecular identity in the broad resonance 
(single-channel) model as the ASF phase is approached, 
the topology of the phase diagram and critical behavior 
(accounting appropriately for constraints) is the same as 
for the two-channel model. 



IX. TOPOLOGICAL EXCITATIONS 

In the previous sections a description of a (s-wave) 
resonantly interacting atomic Bose gas was presented. 



formulated in terms the atomic and molecular super- 
fluid (condensate) order parameters ^lo and \E'2o, and 
corresponding fluctuations in the two ordered, ASF and 
MSF states were studied, and characterized the associ- 
ated T — and finite T phase transitions. 

In this section, a complementary description of this 
two-component (atoms and diatomic molecules) Bose gas 
is presented, in terms of topological excitations in the 
MSF and ASF phases. These will be shown below to be 
vortices and domain walls. Descriptions of phases and 
phase transitions in terms of topological excitations has 
a long and successful history, with ordinary vortices in 
superfluids and superconductors, dislocations and discli- 
nations in crystalline solids, and domain walls in Ising fer- 
romagnets being only a few most prominent examples.^— 

The importance of this description is two-fold. Firstly, 
topological defects are true nonlinear excitations of the 
system and thus are essential for a full characterization 
of the response of an ordered state to an external per- 
turbation. For example, a rotated neutral superfluid (or 
a superconductor in the presence of a sufficiently strong 
magnetic field) responds by nucleating quantized vortices 
that carry discrete units of fluid's angular momentum 
(magnetic flux) . Secondly, fluctuation- induced (quantum 
or thermal) topological defects provide a dual charac- 
terization of phases and transitions between them that 
complements their description in terms of a Landau-type 
order-parameter. For example, a superfluid-to-normal 
transition can be understood through a dual model of 
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fluctuation-induced proliferation of vortices, with the su- 
perfluid state acting as a vortex vacuum (or an insulator) 
and the normal state as a vortex condensate j2& 

In addition to simply playing a complementary role, 
such dual vortex "disorder parameter" descriptions are 
also a powerful way to characterize subtly ordered phases 
that do not allow a direct Landau order parameter de- 
scription. The most prominent examples of this are 
2d ordered phases with a continuous symmetry that 
are (usually^) "forbidden" by the Hohenberg-Mermin- 
Wagner theore m^^'^^ to exhibit true long-range order 
and thus cannot be characterized by a condensate or- 
der parameter.^^ Such "ordered" phases (e.g., a 2d su- 
perfluid or a 2d crystalline solid) are in fact disordered, 
only distinguished from the short-ranged (exponentially) 
correlated fully disordered states by a quasi-long-ranged 
(QLR) order with correlation functions falling off as a 
power-law. Descriptions of such QLR-ordered phases and 
their transition to fully disordered states is best done in 
terms of a proliferation of topological defects, e.g., vor- 
tices in 2d superfluidsi^ In higher dimensions, a descrip- 
tion in terms of a proliferation of topological defects can 
also be more effective, as for example found in disorder- 
ing of a 3d type-II superconductors by proliferation of 
vortex loops i^i2^ 

Such dual topological defect descriptions, in addi- 
tion to providing added physical insight, provide impor- 
tant complementary computational tools for the study- 
ing these phenomena. With this motivation, topological 
defects in the ASF and MSF will now be considered. 



A. Atomic superfluid 



Since the fully ordered ASF state has two nonzero or- 
der parameters ^io,^20j there are interesting features 
of the topological excitations generated by the (Fesh- 
bach) coupling between them. The thermodynamics of 
the state can be conveniently and equivalently described 
in terms of the local magnitudes and phases of its two 
(atomic, cr = 1, and molecular, a = 2) coherent-state 
fields 



In terms of these, the real-time coherent-state action cor- 
responding to the Hamiltonian (|2.ip takes the form 



iS* — 5*1 -l- 5*2 -l- Si2 
Si = I dtdv 
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where terms involving \Ina , that are less important than 
the finite compressibility g^ terms, have been dropped. 
From the action S the phase diagrams of Sec. IIVI and 
Bogoliubov modes of Sec. IVII can be straightforwardly 
reproduced in terms of the canonically conjugate densi- 
ties ria and phases 0^- 

Note now that the mean field equations of motion 
for the phases 0^, namely the Euler-Lagrange equations 
0, are given by 



ss 
se. 



dtni + —V • (niV6li) 
ra 

dtn2 + T^V • (naV^a) 
Zm 



2^sin(20i-02), (9.3) 



J 



sin(2e'i - 02), (9.4) 



where the internal Josephson coupling between atomic 
and molecular superfiuids. 



J = aniy^, 



(9.5) 



is proportional to the Feshbach resonance amplitude a. 
These can be combined to derived the total boson number 
n = ni + 2n2 conservation (continuity) equation 



9tn + V • j = 0, 



(9.6) 



where the total number current j = ji + 2J2 naturally 
consists of the atomic and molecular contributions, 



Ji 

J2 = 
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h 

2m 



-niW0i 



712^02 



(9.7) 
(9.8) 



Observe that, due to the atom- molecule Feshbach 
resonant interconversion captured by the "current" 
Jsin(26'i - 02) on the right hand sides of ([Q]) and (|01) . 
as expected, ni and 712 are not independently conserved. 
The microscopic action S in (|9.2p is not completely 
generic. A more general model (that can be obtained 
either based on symmetry or by incorporating quantum 
and thermal fluctuations) includes an atomic-molecular 
current-current interaction of the form 



Ipa 
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<55i2 = -i^i2|V(20i-02)r, 



(9.9) 
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arising from coarse-graining of the action S in the pres- 
ence of the Feshbach resonance cosine nonhnearity. The 
form of this term ensures that total atom conservation 
embodied in the continuity equation (j9.6p remains satis- 
fied. 

Fluctuations lead to corrections to the mean field equa- 
tions of motion. At the hydrodynamic level, in which 
only the dynamics of slow, large scale distortions of the 
fields are considered, these corrections may be embodied 
simply in renormalization of the terms appearing in S. 
Most significantly, in this limit phase fluctuations domi- 
nate, and fluctuations in Ua- may be subsumed into renor- 
malized stiffness coefficients. Thus, the squared phase 
gradient terms undergo the replacement 



2rn, 



-n,|V0,| 
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:KJve„ 



K„ 



(9.10) 



which replaces the fluctuating number density n^ by the 
atomic and molecular superfluid (number) densities, risa-, 
to be distinguished from the corresponding, quite dis- 
tinct (see Sec. IVI[) condensate fractions noo-. As in a 
single-component superfluid at T = 0, Galilean invari- 
ance enforces the condition 



2n, 



n, for T = 0, 



(9.11) 



namely that the total superfluid density equals the total 
density. 

With this preface, the focus will now be on the finite 
temperature classical limit, ignoring quantum dynamics 
that are left to future investigation. The model to be 
studied is defined by the "hydrodynamic" energy density 



£ = -K,(ye,f 



^K.iye^f 



-ifi2|V(20i 



Jcos(26ii-6'2), 



(9.12) 



and is valid in any region where the phases 9a are well 
defined, but must be supplemented by core energies in 
regions where an order parameter magnitude vanishes. 
The four coefficients are all renormalized hydrodynamic 
parameters that depend on the chemical potentials, and 
other microscopic parameters, and have their own non- 
trivial critical behavior,^?- The detailed knowledge of 
their exact values (some of which have been computed 
in the dilute limit in earlier sections of this paper) is not 
required to understand general features of topological ex- 
citations. 



1. Vortices in the ASF 

In the absence of the atom-molecule couplings, K12 
and J, the superfluid admits independent atomic and 
molecular vortices — pictured in Fig. 1131 Focusing for 
simplicity on d = 2, these are point defects in the atomic 
and molecular superfluid order parameters, around which 
their respective phases wind by an integer-multiple of 2tt 
as the point is encircled. As usual, this quantization con- 
dition arises from the single- valuedness of the superfluid 
order parameter away from the vortex core (located at 
position rocr): 



the corresponding vortices 



V6a-dr^2TTpa, 
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with "charge" Pa- 

Imposing this topological constraint and minimizing 
the energy £ at K12 = ./ = 0, one obtains indepen- 
dent atomic and molecular superfluid velocities around 
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with integer charges Pa. Equivalently, the superfluid 
phases 9a- are given simply by integer multiples of the 
azimuthal coordinate angles ipa (measured with respect 
to an origin chosen at the vortex core positions roo-), with 
9a = Pa'fia- As usual in 2d, vortex energies grow loga- 
rithmically with system size L, 



E(^"'>^plKa7rHL/^a), 



(9.15) 



where S,a is the vortex core size set by the corresponding 
coherence lengths. 

In the presence of the inter-superfluid couplings K12, 
J no general solution is available. However, considerable 
insight can be obtained by analyzing limiting regimes. 
It is clear from the energy density £, Eq. (|9.12|) . that 
to avoid extensive (scaling with system size) energy cost 
proportional to J (and/or K12), the two phases are on 
average locked together according to 



^29i. 



(9.16) 



Hence the energy is minimized when positions of the 
atomic and molecular vortices coincide, and their topo- 
logical charges (winding numbers) are related by p2 = 
2pi. Thus an elementary pi = 1 vortex in the atomic 
order parameter will be accompanied by a spatially co- 
incident molecular vortex of topological charge that is 
double its elementary value, p2 — 2. 

In contrast, an elementary P2 = 1 molecular vortex is 
energetically significantly more costly due to incompat- 
ibility of the atomic order parameter single-valuedness 
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(a) 



(b) 




FIG. 13: (a) A unit vortex in the order parameter field ^io(r) 
represented as a 27r-rotation in the vector field Mo(r). (b) In 
the ^10 complex plane, a unit vortex in ^20 (r) is represented 
by a TT-rotation in the double headed vector field Qo{r). One 
identifies the Feshbach resonance coupling — aRe[^2o'iI^io] BiS 
the energetic tendency for spatial alignment of Mo(r) and 
Qo(r)- Details of the mapping between these two representa- 
tions are provided in Appendix El 



constraint and the Feshbach resonance constraint (j9.16p . 
It is clear that energetically there are two competing, 
least costly, ways to accommodate this frustration. One, 
illustrated in Fig. [Hf a). is with a spatially coincident 
half-integer (pi — 1/2) atomic n vortex, that requires 
that ^10 vanish (atomic component of the gas is normal) 
along a ray emanating from the location of the vortex 
core, and across which 61 exhibits a tt jump discontinu- 
ity. The cost of such a defect clearly scales linearly with 
the length, L of the defect ray (more generally, as L'^~'^ 
in d dimensions) and is dominated by the loss of the con- 
densation energy along the linear defect. 

Another competing possibility, illustrated in Fig. 
I14r b). is an atomic pi = 1 vortex localized on the elemen- 
tary molecular vortex, but (in contrast to the noninter- 
acting case where the vortex is isotropic with Oi{(p) = ip) 
the atomic phase winding is highly anisotropic, with 



lM 



02(^)/2, 



(9.17) 



outside a narrow domain wall strip. The atomic phase 
makes up the remaining tt deficit angle (required by 
single- valuedness of the atomic order parameter), by 
rapidly winding across the domain wall of width set by 
J and a combination of phase stiffnesses K^r- 

It is clear that the first scenario is the limiting case of 
the second configuration, with large J and small K^, such 
that the wall width is microscopic (formally smaller than 
^a- ) and the corresponding energy comparable to conden- 
sation energy, thus driving the discontinuity ray normal. 
The resulting energy in both cases clearly grows linearly 
(as L'^^^ in d dimensions) with the length of the domain 
wall ray, and the energy scale is set by the minimum of 
the condensation energy or line tension, with the latter 
given by the geometrical mean of J and a combination 
of the Kcr (see below) . 
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FIG. 14: A schematic iUustration of a 27r (elementary unit) 
molecular vortex that induces a tt (half-unit) atomic vortex, 
that in turn induces a domain- wall ray. In (a) the wall width, 
^ is smaller than the coherence length and the energy cost per 
unit of wall length exceeds that of the condensation energy, 
thus leading to a "normal" ('I'lo = 0) domain wall. In (b) the 
superfluid stiffness is large and Feshbach resonance is narrow 
(small a) leading to a wide domain wall (width ^ exceeding 
the coherence length), with an interface that is in the ASF 
state and "tio only slightly suppressed below its bulk value. 
The wall structure is illustrated in more detail in Fig. 1151 



2. Domain walls in the ASF 

It has been argued that the domain wall with energy 
linear in its length L (more generally, growing as its sur- 
face area L"*^^ in d dimensions) is another type of a topo- 
logical excitation in the ASF. Although in the previous 
subsection it emerged as a necessary string component of 
a p2 = 1 molecular vortex, the existence of a domain wall 
excitation can be understood on more general grounds. 
Quite similar to domain walls in an Ising ferromagnet, 
here too it is a defect that separates ordered domains 
associated with two physically distinct configurations of 
the Ising order parameter in the ASF that spontaneously 
breaks the Z2 symmetry of the MSF state. In terms of 
phases 9a- the two domains correspond to two solutions 

6*^"^ = 6I2/2 + TiTT, n = 0, 1, (9.18) 

of the constraint in (|9.16p , that are associated with two 
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values of the atomic order parameter ^]^ 
±e**"/^, pictorially illustrated in Fig.fTSl 

A detailed solution for a domain wall can be straight- 
forwardly worked out. To this end, the key (internal) 
Josephson nonlinearity associated with the Feshbach res- 
onance is isolated by a convenient change of phase vari- 
ables to new phase fields 6 and 0: 



1(20, 



(9.19) 
(9.20) 



corresponding to the in-phase and out-of-phase fluctua- 
tions of 61 = {9 + 4>)/2 and 02 = 9 — (j) phases. 
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(j)=7I 



# 



(t)=0 



FIG. 15: Details of the domain wall structure, separating 
'I'j^ = ±1 Ising domains of the ASF corresponding to Sj = 
S2/2 and 6*^ = ©2/2 + vr, respectively (illustrated here for 
Kerjy/Ke = 1). Across the wall the atomic condensate phase, 
01, winds by vr relative to the molecular condensate phase 
(double-headed arrow angle), 62- 



In terms of 9 and 0, the energy density is given by 
£ = \Kg{Vef + \K^{Vd^f - Kg^Ve ■ V(b ~ Jcos(20), 



with 






(9.21) 

(9.22) 
(9.23) 
(9.24) 



The corresponding saddle-point equations 5E/59 = 
58/54) = are given by 

-KeV^e + Kg^V'^dp = 0, (9.25) 
-K^V'^(j) + Kg^V'^e + 2Jsm2(j) = 0. (9.26) 



Eliminating 6 via (|9.25p reduces (|9.26p to the well-studied 
sine-Gordon equation for 0, 



where 



i\:v2(/) + 2Jsin2<^ = 0, 



K^K^- KljKe. 



(9.27) 



(9.28) 



For a straight domain wall oriented along x, defined 
by the boundary conditions 4>dw{y — > —00) = and 
(j)dw{y -^ +00) = TT, the solution is given by 

0dw(2/) = 2 arctan (e^?^/^) , (9.29) 

illustrated in Fig. [151 The domain wall width is given by 



e = v/W- 



The associated 9dw{y) and corresponding 0'^^{y) solu- 
tions can now be easily obtained from (|9.25p . (|9.19p . and 
(lOOll : 
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1) </>d^(y), (9.32) 
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Using the above expressions inside the energy den- 
sity £, Eq. (|9.2ip . and taking advantage of the Euler- 
Lagrange equations (|9.25p . (|9.26p one finds that the do- 
main wall energy is given by 



Edw = / dxdy K{y(j)dwf 



(9.34) 

dxdy 2J[l - cos(20d^,)], (9-35) 

= o-dw^x, (9.36) 

with domain wall line tension (energy per unit of length) , 

CTdw = 4V7^. (9.37) 

3. Point-to- "dumbbell" atomic vortex transition 



As is clear from the discussion in Sec. IIX All in the 
ASF a 27r (pi = 1) atomic elementary vortex is driven 
by the Feshbach resonance (internal Josephson) coupling 
J to be accompanied by a 47r {p2 = 2) molecular double 
vortex. In the limit where J ^ K, corresponding to a 
broad Feshbach resonance, and deep in the ASF state, the 
two supcrfiuids arc strongly coupled, and will behave as a 
single-component conventional atomic superfluid. Thus 
the 27r atomic and 47r molecular vortices must spatially 
coincide. This leads to an isotropic topological defect 
with energy (measured relative to the background energy 
— J of the uniform state) given in 2d by 

^(point) ^ ^(point) ^ ^j^^ ln(L/a) + 47ri^2 ln(L/6), 

(9.38) 
where 



^^point) ^ ^(2-k) ^ ^ 



(477) 



(9.39) 



(9.30) 



consists of the atomic (27r) and molecular (47r) vortex 
core energies. 

Such a concentric, isotropic vortex configuration mini- 
mizes the Feshbach resonance energy. However, because 
it involves a 47r {p2 = 2) molecular vortex that is double 
the elementary charge, it raises the large-scale part of 
the kinetic energy by 27ri4r2 ln(L/^2) over the energy of a 
topologically equivalent vortex configuration consisting 
of two elementary 27r (p2 = 1) molecular vortices — see 
(|9.15p . As will be shown below, this can drive the split- 
ting of the 47r double molecular vortex into its elementary 
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27r constituents. This is driven by the fact that, in the ab- 
sence of the atomic component, two elementary vortices 

repel via a potential V.} ^ ^ (R) — —2-!TK2ln{R/^2), 
where R is the separation. On the other hand, as shown 
in Sec. lIX ATI in the ASF phase an elementary 2tt molec- 
ular vortex is driven by the internal Josephson coupling 
J to be accompanied by a (pi = 1/2 fractional) atomic 
IT vortex and a domain-wall string defect. Thus, in the 
ASF state the logarithmic repulsion of two 2tt molecular 
vortices is arrested by the confining (Josephson coupling) 
domain- wall energy that, according to (|9.36p . grows lin- 
early with separation R. Hence, the energy E^^" of a 
point vortex (consisting of coincident 27r atomic and An 
molecular vortices) must be compared to a topologically 
equivalent "dumbbell" configuration split by a separation 
R into two units, each consisting of coinciding tt atomic 
and 27r molecular vortices — see Fig. [5l^^ 

The energy of the dumbbell configuration is estimated 
as 

^(dmbi) _ Ei'^"'^^'^ +nKi\nL/R + AnK2\nL/R 

- ttKi ln(i?/a) - 47ri^2 ln(i?/6) + ^dwi? 

(9.40) 
« £;(dmbi) ^ ^^^ ln(i/a) + 4ttK2 ln(L/6) 

- 2TTKi \n{R/^i) ~ 8ttK2 ln(i?/6) + a^^R, 

(9.41) 



where the core energy is 

^(dmbl) ^ 2 \e[1^ + eH^^ 



(9.42) 



The first two logarithmic terms in (|9.40p estimate the 
energy associated with the total 2tt atomic and total An 
molecular topological charges outside the dumbbell of 
size R. The second two logarithmic terms in (|9.40p give 
the energy of two tt atomic vortices and two 2n molecular 
vortices, including their repulsive interaction at separa- 
tion R. Finally, the last term accounts for the energy 
of the domain wall in the atomic superfluid. Accuracy of 
the estimate requires that the system size be much larger 
than the dumbbell size, which in turn must be larger than 
the microscopic scale: L :^ R^ ^i,S,2- 

To estimate the optimum size Rq of the dumbbell 
vortex configuration, one minimizes Ev ™ (i?) over R, 
yielding 



Rq 



TT Kl+ AK2 




(9.43) 
(9.44) 



where (|9.37p has been used for ddw, and the estimate in 
the last line is made by ignoring depletion effects, valid 
in the dilute, weakly interacting limit. 

It is clear that there is a range of parameters K„ and J 



dumbbell vortex is lower than that, Eq. (|9.38p . of the 
point vortex configuration. In this case, a 2d system will 
undergo a transition to a state in which the 2tt atomic 
vortices are in a dumbbell tt—tt configuration (locked to a 
2tt — 2tt molecular vortex pair). Because there is no sym- 
metry change associated with a transition into a state 
where vortex dumbbells are randomly oriented, one ex- 
pects this transition to generically be first order J^ Since 
the energy balance between the two competing states in- 
volves core energies, a more detailed microscopic anal- 
ysis (that is not pursued here) is necessary to pinpoint 
the location of the transition. However, it is clear from 
the structure of M "™ ^ (^0) and i?o, that such transi- 
tion takes place for a sufhciently large domain- wall width 
^J K I J 3> ^CT, where the line tension is small and same- 
sign vortex repulsion is large. One thus expects such a 
transition in narrow (small a) Feshbach resonance sys- 
tems. 



B. ASF-MSF as a confinement-deconfineinent 
transition 

Observe that the atomic condensate density nio van- 
ishes as the detuning v decreases towards a critical value 
Vc- Thus the 2tt dumbbell length 



Rr, 



2tt 



1/2 1/4 
ho ^20 



for V 



(9.45) 



such that the energy Ey 



(dumbbell) 



(Rq), Eq. dnH]), of the 



diverges along with the associated domain-wall width, 
£,. Therefore, the ASF-MSF transition in d = 2 can 
be complementarily described as a 2tt molecular vortex 
deconfinement transition. While 2tt molecular vortices 
are confined by a /mear potential inside the ASF state, in 
the MSF state this confining potential (in 2d) is replaced 
by a much weaker logarithmic potential, that binds each 
2tt molecular vortex to its oppositely charged partner i^Si 



C. Molecular superfluid 

Since molecular 2tt vortices only appear as neutral 
dipoles in the MSF phase, the state is characterized 
by long-range order in the molecular order parameter, 
^20 ~ (e'^^) (however, as usual in 2d, at finite T, 
^20 itself vanishes, while the molecular helicity modu- 
lus, or superfluid density, 71^2 remains finite). On the 
other hand, a deconfined domain wall (across which the 
atomic phase 61 jumps by tt) leads to a vanishing of the 
atomic order parameter, *i -- (e'^i) = 0}°^ The MSF 
state exhibits ordinary molecular 27t point vortices, along 
with the atomic and molecular Bogoliubov quasiparti- 
cles discussed in Sec. |VIl It is easy to see that, analo- 
gous to the conventional BCS superconductor, here too 
an atomic Bogoliubov quasiparticle (that is gapped) ac- 
quires a phase of tt upon encircling a molecular 2tt vortex. 
Thus, these two excitations interact strongly with each 
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other, with statistical-hke interactions that can be cap- 
tured by a Chern- Simons field theory J^ 



X. SUMMARY AND CONCLUSIONS 

In this paper the thermodynamics of a resonant atomic 
Bose gas has been studied. Working within a two-channel 
model, formulated in terms of bosonic atoms and their 
diatomic molecules, the complete phase diagram of the 
system has been worked out as a function of temperature 
and detuning, the properties of the phases, and the na- 
ture of quantum and classical phase transitions between 
them, studied. This analysis was supplemented by a vari- 
ational calculation on a one-channel model, whose salient 
results appear in Sec. lIBI 

A most notable feature is the appearance of two dis- 
tinct superfluid phases, ASF and MSF, separated by an 
Ising type transition. These are distinguished by the 
respective presence and absence of atomic off-diagonal 
long-range order, atomic (gapped and gapless) Bogoli- 
ubov spectra, and the nature of topological excitations. 
In addition to a distinction based on the atomic mo- 
mentum occupation number, these phases can be dis- 
tinguished through the domain wall excitations in the 
ASF (which separate regions in which the atomic phase 
aligns with the molecular phase in the two possible dif- 
ferent ways), characteristic of the broken discrete Ising 
symmetry. 
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APPENDIX A: DETAILS OF CONNECTION TO 
POLAR-NEMATIC ORDERING 



In this Appendix a connection, stated in Sec lIIIi be- 
tween the two complex scalar order parameters, charac- 
terizing phases of a resonant atomic Bose gas, to those 
of a vector-tensor model of polar, nematic liquid crystals 
is elaborated. 

In a thermodynamic description, the order parameters 
are derived from the full free energy density /am via 
derivatives with respect to their conjugate fields: 



*ao(r) = -2 



fdl 



AM 



V 9h* 



= (Mr))- 



(Al) 



This prescription is unambiguous for the two-species 
model (|2.ip . but in the one-species model (|2.15p or (|8.ip . 
^20 remains to be defined. One expects a molecular com- 
posite operator of the form 

i^l (r) = / dr0o {r)^l (r + r/2) V'J (r - r/2) . ( A2) 



to play the role of V'2 , where (f>Q (r) is the molecular wave- 
function. Therefore, to investigate molecular superfluid 
ordering, one is motivated to look at the anomalous cor- 
relation function^^^ 



<i>2o(r,r')-(^i(r)^i(r')). 



(A3) 



Since the fundamental object is a two-point function, 
one generally lacks a unique definition of the one-point 
quantity \E'2o(r). However, if the molecular size is much 
smaller than their separation, in the spirit of the coarse 
graining picture of the molecular operator -02 , one could 
define 



*2o(r)= /"dr>o(r-r')$2o(r,r') 



(A4) 



in which the two-particle molecular wavefunction (f>o is 
used to weigh the local volume average. 

A complex-scalar atomic superfluid order parameter is 
clearly isomorphic to a 2d vector order parameter Mq = 
%/2(Re^io, Im^t^io) whose components are the averages 

Afo,i(r) = (Q(r)), MoAr) = (^(r)) (A5) 

of the corresponding conjugate Hermitian operators, 
1 



Q(r) = — Mr)+i^lir) 



P(r) 



1 

1^2 



-01 (r) -■i/iUr) 



(A6) 



obeying the commutation relation [Q{r), P{r')] — iS{r — 
r'). Phase symmetric combinations of ^I^io (i.e., prod- 
ucts of its complex conjugate with itself) correspond to 
rotation invariant combinations (i.e., dot products) of 
Mq. Similarly, one defines the 2d conjugate field vec- 
tor H = -y={Rehi,linhi) . 

Since the molecular order parameter is fundamentally 
a two-point correlation function, one is motivated to ex- 
amine the tensor 



Q(r,r') = 



(g(r)Q(r')) (O(r)P(r')) 
(P(r)O(r')) (P(r)P(r')) 



(A7) 



1. 



= Qs{r,r') + Qa{r,r') + -Tr[Q(r, r')]l 
in which the symmetric traceless part is given by 
Q,(r,r') = i{Q(r,r') + Q^(r,r')-Tr[Q(r,r')]} 



Re$2o(r,r') Im$2o(r,r') 
Im$2o(r,r') -Re$2o(r,r') 



(AS) 
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and is therefore seen to be precisely equivalent to the 
anomalous correlation function $2o(r, r')- The antisym- 
metric part 



Sa(r,r') = i[Q(r,r')-Q^(r,r')] 

-ImG(r,r') 

ImG(r,r') 



and the trace 



Tr[Q(r,r')]=ReG(r,r') 



(A9) 



(AlO) 



are given by the real and imaginary parts of the usual 
(non-anomalous) two point correlation function 



G(r,r') = (V^i(r)^I(r')), 



(All) 



and are therefore phase invariant scalars. The informa- 
tion about the molecular order parameter therefore lies 
entirely in Q^. One similarly defines the conjugate field 
tensor 



'Hq — - 



1 / Rc/i2 Imft,2 



2 \ Im/i2 — Rc/i2 



(A12) 



By substituting Qs(r,r') into (|A4[) one obtains a def- 
inition of the local tensor order parameter Qo{r). The 
following identities now follow: 







Im,^ 



Ho -Mo 

^Tr(Qg) 

Tr(7^oQo) 

iM;^QoMo 



|*ior 
Re[/ii*io] 
- det [Qo] = 
Re[/i;*2o] 
Re [vl'^o^io] 



*2 



(A13) 



This demonstrates that the phase invariant combinations 
in the Ginzburg-Landau free energy correspond to 2d ro- 
tation invariant combinations in the vector-tensor repre- 
sentation. 

A symmetric traceless tensor order parameter is famil- 
iar from a theory of nematic liquid crystals, encoding the 
headless-arrow nature of the nematic state of anisotropic 
molecules i§"i^^ Using the "dictionary", Eqs. (|A13[) . in the 
vector-tensor representation, Mq, Qq, the mean field 
Hamiltonian (|4.ip takes the form 

Z o 

- ^/i2Tr (Q2) + \g^ [Tr (Q^)]^ _ Tr(7ioQo) 

+ iffisM^Tr (Qg) - |mJQoMo, (A14) 

representing the theory of an interacting vector and ne- 
matic order parameters. In the current scalar superfluid 
context both have two components, but (|A14[) is clearly 



not limited to this case. The conjugate field Hq is the 
analogue of a magnetic field, while Tip is the analogue of 
a nematic liquid crystal polarization field. 

The "double-headedness" of the nematic order param- 
eter is exhibited via the eigenvector decomposition 



Qo(r) = 9o[n(r)n(r) - m(r)m(r)], 



(A15) 



in which go ^ is the order parameter and m, fi arc or- 
thonormal unit eigenvectors characterizing the nematic 
order of the MSF. Clearly Qq is invariant under sign re- 
versal of the unit vectors. Thus, although q^ and ii com- 
pletely define Qqj both n and — n characterize the same 
state. 

To make contact with the complex molecular or- 
der parameter, let n = [cos(6'„),sin(0„)], ih = 
[— sin(6'„), cos(6'„)]. From (|A8[) one obtains 



*2 



90 [("-1 
9oe2^«" 



niT^ + i(h\n2 ~ ihifri'i) 



(A16) 



Thus the eigenvalue q^ — |^2o| is the MSF order param- 
eter magnitude. Although 0„ and 0„ -f tt are equivalent, 
02 = 20„ is uniquely defined. 

Inserting (|A15p into (|A14p . the Feshbach resonance 
phase coupling (a) term takes the form 



-ago 



^i^2 



(Mo • fi) 



-M^ 



-aq^Ml 



= --agoMoCos(02-20i), (A17) 

where the representation Mq = A/o[cos(6'i), sin(6'i)] has 
been used. The coupling is again clearly invariant under 
^n — *■ ^n + TT, and represents an alignment between the 
nematic and polar order parameters, familiar from the 
theory of polar nematic liquid crystals.—"^ 



APPENDIX B: DETAILED ANALYSIS OF MEAN 
FIELD PHASE DIAGRAM 

In this Appendix a more detailed analysis of the mean 
field phase diagram, summarized in Sec. IIVBI especially 
Figs.fTOlandfTl] is presented. To this end, it is convenient 
to introduce the following scaled quantities: 



n 



912 

— 2"Mi, r2 



91 31.92 

-— M2, 7= ^- 
a 5i2 



(Bl) 



7 V5i5i2 - .912,^, ^ 9i9i2nj 

■01 = *10, ■02 = *20, ^ == 2— H 



Substituting these into the mean field Hamiltonian (|4.ip . 
one obtains the dimensionless form 



n = n\i;i\' + l\i'i\^ + rM2\' + li\^2\^ 

+ |Viin^2p-Re[0;Vi?] 



(B2) 
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The scaling has succeeded in reducing the problem to its 
essentials, removing a entirely, and subsuming all inter- 
action strengths into the single parameter 7. The only 
free parameters are the two scaled chemical potentials, 
r^ . One may always choose the phase of 'tpi so that it is 
real and non-negative. It is then clear from the last term 
in (|B2[) that H is minimized by taking ■tjj2 to also be real 
and non-negative, consistent with (|4.3p . With this input 
the scaled extremum equations take the form 

= ^1 (n -^2 + ^1 + {[,1) (B3) 

= -^V^? + V^2(r2+7^2+^D' (B4) 



and lead to the phase diagrams for 7 > 1 and 7 < 1 
shown in Figs. [T7] and [THl respectively (and reproduced 
in Figs. [TOl and fTTj) . These results will now be derived in 
detail. 



Continuous transitions: 



It follows from (|B4| that if V'l = 0, then cither V'2 = 
(normal phase) or -0| = —''2/7 (MSF phase, existing 
only for r2 < 0) . This is the standard result arising from 
7Y = r2'0| -f \"t^2i showing that the N-MSF transition 
must be continuous, and take place at r2 = 0. 

Substituting V'2 = into (|B5|) . one concludes that any 
continuous N-ASF transition must take place along the 
line ri = (with r2 > 0). For ri < 0, it is easy to 
check that 7Y is locally unstable to nonzero V'l- Hence, 
if the transition is first order, representing some global 
instability, it must occur for ri > 0. 

Similarly, substituting f/;! — —r^j^ into (jBSp . one 
finds that any continuous ASF-MSF transition must take 
place along the curve 



n = ric{r2) = \J\r2\h + 'T^h- 



(BIO) 



1. ASF free energy in terms of ?/;2: 

One begins by expressing the energy density in terms 
of -02 alone by using (|B3[) and (|B4[) to eliminate -01 . In 
the ASF phase (the transition to which, from the N and 
MSF phases, is the main focus), where ^i 7^ 0, (jBSp gives 



V'l =02-02 -''I7 



(B5) 



which clearly requires that the right hand side be posi- 
tive. The resulting condition. 



f"! < 02 - 0; 



2) 



(B6) 



will be important in what follows. 

Substituting (JB3| into the first hue of (|B4|) . a cubic 
saddle-point equation is obtained purely in terms of 02 



= 



di>2 



= ri+ 2<202 



301 



2(7. 



^)i'l 



(B7) 



and the corresponding energy density 



T^ASF = ~\rl + ri02 + i20i + i^l + \{l - 1)02, (B8) 



For ri < Tie it is easy to check that 7Y is locally unstable 
to nonzero 0i. Hence, any first order transition must 
occur for ri > ric- 



3. First order N— ASF transition: 

Consider now the possible existence of a first order 
N-ASF transition. For ^2 > all coefficients in (|B8l) . 
aside from the constant term ri, are positive. A positive 
02 root therefore exists only for ri < (vanishing as 
ri —t 0^), and is clearly unique — see the upper panel in 
Fig. [161 The continuous transition at ri = therefore 
takes place for t2 > 0, i.e., r2 > ^■ 

On the other hand, as illustrated in the lower panel 
of Fig. [TBI for ^2 < 0, Hasf has a local minimum at 
some positive value of "02 • Thus, for a range of ri > 
0, there will be two non-negative roots, the smaller of 
which (vanishing at ri = 0) corresponds to a maximum of 
T^ASF- It is the larger root, existing for a range of positive 
Ti , that corresponds to the minimum of T^asf associated 
with the ASF state. For this range of positive ri the 
transition from the normal state to this ASF minimum 
must therefore be first order. These considerations hold 
for both 7 > 1 and 7 < 1. 

To quantify this, define the cubic discriminant 



in which the parameter 



A = {q/2f + {p/3f 



(Bll) 



t? 



r2 -n 



(B9) 



has been defined. 

The apparent instability of TYasf at large |02| for 7 < 1 
is in fact illusory: The condition (|B6p restricts 02 to a 
finite interval max{0, ^(1 - VI -4ri)} < -02 < ^(1 + 
\/l — 4ri) for any given ri. Similarly, even for 7 > 1, 
(jBSp may have its absolute minimum at negative 02 (e.g., 
if Ti > 0), but only minima at non- negative 0^2 are of 
physical interest. 



where 



P 

q 



1 

7-1 


t2 


- 


1 




2(7- 








4(7-1) 
i2 



7-1 2(7-1)2 



(B12) 



For A > 0, T^ASF ^^^ ^^^y '-'^'^ ^^'"^^ root, while for A < 
it has three real roots. "12^ As A ^ 0~ two of the roots 
merge, and subsequently annihilate as A changes sign. 
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Order parameter w 




Order parameter \^ 

FIG. 16: Graphical illustration of the solution to HB7|I for 
7 > 1 and the cases i2 > (upper panel) and t2 < (lower 
panel). The thicker lines correspond to ri = 0. Curves above 
(below) these correspond to ri > (ri < 0). For f2 > and 
ri < there is a unique positive root (stars) that approaches 
the origin as ri — > 0~ , and then crosses to unphysical negative 
values. As discussed in the text, this signifies a second order 
N-ASF transition at r\ — 0. For t2 < there is a range of 
negative r\ for which there are two positive roots. However, 
the one (triangles) that vanishes at ri = is a local maximum 
of the free energy HB8|) . and therefore unstable. The physical 
root (stars) is the larger one, which never reaches the origin, 
but annihilates with the unstable one at the positive value 
at which the cubic discriminant A, Eq. (IB11|I . vanishes. The 
transition must therefore be first order, taking place at some 
intermediate point where the free energy itself vanishes, but A 
is still positive. For 7 < 1 the result is qualitatively identical. 
For 7 < 1 the far left root moves to the far right, but remains 
unphysical since it corresponds to a local maximum in the 
free energy. For 7 = 1 it disappears entirely. In either case, 
the evolution of the other two roots remains as described. 



There are two branches to the A = curve (shown as 
dashed hnes in Figs. [T71 and fTB]) given by 



3/2 



ri±(t2) 



±[1-1(7-1)^2] ^- 1 + 2(7- l)t2 
2(7-1)2 



(B13) 

of which r'i_|_ corresponds to the merging and subsequent 
disappearance of the two roots of interest. The two 
branches meet and terminate at a cusp (denoted by a 
star in Figs. fTTIandfTSl) at the point 



^2A 



4(7-1) 



riA = 



1 , _ 7(27-1) 

4(7-1)2 '^'^ (7-1)^' 

(B14) 
where all three roots coincide. The first order transition 
must therefore take place at a point ric{t2) in the interval 
< Tic < ri+ where the energy density (|B8|) itself van- 
ishes. For small ^2 one finds ri+ = j^ijl + *^[(7 ~ 1)^2]} 
(the leading quadratic form being exact to all orders 



for 7=1). and the very small deviation ri+ — ric 



at which the transition 



The point ri — 0, r2 ~ 5, 
line turns from second order to first is a tricritical point, 
labeled Can in Figs. [T7] and [THl Numerical results for 
Tie are shown using 7 = 2 and 7 = 0.2, respectively. 



First order MSF-ASF transition: 



Along the putative second order line (jB10[) . the ex- 
tremum condition (jB7p may be factored in the form 



WaSF = 2(7 - 1)(^2 - V'20)(^ - V'2+)(V' - ^-2-) (B15) 



in which '1IJ20 — VP2|77 i^ the MSP state value, and the 
other two roots are given by 



V-: 



2± 



4(7-1) 



3-2(7-1)^: 



20 



(B16) 



± V [2(7 - 1)^20 + 1]2 + 87 



The argument of the square root is positive, so all three 
roots are real, and a picture similar to the lower panel of 
Fig. [TBI obtains. 

a. The case 7 > 1 ; In order that -020 be the physical 
root, the condition 7/120 > 4'2+ is required, leading to the 
inequality 



V'20 > 



1 



2(V7+1)' 



which remains well defined for 7 = 1. Thus, for 



r2 < r2T 



7 



4(1 + V7)2 



(B17) 



(BIS) 



the transition is indeed second order. For r2 > r2T, 
■02+ becomes the physical root: the transition turns first 
order, and takes place when the two energy densities 
match: Hasp = ^msf = -^2/27- The point r2 = r2T, 
Ti = riT = {2^ + l)/4(.y7 + 1)2 is a second tricritical 
point, ■^'^ labeled Cam in Fig. [TTl At r2T, (|B15p has a co- 
incident pair of roots, and the line ri+ and the transition 
line must therefore osculate — see Fig. [T7l 

b. The case 7 < 1; For 7 < 1, 02o must be the 
intermediate root: 1/2- < ^20 5: 02+, leading to the 
inequality 



1 7 1 

< 020 < 



Thus, for 
r2T- = 



2(1 + V7) 



4(1 - V7)^ 



2(1 -V7)' 



(B19) 



<r2 < r2T+ = — 



4(1 + V7)2 

(B20) 
the transition is indeed second order. For r2 > ?'2T+, 
02+ becomes the physical root, while for r2 < r2T-, 0'2- 
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Mean field phase diagram: global {y= 2) 



Mean field phase diagram: detail (y = 2) 




0.1 



-0.1 



-0.2 
0.05 





'X^ 


MSF 


ASF 


^Nt^AM 




^ 



0.1 



0.15 



0.2 



FIG. 17: Mean-field phase diagram for 7 > 1, in the scaled coordinates defined by (|B1|) . with second order phase boundaries 
given by the thinner solid lines, and first order phase boundaries the thicker solid line. The right plot shows an expanded 
detail near the labeled points E and Cam- The value 7 = 2 is used in the numerical computation here, but the basic structure 
remains unchanged for other values. The second order N-ASF line along the positive r2 axis encounters a tricritical point Can 
at r2 = 1/2, below which the transition turns first order and the line bends into positive ri. The second order N-MSF transition 
line along the positive ri axis terminates on this first order line at a critical endpoint E. The first order line continues below 
E, now separating the ASF and MSF phases, but turns second order at another tricritical point Cam, and follows the line 
defined by (|B10|) . For large ri,r2 this line asymptotes to r2 = 7ri, which agrees with (|4.11|l in unsealed units. Its unphysical 
continuation toward the origin is shown by the dash-dotted line. The dashed line is the curve A = 0, and therefore provides a 
bound on the ASF phase boundary. All three lines osculate at the point Cam- The N-ASF and dashed lines also osculate at 
Can- For 7 = the cusp together with the left hand branch of the A = line is pushed off to infinity, but Can, E and Cam 
remain finite and well defined. 



becomes the physical root: in either case the transition 
turns first order, and takes place when the ASF and MSF 
energies match. The point r2 = ?'2T±, ri = riT± = 
[2(^± 1) - l]/4{y/j ± 1)^ are both tricri tical points, 
labeled Caai± in Fig. [THl At both points (JBISP has a 
coincident pair of roots, and at both points the (two dif- 
ferent branches of the) A = curve and the transition 
curve must therefore osculate, as illustrated in Fig. [TH] 
Notice that for 7^1—, Caai+ remains well defined, and 
coincides with Caai, while Cam- is pushed to infinity. 
The semi-infinite second order line is therefore recovered 
in this limit, consistent with the 7 > 1 result. 



5. Critical endpoint for the N— MSF transition: 

Since Hmsf ^ as ra ^ 0, the N-ASF and MSF- ASF 
first order transition curves must meet at r2 = 0. The 
second order N-MSF curve also terminates at this point, 
which therefore represents a critical endpoint, labeled E 
in Figs, [m and [m 

Figures [T71and[T51show complete phase diagrams, com- 
puted for the cases 7 = 2 and 7 = 0.2, respectively, 
and correspond to the two phase diagram topologies il- 
lustrated in Figs. [TOl and [TT] of the main body of the pa- 
per. Here the first order transition curves are computed 
numerically. The osculation with the A = line at the 
various tricritical points is also shown. 



APPENDIX C: THE TWO-BODY MOLECULAR 
BINDING PROBLEM 

The center of mass Schrodinger equation for two par- 
ticles of mass ruA interacting via an attractive potential 
v(r) with microscopic range do is given by 



mA 



-V^V'W + v{r)^P{r) = E^j{r). 



(CI) 



Of interest is the regime where the particle is very weakly 
bound, with binding energy E < obeying \E\ <^ \v\- 
The wavefunction ^ of such a weakly bound state will 
extend a distance much greater than do outside the po- 
tential. A good approximation is then to treat ip ^ ipi^O) 
as essentially constant over the potential region r < dp 
(this notion may be made rigorous using effective s-wave 
scattering parameters). A Fourier analysis of (jCip then 
yields 



^(k) 



mA\vo\ v{k) 

h? fc2 -1- k2 



V'(O) 



(C2) 



where u(k) = wo^(k) is the Fourier transform of the po- 
tential, uq = w(0) < is the area under the potential, 
and K^ — mA\E\/li^. The self-consistency condition de- 
termining E is therefore 



mA\vo\ 



d'^k w(k) 

{2'kY fc2 + k2 



(C3) 
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Mean field phase diagram: global {j= 0.2) 






1 
1 
t 
/ 
1 
/ 
t 


0.5 


ASF 


C.N N 







^*i^E p_ 








^ ' 


^AM- 




t^r^^'^'^^"'^ 


MSF 



Mean field phase diagram; detail (y= 0.2) 



-0.5 



0.5 




FIG. 18: Mean field phase diagram for 7 < 1, in the scaled coordinates defined by (|B1|) . with second order phase boundaries 
given by the thinner solid lines, and first order phase boundaries the thicker solid line. The right plot shows an expanded detail 
near the labeled points E, Cam+ and Cam-- The value 7 — 0.2 is used in the numerical computation here, but the basic 
structure remains unchanged for other values. The upper part of the phase diagram is very similar to that for 7 > 1. The second 
order N-ASF line along the positive r2 axis encounters a tricritical point Can at r2 = 1/2, below which the transition turns 
first order and the line bends into positive ri. The second order N-MSF transition line along the positive ri axis terminates on 
this first order line at a critical endpoint E. The first order line continues below E, now separating the ASF and MSF phases, 
but turns second order at another tricritical point, now labeled Cam+, and follows the line defined by HB10|I . However, now 
the second order ASF-MSF line turns first order again at a new tricritical point Cam-- The unphysical continuation of this 
line above Cam+ and below Cam- is shown by the dash-dotted line. When 7 — * 1~, the latter is pushed out to infinity, while 
the former remains finite. For large ri,r2 the first order line asymptotes to r2 = ^/^ri, which agrees with (|4.14p in unsealed 
units. The dashed line is the curve A = 0, and therefore provides a bound on the ASF phase boundary. All three lines (though 
two different branches of the A = line) osculate at the points Cam+ and Cam-- The N-ASF and dashed lines also osculate 
at Can- 



For d < 2 the right hand side diverges at £■ = 0. There- 
fore a bound state solution exists for arbitrarily weak 
potential. For d > 2 there is a critical potential strength 
|wo,c|, below which the potential fails to bind, given by 

and one may write 

n^r 1 1 

ruA 



m,, 



yo 



d'^k 



C(k) 



{2t:Y k^k^ + K^) 



(C5) 



Close to the critical point the integral is dominated by 
the small k region, k <C vr/a. One may then ignore the 
fc-dependence of v (the integral remains convergent for 
2 < d < 4), to obtain 



mACd \\vo 



1 



1 

N 



where 



Cd^ 



d'^i 



{2ttYu^{u'^ + 1) 



(C6) 



(C7) 



in particular C3 — l/Air. One therefore obtains the power 
law relationship \E\ ^ {\vo\ — |i'o,c|)^^^'^~^'' describing the 
vanishing of the binding energy upon approach to the 
critical point. 



APPENDIX D: THERMAL AND BOSE 

CONDENSATE DENSITY PROFILES IN A 

HARMONIC TRAP 



The most direct probe of a trapped degenerate atomic 
gas is through its spatial density profile, obtained from a 
freely expanding cloud. In this Appendix, the details of 
the density profile calculations for a trapped Bose gas are 
presented. For the weakly interacting case, the profile is 
well approximated by the noninteracting expression 



n{r) = J2 \Mr) 



(Dl) 



n=0 



with the occupation number nn, for a Bose gas given by 
the Bose-Einstein distribution 



n„ = (alfln) = 



1 



=/3(e„-M) _ 1 ■ 



(D2) 



The single-particle spectrum £„ and normalized wave- 
functions (/'n(r) are solutions of the single particle 
Schrodinger equation h(j)n = [— (ft^/2r7i)V^ -|- V{r))(j)n — 
£n(t>n appropriate to a trapping potential ^(r). For a 3d 



harmonic potential ^(r) 



^muiQ r 



jhujQ (that for 
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simplicity is taken to be isotropic) 

£„ — hajQ^n^ + TT-y + n^) = hajQn, (D3) 

</.„(r) = n ( \„. ) i/„^(f,Oe-^-/^ (D4) 

■i^x.y.z 

with Hn{x) the nth Hermite polynomial, a function of the 
normalized coordinates f,j = r^/ro (« = x,y,z), expressed 
in units of the quantum oscillator length vq = yJh/mujQ. 
The chemical potential, /i < (throughout) is deter- 
mined by the total atom constraint 



/oo 
dvniv) = ^ 
n=0 



1 



g/3(e„-Ai) _ I ■ 



(D5) 



In this noninteracting limit, at T = all atoms go 
into the lowest single-particle state (po, forming a Bose- 
Einstein condensate, with 



,r=o 



(r) 



A^ 



7^3/2 j.j^3 



2 / 2 



(D6) 



At finite T, a fraction of atoms is thermally excited to 
higher single-particle states, and 



n(r) 



1 



2/3(;,-p) _ I ' 



2_^e'''^''[r\e '-'""ir) 



pPt^p/ 



~Pph\ 



OO 



^e^'^Ppo.c(r,r;p/3) 
p=i 



(D7) 



is expressible purely in terms of diagonal elements of the 
single-particle density matrix for a harmonic oscillator 



Posc(r,r',/3) = (rle-'^^lr'). 



(D8) 



governed by a single-particle Hamiltonian h. The density 
matrix can be found by solving a diffusion equation in a 
harmonic potential (or equivalently obtained from an- 
alytic continuation of the harmonic oscillator evolution 
operator) with the "initial" condition of posc(r,r';/? = 
0) = (5(r - r') obvious from (IDS]) .^"^ For d = 3 one has 



Posc(r,r;/3) = 
where 

ro'(/3) 



muji^e 



0hujo 



27r?isinh(/3?iaJo) 



3/2 



-rVra^iP) 



(D9) 




coth{l3hujo/2) 



huJo/kBT > 1 
hujo/kBT < 1, 



(DIO) 



is the finite-temperature "osci llator le ngth" that reduces 
to the quantum one tq = ^hjrauiQ at low T and the 



classical (thermal) one, tt = \j2kBT jmuji^ (defined by 
\mJlr\ = kBT) at high T. 

The p sum in n(r)^ Eq. (jDVp . can be evaluated analyt- 
ically in various limits. For hw^jkBT <^ 1 [valid for all 
but extremely low T, where (jD6p holds] the p sum natu- 
rally breaks up into two parts with n{r) = nrir) + no{r). 
The two contributions correspond, respectively, to ranges 
1 < p < Pc = kBT/hujQ ("thermal") and Pc l£ P < oo 
("quantum") with pc determined by PcPhw^ = 1. The 
thermal range is characterized by p such that pfifkuQ < 1, 
within which the finite temperature auxiliary oscilla- 
tor length ro(/3) and the corresponding density matrix 
Posc(r, r;p/3) can be approximated by their thermal clas- 
sical forms, giving 



nrir) 



f ksT 
\2iThwo 



\2'k1VjJq 



where a "cutoff" extended zeta function. 



3/2 


1 

l3 
'^0 


Pc-1 

>; 
p=i 


1 
p3/2 


^-V(r 


V4+IM 


/kBT) 


3/2 


1 

l3 
'^0 


ff3/2 


\ 


V4- 


iMl/fcsT 


kBT\ 
hjjQ ) ' 
(Dll) 



ga{x,Pc) 



'^\P 



p=i 



pu 



(D12) 



has been defined. 

The quantum density contribution ^.^{r) is character- 
ized by values of p such that pPhujQ > 1, and thus 
by a zero-temperature oscillator length ro(p/3) « vq 
and the density matrix is given by posc(r7r;p/3^o) ~ 

2/2 

-r /rg^ The resulting sum is then easily com- 



-3/2^-3^ 



puted, yielding 






jO / -^ ip-J 



with amplitude factor 

No{T) « £ 



p\tl\/kBT 



-{Pa-l)\t^\/kBT 
g\^i\/kBT _ I 



(D13) 



(D14) 



being the number of bosons occupying the single-particle 
ground state, reduced by the factor e-(p<:-i)l'^l/'=B^ r:; 

Eliminating pc = ksT/huJo, it is clear that nxir) and 
no(r) depend strongly on the ratio of the chemical po- 
tential /i to the trap level spacing fkoo, with the former 
determined by the temperature through the total atom 
number constraint. For |/z|/Swo ^ 1 (which corresponds 
to T > Tc), the sum in the expression for riT^r), Eq. 
(jDlip . can be extended to infinity, introducing only ex- 
ponentially small error 0{e~'^^^^°) <C 1. The resulting 
thermal density nxir) is then given by the extended zeta 
function 



nrir) 



(. 



knT 



V 2Tlhwr 



3/2 



r33/2 



,-rV4-|Ml/fcBT 



for T > Tc, (D15) 



As expected for T > Tc, the condensate spatial distribu- 
tion is still given by the Gaussian expression (|D13[) . but 
with an exponentially small condensate 



NoiT) 



-\fj.\/hujo 



0, for T > Tc 



(D16) 



At high T :$> Tc = huJoiN/C{3)y/^ (where 
the gas is nondegenerate) , such that < |/i| ~ 

-kBTln[(^'\ N] w 3fcsrin(T/rj > fcsT, the ther- 
mal (and therefore total) density (|D15p reduces to a pure 
Gaussian with thermal width r^, reflecting the high- 
temperature Boltzmann statistics 



n{r) 



2'KhujQ 

N 



3/2 



'O 



7r3/24 



'■t, for r>Tc 



(D17) 
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in which the relation TV = {ksT /tiuJo)'^e^'^^/^^'^ has been 
used to satisfy the particle number constraint. 



As T is lowered, approaching Tc from above, the mag- 
nitude of the chemical potential drops below T (remain- 
ing negative), and (while the condensate fraction remains 
vanishingly small) the boson density profile (IDISP devel- 
ops a non-Boltzmann peak structure even above T^: 



n{r) 



keT 
2'Khujo 



3/2 



2 / 2 



-iMl/fcsT 
27 

r 



^C(3/2)-2.V2(^ + J|L) 



1/2 






(D18) 



retaining a Gaussian falloff at large r/rx- The small- 
r cusp in n(r), that develops as Tc is approached from 
above, is rounded on the cutoff length rdT), which for 
T > Tc is given by rdT) w r^ = rQ^y2\^\/hua, with 
ro < rc{T) < rr- 

From the prefactor in n(r) it is clear that this thermal 
form cannot persist to low temperatures, as the volume 
under n(r) drops with T. Thus, for even lower tem- 
perature, T < Tc ~ fuvN-^/^, to accommodate all N 
particles |/i| is forced to drop below the level spacing, 
ImI ^ f^Q- In this region, the high p terms in (|D11[) . 
estimated to be 0(e-P-l''l/'=«^) = 0{e-\''\^'^°), are no 
longer exponentially small and the sum can no longer, in 
general, be extended to infinity. However, for pc ^ 1, 
on length scales longer than the cutoff length Vc (T) , the 
thermal density nxir), Eq. (JD11[) . is still well approxi- 
mated by the extended zeta function in Eq. (jD15|l . but 
with ji/UbT ~ and its small r cusp smoothed out on a 
length scale rc{T < Tc) ~ ttI \fpc = 2ro [following from 
the condition Pcir/rxY — 0{\)\. 

Correspondingly, the condensate contribution no('') 
begins to grow for T < Tc through the growth of the coef- 
ficient No{T), Eq.lml In fact for T < T^, 1^1/^5^ < 1 
and \^i\/fuvo < 1, one has No{T < Tc) w ksT/l^i] > 1, 
and the finite fraction of the bosons condensed into the 
lowest single-particle oscillator Gaussian state (with a 
narrow width ro) are precisely what is required to make 
up for those that cannot fit into the thermal distribu- 
tion riTir). Thus for T < Tc the total atom density 
profile n{r) changes dramatically, developing an easily 
identifiable bimodal distribution n{r) = nxir) + no(r). 



illustrated in Fig. [T] 

As discussed in Sec. IV Bl this analysis easily general- 
izes to the two-component Bosc gas (bosonic atoms and 
molecules) that is the focus of this paper. 



APPENDIX E: DETAILS OF BOSE-BCS 
MSF-ASF CRITICALITY 



In this appendix, the critical behavior in the vicinity 
of the ASF transition line is derived as a function of di- 
mension d via an analysis of the integral equations (|8.43p 
for small deviations p,T, f . 

For 3 < d < 4 one obtains finite values for the deriva- 
tives 



Ji = -a^J(0,0) 



cf^i 



(27r)'^[u2(u2 + 2)]3/2 
dfiu 1 



J2 = 9^J(0,0) = - , ,„ .^. 2. 2 

2 J (27r)'^ [w^yu^ 

Rl = -drRiO,0) = J2/Id 

i?2 = (9/i?(0,0) = Ji/4/d. 



2)]3/2 



(El) 



Higher order derivatives, however, lead to divergent in- 
tegrals at small u. The subtracted integral has leading 
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behavior 

SJ{rJ) 

Jc(s) 

where 



J{r,f) + JiT-J2f^jc{s) 
s 
2^^~ 



1 f (fu fl 



V2j i^7r)d 



(l-r)2_(l-/)^/-2T 



2(1 -r) 



v^ 



>o. 



(E2) 



(E3) 



With the subtractions, the integral jc(l) converges at 
both large and small u. The vanishing of J requires now 



= J2/- ./lT+jc(l)s 

and leads to 



f 



Jj_ Jc(l) 

J2^ J 2 



2Jo 



- 1 



+ 0(/2,r2), 



+ 0{t' 



(E4) 



(E5) 



It is easy to see from (IE1[) that Ji > 2J2 and that this 
solution is consistent with the requirement that s > 0. 
Similarly, for the density deviation one obtains 



P 



R2I - RiT 



2/rf 



Jl 



J2 

Id V±J2 

2 



^'w^-''^- 



Jc(l) 



2/d V2J2 



o{f- 

Ji 



- 1 



+ O(r^). 



(E6) 



For d < 2) integrals in (jEip diverge at small m indicating 
singular dependence on r and /. This infrared singularity 
has been isolated by writing 



J{rJ) - Jc(t,/) + <5J(t,/) 
R{T,f) = i?e(r,/) + <5i?(T,/) 



(E7) 



where 



Jc{tJ) = 

Rc{r,f) - 

jcis) = 



VI- 
Vi- 


=Jc[s) 

T 

'''1 ('«'! 


2/rf 

' 1 


Jc{S) 



V27 (27r)" 



V^ 



(E8) 



where jc(s) now requires only the single subtraction for 
convergence at large u in this lower dimension. It can 
be checked that for 1 < d < 3 the derivatives of the 
subtracted integrals 5 J and 5R are finite at the critical 
point, and so have leading linear dependence on t and /. 
The corresponding coefficients are 



Ji 



J2 
Ri 



(2^ 
1 /■ d'^u 
2 



1 



1 



[7/2(^2 + 2)]3/2 2^27/3 



1 



{2ttY i[u2(M2+2)]3/2 

h/h, R2 = Ji/4/d. 



2V2u3 



(E9) 



One obtains again 



3c{s) = 3c{l)s~ 



(ElO) 



in which jc(l) is finite for 1 < rf < 3 (note that d = 1 
is the lower critical dimension, below which the phase 
transition ceases to exist, and so one expects special be- 
havior here). This dominates the linear behavior, and 
the constraint J{t, /) = becomes 

= ~ic{^)s^ - JlT + J2I + 0{t\ /2), (Ell) 

which, to leading order, has the solution s = 0, and hence 



/ = 2r + 2 



Ji - 2 J2 



0{t^ 



(E12) 



From ((E9| it is easy to see that Ji > 2J2, so that (|E12p 
is again consistent with s > 0. Using (|E7|1 and (jESp . the 
density deviation is 



p = R2f-RlT 

Ji - 2 J2 
2/rf 



2id 



2t- 



Ji - 2 J2 

Jc(l) 



7- d-1 




The singular corrections to the linear behavior may be 
identified with the energy exponent 1 — a, where the 
Gaussian specific heat exponent is a = |5f (not to be 
confused with the Feshbach resonance coupling) . In a full 
theory this form would be replaced by the exact Ising ex- 
ponent. 

At the upper critical dimension d — 3 itself there will 
be logarithmic corrections. Since d = 3 is the physically 
important dimension, it is worth presenting the results 
for this case as well. One may now express the results in 
terms of elliptic integrals (see Ref.lisl pp. 232, 235, 905): 
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(E15) 



At small T, /, hence g — > 1, using the asymptotics of 
the elliptic integrals (see Ref. H^, p- 906), the J — 
constraint takes the form 

/ 64e~ 
+ 0[t^ ln(/ - 2r), f ln(/ - 2r)], (E16) 



with solution 
/ = 2r(l- 



ln(8/e3r) 
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(E17) 



The density deviation takes the form 
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